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CAPITOLO 1 



Introduzione 

1. Il problema di una descrizione microscopica 

1.1. Equilibrio e non equilibrio: definizioni. In questa tesi vogliamo discutere il problema della 
costruzione di una meccanica statistica dei sistemi fuori dall'equilibrio, ovvero di una teoria che, partendo 
dalla dinamica microscopica, permetta di predire relazioni tra valori medi di osservabili fisiche in uno stato 
di non equilibrio. 

Innanzitutto, dobbiamo specificare cosa intendiamo per equilibrio e non equilibrio. Un sistema all'e- 
quilibrio è un sistema di particelle in uno stato stazionario che evolve risolvendo le equazioni del moto di 
Hamilton; chiamiamo stato stazionario uno stato fisico che non muta nel tempo, ovvero se la termodinamica 
del sistema è descritta da un insieme di osservabili A — {Ai} i valori medi delle Ai sono costanti 

Un sistema fuori dall'equilibrio può essere tale per due motivi: 

(1) le particelle evolvono risolvendo equazioni del moto che non sono le equazioni di Hamilton, e/o 

(2) il sistema non si trova in uno stato stazionario. 

In seguito, considereremo esclusivamente stati stazionari fuori dall'equilibrio: in generale, il nostro siste- 
ma sarà soggetto a delle forze non conservative - ovvero che non possono essere espresse come il gradiente di 
un potenziale - e, per permettere l'esistenza di uno stato stazionario, il calore prodotto dovrà essere assorbito 
da un termostato. 

1.2. Caratteristiche di una teoria microscopica. E chiaro che, all'equilibrio o fuori dall'equilibrio, 
lo studio di 10^"^ equazioni del moto è un problema intrattabile. Ad ogni modo, è ragionevole pensare 
che ai fini del calcolo di "poche" osservabili fisiche i dettagli della traiettoria della singola particella siano 
inessenziali; prima di scrivere una ad una le equazioni del moto su un calcolatore e aspettare il tempo 
necessario affinché vengano risolte, dobbiamo chiederci se sia davvero necessario avere così tanta informazione. 
Quali sono le caratteristiche del moto che hanno una rilevanza da un punto di vista macroscopico, ovvero 
che giustificano la termodinamica? 

Si potrebbe provare a rispondere a questa domanda assumendo le proprietà del moto che riteniamo 
plausibili, e verificando l'accordo tra le nostre predizioni teoriche e la realtà fisica. Questo è esattamente ciò 
che accade nello studio dei sistemi all'equilibrio, dove l'assunzione sulla dinamica microscopica che permette 
di giustificare la termodinamica prende il nome di ipotesi ergodica; nonostante nelle applicazioni questa 
ipotesi sia praticamente impossibile da dimostrare, le sue conseguenze fisiche sono largamente verificate. 
Ciò suggerisce che non abbia senso, o perlomeno che sia superfluo, verificare la stretta validità dell'ipotesi 
ergodica ogni volta che si vogliono usare le tecniche della meccanica statistica dell'equilibrio; piuttosto, 
l'atteggiamento corretto è quello di assumere che al fine di calcolare quantità macroscopiche un sistema 
all'equilibrio si comporti come se fosse ergodico. 

I sistemi fuori dall'equilibrio sono incompatibili con l'ipotesi ergodica. Per derivare la termodinamica 
partendo dalle leggi del moto è necessario fare delle nuove assunzioni sulla dinamica microscopica, e verificare 
che la teoria che si ottiene sia consistente con i risultati sperimentali e con quanto già noto all'equilibrio; 



Ad esempio, possiamo immaginare che A- sia formato dall'energia potenziale media <I>, dall'energia cinetica media T, 
dall'energia totale U, dalla pressione P, dalla densità p e dal volume V. 
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potremmo dire di aver raggiunto l'obiettivo piìi ambizioso se riuscissimo a costruire una teoria che permetta 
di descrivere la statistica di un sistema allo stato stazionario indipendentemente dal fatto che si trovi o meno 
all'equilibrio, e che dunque si riduca all'usuale meccanica statistica dell'equilibrio nel limite di forze non 
conservative nulle. 

Negli ultimi 70 anni molte proprietà di sistemi "quasi" all'equilibrio, ovvero sottoposti a forze non con- 
servative con "piccole" intensità, sono state efficacemente descritte con le tecniche della meccanica statistica 
dell'equilibricj^ chiaramente, c'è un'ambiguità in cosa si intende per "piccolo" , e questa ambiguità non può 
che essere risolta caso per caso. In queste situazioni si assume che tra equilibrio e non equilibrio "le cose 
cambino poco" in prima approssimazione, e con questo approccio si riescono a dimostrare risultati notevoli 
come la formula di Green - Kubo, meglio noto come teorema di fluttuazione - dissipazione, e le relazioni di 
reciprocità di Onsager (assumendo che la dinamica sia reversibile); parlando di questi risultati ci riferiremo 
alla teoria della risposta lineare. 

Quindi, un'altra caratteristica richiesta ad una teoria microscopica del non equilibrio è che per "pic- 
cole" forze non conservative implichi i ben noti risultati della teoria della risposta lineare, possibilmente 
estendendoli agli ordini successivi. 

2. Un esperimento 

Lo studio dei sistemi fuori dall'equilibrio ricevette un forte impulso nel 1993, grazie ai risultati speri- 
mentali di D. J. Evans, E. G. D. Gohen, G. P. Morriss, [2]; gli autori simularono un fluido bidimensionale 
composto da iV = 56 particelle, che evolve con le equazioni del moto 

(1) = + ny, , = Zj - iiPy, - ap^ , 

du (q ) 

dove p . = {p^.^py.^, 7 = "f{q .) = — }:^'' se u{q .) = {ux{q .), Wj,((7 .)) è la velocità della particella j-esima 

3 3 ^ 3 3 3 

nel punto g^., iè il versore nella direzione x, a è determinato con il principio di minimo sforzo di Gauss {cf. 
capitolo [3|, e è la risultante delle forze d'interazione tra la particella j e le rimanenti — 1 = 55. Le 
particelle interagiscono mediante un potenziale 
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per r < 26(7 
altrimenti. 



Con questo modello Evans, Cohen e Morriss verificarono numericamente che, chiamando a il tasso di 
produzione di entropicj^ a-r la sua media su un tempo r e indicando con Freq{B) la frequenza delle occorrenze 
dell'evento B, 



1 1 Freq{ar = P) . 

a meno di correzioni tanto piìi piccole quanto piià è grande r; l'evento {ar = —P} può essere visto come 
una "violazione" della seconda legge della termodinamica, vera solo nel limite t — >■ oo. Nello stesso lavoro 
Evans, Cohen e Morriss provarono a giustificare la ([S]) congetturando che la probabilità che il sistema si trovi 
a percorrere in un tempo r il "segmento di traiettoria nello spazio delle fasi" i abbia la forma 



■^Vedere [T] per una esauriente esposizione di questi risultati. 

''Come discuteremo nel capitolo [s] non è ovvio chi siano l'entropia e il suo tasso di produzione in un sistema fuori 
dall'equilibrio; ad ogni modo, nel caso di termostati gaussiani il tasso di produzione di entropia può essere identificato con il 
tasso di contrazione dello spazio delle fasi, proporzionale alla quantità a che appare nella 



dove è Tn-esimo esponente di Lyapunov positivo calcolato sul segmento i-esimo, e Z è un'opportuna 
normalizzazione. 

Come abbiamo accennato poco fa, il parametro a che appare nella ([T]) è determinato con il principio di 
Gauss; nella sezione |3.1| del capitolo [3] vedremo che questa scelta renderà reversibili le equazioni del moto, 
e inoltre ci permetterà di stabilire un'identificazione tra il tasso di produzione di entropia e la somma di 
tutti gli esponenti di Lyapuno\j^ Grazie a questa corrispondenza e alla reversibilità del moto, assumendo la 
Q come misura di probabilità invariante Evans, Cohen e Morriss riuscirono a dedurre la ([3|, e conclusero 
interpretando il risultato sperimentale (|3| come una verifica della loro congettura Q. 

3. Il teorema di fluttuazione di Gallavotti-Cohen 

L'argomento di Evans, Cohen e Morriss, però, è troppo vago per essere considerato una dimostrazione del 
risultato (|3]); data l'importanza della questione, è necessario capire perché la misura di probabilità invariante 
di un sistema fuori dall'equilibrio dovrebbe avere la forma Q. 

Nel 1995 G. Gallavotti e E. G. D. Cohen proposero, c/. [3], [4], un punto di vista che, se accettato, 
permette di dire con rigore matematico quale sia la misura di probabilità invariante di un sistema fuori dal- 
l'equilibrio; in questo contesto la formula ([s]) diventa un teorema, [5, , ,5^, noto come teorema di fluttuazione 
di Gallavotti - Cohen. La proposta di Gallavotti e Cohen consiste in un'assunzione sulla dinamica delle 
particelle che generalizza l'ipotesi ergodica, perché implica che all'equilibrio lo stato del sistema sia descritto 
dalla misura di Lebesgue sulla superficie dello spazio delle fasi ad energia costante, ovvero dall'ensemble 
microcanonico; l'assunzione consiste nel dire che allo scopo di misurare proprietà macroscopiche un sistema 
macroscopico può essere pensato appartenente ad una speciale classe di sistemi dinamici caotici, chiamati 
iperbolici transitivi o Anosov, per i quali il valor medio di un'osservabile O nel futuro può essere scritto 
informalmente comejf] 



(5) (O), = 



^traiettorie nello spazio delle fasi L V Jì 



z 

dove x{t) è lo stato del sistema lungo la traiettoria nello spazio delle fasi al tempo t, A„ è il logaritmo 
del coefficiente di espansione^ ovvero la somma di tutti gli esponenti di Lyapunov positivi del sistema di 
particelle, e {x{tk)} è un "campionamento opportuno" dell'evoluzione del sistemgj^i — >■ x{t). Sotto le stesse 
ipotesi che implicano la ([s]), invertendo la freccia del tempo il valor medio di un'osservabile nel passato 
può essere ottenuto sostituendo —Au[x{t)] con il logaritmo del coefficiente di contrazione As[x{t)]. Infine, 
dal momento che la misura è invariante possiamo sostituire l'argomento x{0) dell'osservabile O con x (f) e 
t' generico. L'assunzione che implica la ([s]) è chiamata ipotesi caotica, e verrà discussa ampiamente nella 
sezione [2] del capitolo [3j 



4. Organizzazione del lavoro di tesi 

Il lavoro di tesi è organizzato in 3 capitoli e 4 appendici; ne riassumiamo brevemente il contenuto. 



Più precisamente, potremo identificare il tasso di produzione di entropia con il tasso di contrazione dello spazio delle fasi. 

^Per dare senso alla ^ è necessaria una discretizzazione dello spazio delle fasi; come vedremo, per la classe di sistemi 
"molto caotici" che studieremo ciò potrà essere fatto in modo relativamente facile. 

^Se la dinamica avviene ad intervalli di tempo discreti, come assumeremo sempre in questa tesi, non avremo bisogno del 
campionamento. 



Nel capitolo [2] presentiamo qualche concetto della teoria matematica che ci permetterà di com- 
prendere l'ipotesi caotica, con particolare attenzione ad una semplice applicazione nella quale è 
possibile calcolare esplicitamente i potenziali della misura di probabilità invariante. 
Nel capitolo [s] riprendiamo il problema del non equilibrio accennato in questa introduzione, in 
particolare discutendo l'identificazione tra tasso di produzione di entropia e tasso di contrazione 
dello spazio delle fasi, e presentiamo una dimostrazione rigorosa del teorema di fluttuazione di 
Gallavotti-Cohen. In seguito mostriamo una generalizzazione del teorema, e riassumiamo qualche 
altro risultato noto fuori dall'equilibrio. 

Nel capitolo |4] discutiamo alcune conseguenze del teorema di fluttuazione; il lavoro originale 
consiste in una generalizzazione della teoria della risposta lineare, ricavando a partire dal teorema 
di fluttuazione alcune relazioni tra funzioni di correlazione a tutti gii ordini. In particolare, saremo 
in grado di scrivere il valore di aspettazione di qualunque osservabile O dispari sotto inversione 
temporale attraverso uno sviluppo in funzioni di correlazione O — cr, se cr è il tasso di produzione di 
entropia; come verificheremo, all'ordine piìi basso questo risultato implicherà la formula di Green 
- Kubo e le relazioni di reciprocità di Onsager. In seguito, applicheremo le relazioni ottenute per 
verificare il teorema di fluttuazione in un semplice caso in cui, da alcune simulazioni numeriche, 
sembra essere vero nonostante non ne siano verificate le ipotesi. 

Nell'appendice [A] discutiamo una derivazione euristica della misura di probabilità invariante che 
descrive i sistemi fuori dall'equilibrio nel contesto dell'ipotesi caotica. 

Nell'appendice [B] calcoliamo esplicitamente il funzionale di grandi deviazioni del tasso di produ- 
zione di entropia mediato su un tempo finito. Il contenuto di questa appendice è originale. 
Nell'appendice [C] discutiamo una tecnica matematica che ci permetterà di studiare le proprietà 
di analiticità della misura di probabilità invariante che descrive i sistemi fuori dall'equilibrio nel 
contesto dell'ipotesi caotica. 

Infine, grazie alla tecnica esposta nell'appendice [C] nell'appendice [D] discutiamo le proprietà di 
analiticità dei valori medi di osservabili analitiche. 



CAPITOLO 2 



Cenni di teoria ergodica 



1. Sistemi dinamici astratti 

In questa sezione vogliamo presentare una breve introduzione alla teoria astratta dei sistemi dinamici, e 
ricordare qualche concetto fondamentale di teoria della misura. 

Definizione 1.1. (Sistemi dinamici discreti) 

Sia n uno spazio metrico compatt^ e separabili e sia S una mappa continua di fi in sé stesso. La 
coppia (Q, S) sarà chiamata sistema dinamico topologico discreto su il. (f2, S) sarà detto invertibile se S 
ha un'inversa continua S^^ . 

Sia fi una misura di probabilità compiei^ definita sulla a-algebra dei sottoinsiemi di Vi, dove N E B 
è un insieme di ^-misura nulla, e supponiamo che S sia una mappa misurabile rispetto a B al di fuori di 
ÌV0 Se 



(6) fi{A) ^ ^{S-^A) per ogni A e BnN", 

dove N'^ — Cl/N è il complemento di N, diremo che {fi, S) è ^.-invariante. La tripla {^,S,^) sarà 
chiamata un sistema dinamico discreto mod definito al di fuori di N. Se N = $ la tripla {Q,S,fi) sarà 
detta un sistema dinamico metrico. 

La definizione che abbiamo appena dato coinvolge sistemi dinamici che evolvono ad intervalli di tempo 
discreti; l'analogo nel continuo della definizione 1 1 . 1 1 può essere formulato nel seguente modo. 

Definizione 1.2. (Flussi) 

Sia fi uno spazio metrico compatto e separabile. Sia {St)t<£R o {St)t£R^ un gruppo o un semigruppo, 
omomorfo a M o a di mappe che agiscono con continuità su fl^ 

Le coppie {fi, (S't)tgE), {fi, (»S't)tgR^) saranno chiamate rispettivamente un flusso topologico invertibile o 
un flusso topologico su fi. 

Sia /i una misura di probabilità completa su una a-algebra B e sia {St)teK. un gruppo di mappe /i- 
misurabili omomorfe a M e che conservano fi. Supponiamo che la funzione {t,x) — )■ SfX definita in R x fi 



Un insieme compatto G £ H è un insieme chiuso e limitato; un insieme chiuso è un insieme che contiene i limiti delle sue 
successioni, mentre un insieme è detto limitato se può essere contenuto in una "palla" di raggio finito, ovvero se esiste x £ fi 
tale che, indicando con d una distanza su H, per ogni g & G d{x, s) < r, r > 0. 

^Un insieme separabile G è un insieme che contiene un sottoinsieme denso e numerabile; ciò vuol dire che l'insieme G può 
essere " approssimato arbitrariamente bene" da una famiglia numerabile di suoi sottoinsiemi, ovvero che possono essere messi 
in corrispondenza uno-a-uno con l'insieme dei numeri interi Z. 

■^Una misura fi è detta completa se ogni sottoinsieme di un insieme di /i-misura nulla ha ^-misura nulla. 

'^Vn' algebra è una collezione di insiemi contenente l'insieme vuoto chiusa rispetto alle operazioni di complemento e di 
unione finita; se è anche chiusa rispetto all'operazione di unione numerabile, la collezione viene detta una a-algebra. 

^Una funzione S è detta misurabile rispetto alle cr-algebre B, B' se l'immagine inversa S~^{E) di un insieme E a B 
appartiene a S'; in particolare, se B = B' = Bf^, dove è la cr-algebra degli insiemi /j- misurabili, la funzione è detta 
/^-misurabile. 

^Questo significa che la funzione (t, x) — > Stx di R X fi, o di M^. X f2, in f2 è continua. 
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con valori in f2 sia ^-misurabile rispetto alla a-algebra generata dagli insiemi ir{]S(M) X B. Allora il nusso 
(ri, (>5't)tgR, /x) sarà chiamato un flusso metrico invertibile su iì. 

Se sostituiamo (S't)tgR con un semigruppo {St)teR+j '^^(1^+) ^ B-misurabile, otteniamo la nozione di 
flusso metrico. 

Nel seguito ci occuperemo sempre di sistemi dinamici discreti, dal momento che un sistema dinamico 
continuo può essere ricondotto ad un sistema dinamico discreto mediante la tecnica della sezione di Poincaré, 
cf. [6]. Inoltre, da un punto di vista fisico i sistemi dinamici continui sono piuttosto inusuali, perché in un 
esperimento è impossibile controllare la dinamica di un sistema fisico in tempo continuo; quello che si fa, 
invece, è acquisire informazioni sullo stato del sistema a determinati intervalli di tempo, che, a causa delle 
limitazioni sperimentali, non possono essere arbitrariamente piccoli. 

Parlando di sistemi dinamici discreti, qualifiche tipo topologico, metrico ecc. possono essere omesse 
riferendosi semplicemente a sistemi dinamici astratti (fi, S), dove fi è lo spazio sul quale agisce la mappa S. 

Per concludere, introduciamo la nozione di isomorfismo, o coniugazione, tra sistemi dinamici. 

Definizione 1.3. (Isomorfismi) 

S'è {^Ì,S) e (r2',S") sono due sistemi dinamici astratti diremo che sono isomorfi, o coniugati, se esiste 
una mappa invertibile / : il n' tale c/i(0 

(7) IoS = S'oI. 

2. Dinamica simbolica 

L'evoluzione delle traiettorie di un sistema dinamico {Q, S) può essere descritta misurando ad ogni 
instante di tempo prefissato delle quantità di interesse, e compilando una lista con tutte queste caratteristiche 
in funzione del tempo; matematicamente, sintetizziamo queste informazioni nella storia g^{x) di un punto x. 
Ovviamente, un punto x identifica univocamente la sua storia, mentre è naturale chiedersi se da una storia 
sia possibile risalire al punto iniziale; approfondiremo questa questione tra poco. 

Proviamo a formulare il concetto di storia di un punto in modo astratto. Considerando una partizioni 
V = {Po, Pi, Pn} di ri, n > 1, diremo che la sequenza <7{x) = {o'ì{x)}ì^z G {0, n}^ è la storia simbolica 
di un punto x se 

(8) S'x e P^^^) Vi e Z; 

l'insieme Pi, i = l, 2,3, ...,n, è la collezione di punti di il che ha le caratteristiche indicate dall'etichetta 
i, cosicché l'unione Uf^iPi è l'insieme dei punti che hanno almeno una delle caratteristiche i = l,2,...,n 
mentre Pq — ft \ Uf^^Pi è l'insieme dei punti che non ne ha nessuna. Richiederemo sempre che gli elementi 
della partizione V (gli atomi della partizione) siano insiemi boreliani in B{n); partizioni di questo tipo sono 
chiamate partizioni boreliane. 

Possiamo pensare le informazioni acquisite sull'evoluzione del sistema come una successione a e {0, n}'^, 
generata prendendo un grande tempo di osservazione N durante il quale il moto visita in successione i sot- 
toinsiemi P^., i — —N,...,N; la sequenza o;^ € {0, che risulterà dalle osservazioni avrà la 
proprietà 



^Con la notazione S(fì) indichiamo la cr-algebra dei boreliani di Q, ovvero la cr-algebra generata dagli insiemi aperti di Q. 
Se n = R è possibile dimostrare, cf. \7\, che la cr-algebra dei boreliani può essere equivalentemente generata dagli insiemi chiusi 
di R, oppure dagli insiemi semi-aperti di R ecc. 

*^Con il simbolo "o" indichiamo la composizione tra mappe, ovvero {So I){x) = S{I{x)). 

^Una partizione, o un pavimento, V = (Pi , Pn) di H è una collezione di sottoinsiemi di fi tali che [J^Pi = Q e PiDPj =0 
per i 7^ j. 



(9) n^=-N S-'P., + 0, 

dal momento che S^x e P^^, Vfc (e quindi l'intersezione nella ([9| conterrà almeno un punto). 
Per quanto detto, è naturale pensare che un moto, o piuttosto la sua storia, abbia la proprietà di 
intersezione finita, ovvero che 

(10) n^ej S-^ Pa, ^ <è VJcZ,|J|<oo, 

dove I J| indica il numero di elementi di J. Quindi, indentifichcrcmo lo spazio dei moti osservati in V 
con l'insieme delle sequenze infinite 

(11) ^ = {<l\<ie {0,...,n}^,njejS-^Pa, 7^0VJcZ, |J| <oo}. 

Inoltre, diremo che una partizione V è S-separante se la storia a indentifica univocamente il punto x; 
questa proprietà equivale alV espansività di S in V, ovvero al fatto che, per ogni q_ € fij^ 

(12) Jim (diamn^_^ S'^^P^.^) = 0. 
Riassumiano questi concetti nella seguente definizione. 

Definizione 2.1. (Moti simbolici) 

Sia (ri, S) un sistema dinamico topologico invertibile e sia V — {Pq, Pn} una partizione di fi in n + 1, 
n > l, insiemi boreliani. 

(1) Consideriamo l'insieme 

(13) ^ = {q_\a € {0, n}^ nf_^5-JP,^. ^ ViV} ; 

chiamiamo 17 l'insieme delle (V , S)-stoTÌe dei moti simbolici generati da 5 in fi e visti da V. 
Quando non ci sarà ambiguità chiameremo fi semplicemente insieme dei moti simbolici. 
Se X € n chiameremo la (V, S) -storia di x l'elemento a{x) di O tale che 

(14) a;e5-^P,^(,) Vj e Z. 

(2) S'è la relazione (7_{x) = g_{x') implica che x — x' diremo che V è S'-separante. SeV e S verificano 



la (12) diremo che S è P-espansiva o che è espansiva in P. 

(3) La corrispondenza definita dalla mappa E : il — >■ $7 che associa ad ogni x Cz fi la sequenza a;(x) G f2 
sarà chiamata il codice della dinamica simbolica di S rispetto a P . 

(4) Infine, considerando gli insiemi J ^ {jiT--,jq) C Z e /e sequenze (Tj E {0,...,n}'^ useremo la 
notazione 

(15) Pi ^ P?^:X = ^oejS-'Pa, 

per indicare i punti la cui storia in J è specificata da a_j . 



Il diametro di un insieme G C dove fi è uno spazio metrico, è la massima distanza tra duo punti di G. 



Osservazione 2.1. (1) Se iì ^ {0, ...,n}^, S = t è la traslazioncp^ in {0, ...,n}^ e Pi — {(t|(To = i}, 

i = 0, 1, n, gli insiemi saranno indicati anche con e saranno chiamati cilindri di {0, n}^ 
con "base J e specificazione a", ovvero 

(16) C•^ = k'k'€{0,...,n}^a;.^ \fk = l,...,q}, 

se J ^ {ji, ...,jq} eg_^ (ai, (Tg) e {0, 
(2) In generale fi D ^(fì); questo vuol dire che non tutte lo storie sono permesse. 

Se (f2, S) è un sistema dinamico astratto e V è una partizione espansiva possiamo "invertire" la mappa 
codice S, ovvero partendo da una storia a possiamo identificare il punto che la genera; dato un punto (J G il, 
introduciamo l'insieme X[a) G 17 come 

(17) Xia)^n^^_^S-^P.^, 



dove Pj è la chiusura di Pj. L'insieme X{(t) è non vuoto perché Q è compatti 




Proposizione 2.1. (Codice della dinamica simbolica) 

Sia {^,S) un sistema dinamico topologico invertibile; sia V = {Pq^ ...^ P^} una partizione topologica di 
il sulla quale S è espansiva. 

(1) L'insieme X{a) definito nella (j7[ ) contiene uno ed un solo punto cosicché possiamo definire la 
mappa X : Q fi chiamando X{(7) l'unico punto in X{a) oppure come, con un leggero abuso di 
notazione, 



(18) ^^X{a) = nf^_^S-^P,^. 

(2) X è una mappa continua. 

(3) Y^~^X~^(x) = X per ogni x £0,. 

(4) X e Y, sono una l'inversa dell'altra se considerate rispettivamente come mappe tra 17 e Ef7, EÌ7 e 
Vi. 

(5) S{x) — X{tYj{x)), dove {Ta)i — Ui+i è la traslazione di una unità temporale a sinistra della 
sequenza a_. 

Dimostrazione. Vedere [8]. □ 

Intuitivamente, ci aspettiamo che moti "ragionevoli" passino una frazione di tempo definita in ogni insie- 
me "ragionevole" E dVl; tradotto in termini matematici, ciò vuol dire che, se xe è la funzione caratteristica}^ 
dell'insieme il limite 



(20) lim iV-i V xe{.S^x) = v,{E) 

N^oo ^ — ^ 



L'operatore di traslazione nello spazio delle sequenze "sposta" la sequenza verso sinistra, ovvero ra^ = gi^^i. 

^■^Uno spazio è compatto se e solo se ogni sua collezione di sottoinsiemi chiusi {Ca}a£j con la proprietà di intersezione 
finita è tale che n^gjCa ^ (una collezione di insiemi ha la proprietà di intersezione finita se per ogni sottoinsieme finito J' 
di J si ha che Hcc^j'Ca 0); nel nostro caso la proprietà di intersezione finita è garantita dalla definizione di moto simbolico, 
e l'insieme fi è compatto per definizione. 

^■^La funzione caratteristica xe è definita come 



(19) XeÌx) 



1 se X & E 
altrimenti. 



deve esistere, e con esso la frequenza di visita dell'insieme E da parte moti che hanno origine in x. Se 
E = Pj e un elemento della partizione V = {Pi, P,i} di il il limite (316) equivale a 



(21) 



[Pj) = lim N ^{numero indici h tra e — 1 tali che a^{x) = j}. 



N 



Analogamente, definiamo la frequenza di apparizione nel futuro di un blocco di p simboli ci, cr^ in una 
successione a{x) come 

(22) A^~^{numero di indici h tra e — 1 tali che (Tj.^f^{x) ~ (Tì,ì — 1, ...p}, 

oppure, equivalentemente. 



(23) 



(Ti. ..ap 



N-l 



a{x) = hm iV-i V Xp.i....p {S'^x)- 



h=0 



il numero definito tra parentesi graffe nella ( 22 ) è chiamato numero di stringhe omologhe a 



Ji---]p j appaiono in a{x) tra e — 1. e il valore del limite (23) sarà la frequenza di apparizione 
ai...(7p J 



della porzione di storia ^ ^ J in a{x). 

Riassumiamo queste nozioni nella seguente definizione. 

Definizione 2.2. (Sequenze con frequenze definite) 

Siano d e {0, {ji, jp} € Z e tri, CTp e {0, ...,n} 

ii---jp 



(1) Definiamo il numero di stringhe omologhe al ■^^'"■'p ] ^Jig appaiono in ct tra e iV — 1 come il 
numero di indici h tra e A^ — 1 tali che ^j-^^h = '^ii ■■■T^jp+h — '^P' indichiamo questo numero con 





■Jp 






.ap 





(2) Definiamo la frequenza di apparizione in q_ della stringa omologa a ( ^^'"^p j come il limite, se 

\ (Ti. ..CTp 

esiste, 



(24) 



lim N-^Nì 



N- 



N 



3i--]p 

CTi...CTp 



Jl--Jp 

CTi...CTp 



(3) Una sequenza ct sarà detta a frequenze definite se il limite {24) esiste per ogni ji, ...jp e Z, per ogni 
CTi, ...CTp € {0, n}, e per ogni p — 1,2, .... 

Per concludere la sezione, introduciamo il concetto di ergodicità per sistemi dinamici astratti i cui moti 
possono essere rappresentati simbolicamente. 

Definizione 2.3. (Sequenze ergodiche e sequenze mescolanti) 

Una sequenza ct G {0, ...,n}^ è detta ergodica se ha frequenze definite e se 



lim N-^ V p 

fe=0 



(25) 



V 



]i---Jp 



]i--Jp 
a'i...a' 



il + k...iq + k 



p 



P 



ll...lq 

Gì. ..a 



per ogni p,q = 1,2..., cri,...,crp, cr",...,crq G {0, {ji,...,jp}, C Z. Diremo che a è 

mescolante se 



(26) 



k—yoo \ (T]^...(T 



ri...cr^ 



(Ti ...cr„ 



o; P 



ll...lq 
(Ti ...Gq 



per tutte le possibili scelte di indici. 

Osservazione 2.2. Ovviamente ogni seguenza mescolante è anche ergodica. 

Con la prossima proposizione capiremo l'importanza del concetto di ergodicità; vedremo che l'esistenza 
delle frequenze di visita delle {V, S')-storie dei punti di un sistema dinamico è un fatto molto generale. 

Proposizione 2.2. (Teorema di Birkhoff) 

Sia (Q, S, 11) un sistema dinamico metrico discreto invertibile e sia B la a-algebra .sulla quale fi è definita. 
Se f è una funzione /i-misurabile il limite 



(27) 



N 



N-1 

lim N-^ y f{S^x) = J{x) 

3=0 



esiste fi-quasi ovunque, e se f ii(/i) il limite [21) è raggiunto in ii(/i) 



Corollario 2.1. Sotto le ipotesi della proposizione (2.2), il limite (21) può essere riscritto come 



(28) 



N-l . 

lim V fiS^x) = / fi{dy)fiy) per ogni / e Li (/i) . 

A*— fC30 ^ — ^ J 

3=0 



Dimostrazione. Per una dimostrazione di questo teorema e del suo corollario vedere [9], [8], [10 ]. □ 
Osservazione 2.3. 

SeV = {Pi, ...,P„} è una partizione di $7 in insiemi misurabili il teorema di Birkhoff implica che, scegliendo 
come funzione f la funzione indicatrice Xn^ s-^^p = XnJi jp , 

k = l ^ k -I <T 1 . . . cr ^ 



(29) 



lim iV-i^Xpn..o,(5'=a:)=p ^'-^^ 

N=fOO ^ ' t^ai...ap \ (Tl...CT„ 

fc=0 ^ ^ 



a meno di un insieme di misura nulla; ciò vuol dire che le (V, S)-. storie di fi-quasi tutti i punti di fi 



hanno frequenze definite. Inoltre, grazie al corollario 2.1 
(30) 



lim iV-l Xpn-.. (S'x) = fi (nl^iS-^-P,,) , 

J ' ' CTI ...CT^ 



k=0 



ovvero per fi-quasi tutti i dati iniziali la misura di un insieme è asintoticamente uguale alla frazione di 
tempo che i moti trascorrono visitando l'insieme. 

Questo risultati ci permettono di introdurre due definizioni più naturali di sistema ergodico e sistema 
mescolante. 

Definizione 2.4. Se (il, S, fi) è un sistema dinamico metrico invertibile diremo che è ergodico se per ogni 
E,F G B si ha che 



N~l 

(31) lim iV-i ^ fi{E n S-^F) = ^i{E)^i{F)■ 

invece, diremo che il sistema è mescolante se 



(32) lim ^i{E n S-^F) = fi{E)^i{F). 

j-yoo 



Osservazione 2.4. Si può dimostrare che le condizioni {Sii, (30) sono equivalenti, cf. [8], |10) . i.e. che si 

implicano a vicenda. 

3. Pavimentazioni markoviane 

Come abbiamo visto nella sezione precedente, l'evoluzione di un sistema dinamico in uno spazio fi può 
essere schematizzata attraverso l'azione di una traslazione r su delle sequenze di simboli a_. Chiaramente, 
con ciò non abbiamo reso piìi semplice il problema dello studio dell'evoluzione di un sistema dinamico; 
tutta la complessità è racchiusa ora nel codice che permette di associare ad ogni punto a; e il una 

(7^, 5)-storia in V. Come vedremo, per una classe di sistemi "molto caotici" (che definiremo tra poco) 
risulterà relativamente (e inaspettatamente) facile descrivere le traiettorie mediante dinamica simbolica; in 
particolare, saremo in grado di costruire esplicitamente la misura di probabilità invariante del sistema. 

Prima di capire come sia possibile construire in modo concreto una rappresentazione simbolica, intro- 
duciamo la nozione di compatibilità tra simboli G ir. 

Definizione 3.1. (Matrice di compatibilità) 

Una matrice T (n + 1) x (n + 1) con elementi T^a-' uguali a o a 1 sarà chiamata una matrice di 
compatibilità. Questa matrice sarà detta transitiva se per ogni coppia a, a' esiste un intero a„a' tale che 
T^^!^""' > per ogni cr, a' , mentre sarà detta mescolante se esiste a > tale che Tjj",° > per ogni a, a' ; 
chiameremo a il tempo di mescolamento di T. 

Una sequenza a G {0, ...,n}^ sarà detta T-compatibile o ammissibile se e solo se ogni coppia di simboli 
adiacenti che appare in a è ammissibile, ovvero T^a' = 1 Vi G Z, o equivalentemente se e solo se 



(33) 



n 



Chiameremo {0, 
verificano la (33). 



.,n}^ C {0, ...,n}^ l'insieme (chiuso e invariante per traslazioni) delle sequenze a che 



La nozione di compatibilità tiene conto essenzialmente del fatto che non tutte le traiettorie immaginabili 
sono permesse dall'evoluzione del sistema; l'insieme delle sequenze compatibili equivale all'insieme di tutte 
le traiettorie in iì permesse dalla mappa S. La codifica di una traiettoria in una sequenza di simboli è resa 
possibile da una pavimentazione markoviana. 

Definizione 3.2. (Pavimentazioni markoviane) 

Sia (il, 5) un sistema dinamico topologico invertibile. Dato un pavimento Q = {Qi, Qg} di il, i.e. 
una collezioni di insiemi che ricopre il tale che Qi D Qj = ^ se i j, poniamo T^a' = 1 s4^int((5cr) H 
irLt{S^^Q„>) 7^ e T^^' — altrimenti, con ct, cr' G {1, ...,q},q > 2. Diremo che Q è markoviana se le 
seguenti condizioni sono verificate. 



'L'insieme int(E) C E corrisponde, indicando con dE il bordo di E, a i? \ dE. 



(1) L ^insieme 



+00 

(34) Xfe)= n ^"'Q-^ 

k— — oo 

è non vuoto e consiste di un singolo punto X{a) tale che 

+ CSO 

(35) n ^-=-^+1 = 1- 

i— — 00 

(2) La corrispondenza a — > X{(j} tra {1, ...,g}|, e $7 è Haider continua, ovvero esistono delle costanti 
C, a > toZi che 

(36) rf(^(o:), X(o:')) < Crf(CT; a')", 
dove d è la distanza tra a, a' m {1, q}^ definita come 

(37) d{a, a') = exp {-via, a')) , 

con i^((7, (t') pari al più grande intero j tale che Ui = cr[ per \i\ < j. 

(3) // numero di sequenze compatibili mappate in un certo punto x può essere maggiorato con una 
costante M < 00, ovvero la mappa inversa X^^ verifica la proprietà \X~^{x)\ < M, per ogni 
X € il. Il numero M sarà chiamato molteplicità del codice. 

(4) Ponendo di = dQi, « = 1, q e d = YSl^i di esistono due insiemi chiusi tali che 

(38) 9 = (9+U(9-, Sd-cd-, S-^d+cd+, 

ovvero l'insieme dei bordi degli elementi della partizione può essere diviso in due sottoinsiemi, 
che, rispettivamente, vengono contratti ed espansi dall'azione di S. 

Osservazione 3.1. (l) La continuità di X implica che l'immagine inversa di un insieme boreliano 

tramite X^^ sia ancora un insieme boreliano (perché ogni funzione continua è Borei - misurabile, 

cf. m)- 

(2) X permette di codificare la mappa S in dinamica simbolica, dal momento che X{Ta) — SX(a). 

(3) La definizione di pavimentazione di Markov che abbiamo appena dato in realtà descrive una partizio- 



ne S -separante, cf. definizione 2.1: una partizione non S -separante può essere ottenuta eliminando 



la condizione che l'insieme x(o;), cf. formula {34-), consista di un solo punto 

4l5| 



(4) E chiaro che ogni punto appartenente all'insieme Q. \ U^S'^'J sarà l'immagine di una sequenze^ 
mentre è meno ovvio quale sia il modo di rappresentare i punti che appartengono all'insieme dei 
bordi di Q. Se x £ UiS^^d allora esiste a tale che x G Qa (dal momento che = dQa U 



int((5cr) = int(Qcr)|^^[); poniamo a — gq, e deve esistere a' tale che int(Qo.o) ri mt{S^'^Qa') 7^ e 
Sx £ Qa' ■ Quindi poniamo ai — a' , e per costruzione T^r^ai = 1/ continuando questo ragionamento 
si costruisce la sequenza a G {1, q}^ tale che x S S~^Qa^, per ogni fc e Z. Quindi, x — Xi^u). 
(5) La mappa codice X può essere usata per definire una misura di probabilità nello spazio delle se- 
quenze. Infatti, se m è una misura di probabilità nello spazio delle sequenze compatibili {1, ...,q}^ 
la relazione 

^^Basta compilare una lista con tutti gli elementi del pavimento visitati dal punto x evoluto con S*^, con —00 < k < 00. 
^'^L'insieme E rappresenta la chiusura di E, ovvero E completato con i limiti delle successioni in esso contenute. 



(39) Am{E) = m{X-^E) per ogni E e B{n) 

definisce una misura di probabilità Am su Q. La mappa A trasforma misure invarianti per S 
in misure invarianti per t , misure ergodiche rispetto a S in misure ergodiche rispetto a t ecc. (per 
una dimostrazione rimandiamo a [8].) 

Quest'ultima osservazione è molto importante; come vedremo, la codifica in simboli della misura di 
volume di un sistema "molto caotico" assumerà la forma della misura di probabilità che descrive un modello 
di spin unidimensionale, interagente attraverso potenziali a molti corpi che decadono esponenzialmente. 

3.1. Pavimentazioni markoviane per sistemi di Anosov. Dopo aver introdotto la nozione di 
pavimentazione di Markov e di codice della dinamica simbolica, cerchiamo adesso di capire come costruire 
una pavimentazione di Markov; introduciamo una classe di sistemi dinamici, i sistemi di Anosov o sistemi 
iperbolici lisci, per i quali ciò può essere fatto in modo relativamente facile. 

Un semplice esempio di sistema di Anosov è il gatto di Arnold, cf. TO]; lo spazio nel quale questo sistema 
è definito e il toro T^, ovvero ciò che si ottiene identificando i lati opposti del quadrato [0,2tt] x [0, 27r], e 
l'evoluzione discreta è descritta da 



E facile verificare che la matrice S ha due autovalori distinti A+, A_, — \Z^, tali che 

(41) A+ = 

(42) A_ ^ ^^; 
gli autovettori v_ corrispondenti a A+, A_ sono dati da 



(43) v_, = 

(44) ^ ( 1 ) conA = A_. 

Dal momento che i;_ hanno pendenza che non è multiplo razionale di tt, le due rette W"{ip), W^{ip) 
riempiono densamente il toro T^, cf. |10| . e inoltre W"{(p), W{(p) sono (banalmente) covarianti, ovvero 
SW°'{ip) — W'^{S(p) {a = u,s); poiché A+ > 1 e A_ < 1, due punti sulla stessa varietà stabile W, o sulla 
stessa varietà instabile W", si avvicineranno , o allontaneranno, esponenzialmente. Tradotto in termini 
matematici, per ogni rt > si ha che 



(45) d(S'"?A,5V) < A"d(^,V;') per ogni iP, ^^/ e W (ip) 

(46) d{S-"ìp,S-"-i/) < X^d{±,i/) per ogni e W^"(^), 

se d è la distanza misurata sulla varietà stabili e instabili. Quindi, poiché le varietà stabili e instabili 
sono dense in T^, ogni punto del toro sotto l'azione di 5 è fortemente instabile rispetto a qualunque altro; 
questa situazione giustifica l'aggettivo iperbolico. 

Ispirati da questo esempio, possiamo dare una definizione generale di sistema di Anosov (o sistema 
iperbolico liscio). 



Definizione 3.3. (Sistema iperbolico liscio o sistema di Anosov) 

Sia {fi, S) un sistema dinamico definito su una varietà riemannian^^il compatta e connessa di classe 
C°° con S un dijfeomorfismo di classe C°° . Supponiamo che il sistema sia topologicamente transitivo, 
ovvero che esista un'orbita densa. Supponiamo inoltre che esista una metrica riemanniana regolare d (la 
quale può essere diversa da quella definita su fi) tale che misurando le lunghezze con d le seguenti proprietà 
sono verificate. 

(1) Esistono due superfici W'^{x) e ^{x), che supporremo di classe con k > 2, e con piano 
tangente in x dipendente da x in modo Haider continuo. Inoltre, W^{x) e W''{x) sono trasvers^^ 
in X e hanno dimensioni positive complementar^^ 

(2) Chiamando la sfera di raggio 7 centrata in x e ponendo 

W^{x) — {parte connessa di W'^{x) D S^(a;) contenente x} 
W^{x) = {parte connessa di ^{x) D S^(a;) contenente x}, 
esiste 7 > tale che 

(47) sw;{x) c w;{Sx), s-^w:;{x) c w:;ìx) c w:;{s-^x). 

(3) Esiste A < 1 tale che, per ogni n> 0, 

(48) d{S"y,S"z) < X"d{y,z) per ogni y , z £ (x) 

(49) d(S'-"y, S'-"z) < A"d(y,z) per ogm y, z e Wl^ix), 

(4) Esiste £ > 0, e < 7 tale che, se x,y £ Q e d{x,y) < e, l'insieme Wl^{x) H W^{y) consiste in un 
singolo punto [x,y\, che dipende in modo continuo da x e y. 

Sotto queste ipotesi diremo che (fi, S) è un sistema iperbolico liscio o sistema di Anosov e che la metrica 
riemanniana d è adattata alla mappa S . 

Osservazione 3.2. (1) Grazie alla proprietà Q) la varietà stabile e la varietà instabile defi- 

niscono un sistema di coordinate locale. 
(2) La definizione di sistema di Anosov che abbiamo appena dato vale per i sistemi dinamici discreti; 
ad ogni modo, è possibile formulare una definizione analoga per i flussi, e la differenza principale 
è che oltre alla direzione stabile e alla direzione instabile esiste una direzione marginale, lungo la 
quale due punti distinti non si avvicinano né si allontanano esponenzialmente. Questa direzione 
coincide con la traiettoria del punto nello spazio delle fasi. 

Prima di passare ad una dimostrazione costruttiva dell'esistenza di una pavimentazione markoviana per 
i sistemi di Anosov, introduciano il concetto di S -rettangolo. 

Definizione 3.4. (S'-rettangoli) 

Un insieme R (Z fi è un S-rettangolo per il sistema iperbolico {fi, S) se 

(1) R = int(iì), 

(2) il diametro di R è abbastanza piccolo da permettere che se x,y € R esiste 7 tale che l'intersezione 
di W!^(a;) e W^{y) definisce un unico punto [x,y] G R. 



^''Una varietà riemanniana è una varietà differenziabile dove in ogni spazio tangente è definito un prodotto interno, in un 
modo che varia analiticamente da punto a punto; ciò ci permette di definire gli angoli, la lunghezza delle curve ecc. In altre 
parole, una varietà riemanniana è una varietà differenziabile in cui lo spazio tangente in ogni punto e uno spazio di Hilbert. 

^^Ciò vuol dire che si intersecano con un angolo di incidenza non nullo. 

l^Owero dim(n) = dim(VK") + dim(VK=), dim(VK"'^) > 0. 



Osservazione 3.3. Ad esempio, nel caso di un sistema di Anosov bidimensionale (¥^,5*) possiamo "prolun- 
gare" in modo opportuno due "lacci" diretti lungo varietà stabile e instabile in modo da ottenere una griglia 
dove ogni maglia è un S-rettangolo. 

Siamo pronti per dimostrare il seguente fondamentale risultato. 

Proposizione 3.1. (Esistenza delle pavimentazioni markoviane) 

Se (fi, S) è un sistema di Anosov e S > 0, allora esiste una pavimentazione markoviana Q — (Qi, Qq) 
di n dove Qi, i — 1, ■■■,q, sono S-rettangoli con diametro minore di S. 

Osservazione 3.4. Per comodità, dimostreremo questa proposizione nel caso d = 2. L'argomento che 
presenteremo non è il più facile per questo caso, ma ha il pregio di essere facilmente estendibile a qualsiasi 
dimensione. 



Dimostrazione. La dimostrazione della proposizione (3.1) risulta essere molto più semplice di quanto 
non possa sembrare a prima vista se si rappresenta graficamente ciò che si sta descrivendo a parole; per 
questo motivo, cercheremo sempre di illustrare i punti chiave con dei disegni. 

3.1.1. Descrizione di un rettangolo generico. Come abbiamo già accennato in precedenza, varietà stabile 
e varietà instabile definiscono un sistema di coordinate locale. Quindi, se A è un ^-rettangolo e x € int(A), 
ponendo C = W"{x) H A, D — W^{x) n A abbiamo che (vedere figura [l| 

(50) A=[C,D]=Uyec,zeD[y,z]; 

chiameremo x il centro di A rispetto alla coppia di assi C, -C0 e diremo che C, D sono rispettivamente 
l'asse stabile e l'asse instabile di A. Se C, D e C , D' sono due coppie di assi per lo stesso rettangolo diremo 
che sono paralleli, dal momento che o C = C oppure C n C" = 0. Identificheremo con d^A, 82 A i bordi del 
rettangolo giacenti sulle varietà instabili, e con dlA, quelli appartenenti alle varietà stabili. 

Vogliamo dimostrare l'esistenza di un pavimento Q di fi formato da ^'-rettangoli con diametro < 5 
arbitrario, tale che per ogni Q € Q si abbia 

(51) Sd^Q c d'p.Q', S-^d}Q c a^„Q", 
con 13, 13', 13" = 1,2. 



Notare che, in realtà, la proprietà che definisce un pavimento markoviano è (vedere definizione 3.2 1 
(52) S{iJQ^Qd'Q)ciiJQeQd'Q', (Ugee) C UgesS^Q", 



ma è facile rendersi conto che la (51 1, (52) sono equivalenti. Infatti, è chiaro che ponendo d^Q 



d{Q U e 9"Q = d^Q U 92 Q la (51) implica la (52); l'implicazione inversa può essere vista con un 
semplice ragionamento per assurdo. 



Applichiamo la mappa S ad un elemento Q € Q; per la proprietà (52), l'insieme dei bordi stabili 
sotto l'azione di 5* diventa strettamente contenuto in sé stesso, e la stessa cosa accade per 9" sotto l'azione di 
S~^. Assumiamo che la (51) non sia vera; dunque in generale avremo che Sd^Q C Ug/gc^^'Q' dove C C Q è 



una collezioni di insiemi contenuta in Q con \C\ > 1. Senza perdita di generalità, fissato d^Q consideriamo il 
caso in cui C sia composto da due elementi Q' , Q" con un bordo instabile in comune, che chiameremo d^Q' , 
e pensiamo che Sd^Q C d^Q' U d^Q" . Ciò implica che è compreso tra SdfQ e 5*91(5 generando un 

assurdo, perché S^^ manda punti relativamente internj^a Sd'^Q in punti relativamente interni a d^Q e per 
ipotesi 5-19" C 9". 



■^"Con questa definizione ogni punto interni di a; è un centro di A rispetto ad opportuni assi. 

■^^Per punti relativamente interni a 9"Q, o a d'Q, intendiamo i punti compresi tra 9"Q e 92*3' ° dfQ e 9|Q; il punto 
X è compreso in d^Q = d^QUd2Q se è possibile raggiungere 9"Q percorrendo la varietà stabile che passa per x, ed è compreso 
in cJ^Q = 9|Q U 9|Q se è possibile raggiungere d^Q percorrendo la varietà instabile passante per x. 



Figura 1. Un 5'-rettangolo con assi C, D e centro x. 



3.1.2. Un problema equivalente. In realtà, dimostreremo la proprietà (51) per la mappa g — S*™ con m 



"grande" (in un senso che specificheremo in seguito), ovvero troveremo un pavimento B àìVL tale che 

(53) g Ub6B d'B C \JBeBd'B, g-^ UseB C UbebO^'B; 

ma questo problema è equivalente a quello di partenza, dal momento che da B si ottiene un pavimento 
markoviano per S semplicemente ponendcj^ 

(55) Q = r^T=oS''B,,, 

con Bi-^, Bii_ e B, per tutte le possibili combinazioni di indici. 

3.1.3. Costruzione di un ricoprimento di rettangoli. Introduciamo la distanza di Lebesgue di un ricopri- 
mento A'^ come la massima distanza a tale che la sfera di raggio a disegnata prendendo come centro un 
qualsiasi punto di 17 è strettamente inclusa in un elemento di A'^. Nel nostro caso modificheremo questa 



■^■^Infatti, con questa scelta di Q esiste B' tale che 



(54) = UB^^eB-UB^^^_^eBd'n'^-l S''B,^=UQeQd'Q; 

un ragionamento dello stesso tipo può essere fatto per la seconda proprietà nella \52h. 



definizione, sostituendo al posto della sfera le porzioni locali (c/. definizione 3.3) delle varietà stabili e 
instabili W^(x), W^{x) passanti per x; infatti, richiederemo che il nostro ricoprimento A'^ ci permetta di 
trovare per ogni a; e fi un ^^(^) G -4° tale che 



(56) diy,W;'{y)nd'AF(^^)) > a e W^{x) n A%^^^ 

(57) d(z,VF^^(z)n9Mj.(,)) > a e W^fix) n A%^^^, 

dove le distanze d sono misurate rispettivamente lungo W^{y) e W^{z). In pratica, stiamo escludendo 
il caso in cui le varietà locali W^{x), W^{x) si trovino "troppo" vicine ai bordi di tutti gli elementi di A°. 




Figura 2. Per ogni punto x E il esiste un S'- rettangolo Ap(^^-j tale che la varietà stabile e 
la varietà instabile passanti per x sono "abbastanza lontane" dai bordi di Ap(^^y 

Qui e nel seguito assumeremo sempre che il diametro dei rettangoli sia abbastanza piccolo da poter 
considerare W^{x), W^{x) come dei segmenti dritti di lunghezza 2a centrati in x; quindi, assumendo che il 
diametro massimo dei rettangoli di A*^ sia a < 6 /2 (5 è il diametro massimo dei rettangoli del pavimento di 
Markov che vogliamo costruire; la ragione della scelta di a sarà chiara nel prossimo passo della dimostrazione) 
fisseremo < a < a/2. 

3.1.4. Raffinamento del ricoprimento A'^ . Adesso, vogliamo dimostrare che partendo dal ricoprimento 
A'^ appena introdotto è possibile costruire un nuovo ricoprimento A che verifica le seguenti proprietà: 

(1) per ogni l3 = 1,2 e Ae A esistono /3', /3" = 1, 2 e A', A" e A tah che 

(58) gdf;Acidf,,A', g-'d]^Acd}l„A"; 

(2) per ogni z £ d'^Ai, y e d^'Ai si ha che 



(59) W^a/2(5~^y) C per j(z), opportuni. 

Ciò potrà essere fatto nel seguente modo. Per ogni A" G consideriamo il segmento che collega 
due punti ^ € 9f ^' e SfA" giacenti sulla stessa porzione di varietà instabile (la scelta della coppia 

di punti è indifferente); questi segmenti formano un insieme = {C^Y^j=\i dove r è il numero di elementi 
di . Poiché m è "molto grande" (tra poco spiegheremo in che senso), i "nuovi" segmenti .gC°, j = 1, ...,7' 
saranno dei "fili molto lunghi", giacenti sulla varietà instabile. Scegliamo un ricoprimento di 

gC'^ composto da elementi di tali che gCj* li attraversi ad una distanza > a dai bordi paralleli, con a pari 
alla distanza di Lebesgue definita in precedenza. In generale, gC^ non attraverserà completamente alcuni 
rettangoli che lo ricoprono; ciò accadrà per i rettangoli , che contengono le estremità del "lungo filo" 
(7C*, c/. figura [sj In questo caso, "allunghiamo" il filo in modo da far coincidere le sue estremità con un 
punto appartenente S^^ìq e uno in d^iAi^, , ovvero introduciamo un nuovo segmento Cj coma^'^ 




Figura 3. Il segmento Cj diventa un "filo molto lungo" sotto l'azione di g — S™'; nella 
figura, il rettangolo non è completamente attraversato da gCj 

cj = U^'^j^g"^ ^{continuazione di gCj che attraversa tutti 
(60) i rettangoh n <) - u'^^.g-' ([C^ ,gC° n <]) , 

■^■^L'insieme [Cf ,gC^ n A? ] corrisponde a U f ^0 ^o \ l^ih'^ì' dove [Ci,,x] è il segmento passante per x parallelo a 

V J Ih/ 

Cif^ , dunque [C?^ , gC^ fi A?^] è il prolungamento di gCj fi A^^ che interseca sia d"A che 9|A. 



e sostituiamo il vecchio con — [Cj,D^] {Cj è leggermente "più lungo" di Cj, dunque definiamo 
un nuovo A^ che lo contiene). 



Adesso, applichiamo g al segmento Cj definito nella ( 60 1 , e ripetiamo il discorso per questo "nuovo filo" ; 
iterando n volte si ha che 

(61) = ut,g~' m-\gC--' n Al;']) , 

e resta da discutere se questa procedura ha senso nel limite n — oo, i.e. se il filo converge ad un "filo 
asintotico" . 

Ciò può essere fatto con un ragionamento per induzione. Al termine del primo passo, quello che porta 
da C° a Cj, il "nuovo filo" Cj potrà essere espresso come Cj = gi U C° U dove qi, sono i due 
"segmentini" vengono aggiunti alle estremità di Cj per permettere a gCj di attraversare completamente 
tutti i rettangoh del ricoprimento {^^^j, corr^(7i,gj < aA™; definendo AJ = [C",D°], e assumendo che al 
passo n si abbia |g„| < (A™)" a {idem per 5„)7al passo n + 1 il "filo" ffC™, che attraversava completamente 
i rettangoli A^~', dovrà essere "prolungato" ad ogni estremità con dei segmentini q„, q„, ovvero C"^^ = 

g-'qn U U g-'q„ = qn+i U U con |g„+i|, |g„+i| < (A™)"+' a. 

Grazie a questo argomento è chiaro che il limite Cj — U^qC*" avrà lunghezza finita, dal momento che 
< \C°\ + T,n=i 2aA™" = |C°| + < C50, e scegliamo m in modo tale che 

2aA'" 

(62) < a/2 < a. 

Adesso, costruiamo i segmenti {Dj}^^-^ partendo da {D^}j^^ con un procedimento del tutto simile a 
quello usato per {Cj}^^i, sostituendo g con g^'; il risultato è che possiamo definire un nuovo ricoprimento 
A con Aj = [Cj, Dj], diam(74j) < 2a ^ 6, tale che 

(63) g[C,,z]D[C,-,gz], g-'[z, D,] D [g-' z, D,.], 
e inoltre 



(64) A,, D Wiiy) yye[C,,,gz] 

(65) Aj„ D Wli{z) \fz e[g-'z,Dj„]. 



Queste proprietà coincidono con le (58), (59), dunque abbiamo costruito il ricoprimento con le carat- 
teristiche richieste; i rettangoli che formano questo ricoprimento, però, in generale avranno punti interni 
in comunque. Nella prossimo punto della dimostrazione faremo vedere come da questo ricoprimento sia 
possibile estrarre una pavimentazione markoviana. 

3.1.5. Costruzione della pavimentazione di Markov. Per concludere la dimostrazione, facciamo vedere 
come da A sia possibile ottenere un pavimento di Markov B, formato da S'-rettangoli Bi con diam(_Bi) < 
2a = S. 

Partendo da A, possiamo ottenere un raffinamento di A intersecando elementi di A con i complementari 
di altri. Più precisamente, per ogni r C {l,...,r}, possiamo definire Pr = (Higr^i) H (Uj^r^\^jQ per 
costruzione si ha che (come per A, i bordi di Pr appartengono alle varietà stabili e instabili) 



^^Non è detto che il C'j che compare il questa definizione di Cj sia lo stesso di partenza; in generale sarà un segmento 
parallelo al "vecchio" C?, con la stessa lunghezza. 

■^^a è il massimo diametro dei rettangoli A'^ £ e il fattore A™ è dovuto all'applicazione di g~^. 
^^Quindi, se r ^ r' si ha che Pr n P^' = 0- 



(66) 




Figura 4. Non tutti gli insiemi generati dalle intersezioni di elementi di A con i 
complementari di altri sono iS-rettangoli. 



Questo pavimento, però, non è formato da S'-rettangoli, perché non è vero che se x,y Cz P allora [x, y] G P 
(vedere figura [4|. Per ottenere un pavimento formato da quadrati, basta "continuare un po'" i lati degli 
elementi dì Ai A fino a quando il prolungamento del lato non incontra un lato del "tipo opposto" (ovvero 
un lato stabile per un prolungamento instabile, e viceversa) (vedere figura [s]); chiaramente, dal momento 
che il diametro dei rettangoli contenuti in è < 2a, la lunghezza del prolungamento sarà < 2a. Più 
precisamente, ciò vuol dire sostituire d%A, da A con da A, d^A definiti come 



(67) 



ovvero continuiamo d'^A, d^A ad ogni estremità aggiungengo un pezzo della varietà dello stesso tipo 
di lunghezza 7. Poiché 7 > 2a, gli insiemi definiti nella (67) andranno oltre il punto dove intersecano per 
la prima volta un lato di tipo "opposto", dunque dobbiamo considerare solo la parte di d^A, d^A che non 
supera il punto prefissato; con questi "nuovi lati", gli insiemi P^, che chiameremo B € B, saranno degli 
^'-rettangoli con diametro < 2a ~ S. L'ultima cosa che rimane da fare è verificare che gli insiemi B E B 
rispettano le proprietà del pavimento di Markov ( 53 ) ; per fare ciò, notiamo che per ogni punto z appartenente 
al bordo di B, ad esempio al bordo stabile d'^B, esiste un ce G 9*^g(z)j con G(z) opportuno, tale che 



(68) 



Figura 5. Prolungando i lati degli 5-rettangoli Ai ^ A in modo opportuno si ottiene un 
nuovo pavimento B formato da ^-rettangoli che non si intersecano tra di loro. 



e lo stesso vale se z ^ d'^B, sostituendo WIqI^) con W2a{y) per un certo y e d^A^iz)- Grazie alle 



proprietà (59) abbiamo che, per j{x) opportuno, 



(69) ga;G9Mj-(^) e W^/^igx) C A^i^^y, 

poiché il diametro di gW2a{x) è minore di < a/2 (per la scelta di m, vedere eq. (|62])), abbiamo 

che gW2aix) C W^^2Ì9^)- Dunque, grazie alla (691 e sapendo che UAeAd^'A C UBeBd^B (lo stesso vale 
sostituendo s con u), 



(70) gdIsB C UB'eBd'B'; 

allo stesso modo è possibile dimostrare che 



(71) g-'d}^B c UB'eBd^'B'. 

3.1.6. Haider continuità del codice. E facile vedere che la corrispondenza che associa ad ogni sequenza 



(7 un punto X è una funzione Hòlder continua, nel senso della definizione (35|. Infatti, assumiamo che due 



punti xi ~ ^(fl), ^2 = ^(^l') abbiano la stessa storia CTi, ...,(T„ fino al tempo n > 0, e che al tempo n + 1 
Cji+i 7^ cr^i+i! ciò vuol dire che chiamando W]^2(^i) il vettore che congiunge i due punti 



(72) 



Possiamo scrivere il vettore iiy^ come 



(73) ^^12(2^1) = CuiXi)WuiXi) + Cs{xi)WsiXi), 

dove Wu{xi), Ws{xi) sono i due vettori normalizzati tangenti alla varietà instabile e alla varietà stabile 
nel punto xi; grazie alla (73l la (72| diventa, assumendo per semplicità che S agisca linearmente sui vettori 
echfF]c„^0, 



(74) 



|c„(xi)S'"w„(a;i) + Cs{xi)S"-Ws{xi)\ < 6. 



(75) 



Quindi, per la definizione di sistema di Anosov |3.3| la (74) implica che 

< Sd{xi,X2) 



\cu{xi)\X-"dixi,X2) ( i-A^""^^^""!^ 



ovvero, infittendo opportunamente 
C > tale che 



28 



Cuixi) 

la partizione in modo da avere A^" 7^ 1, esiste una costante 



(76) 



d{xi,X2) < A"C; 



infine, ponendo c — — log A, possiamo riscrivere la ( 76 ) come 



(77) d{xi,X2) = d{X{a,X{a'))) < Ce"^" = Cd{a,a'). 

Inoltre, scegliendo una partizione abbastanza fitta l'intersezione tra un elemento di P e l'immagine sotto 
S di qualunque altro sarà un insieme connesscj^ dunque in questa situazione non può avere molteplicità 
infinita; ciò conclude la dimostrazione. 

□ 

4. Codifica della misura di volume 

In questa sezione mostreremo che, per un sistema di Anosov, la codifica in simboli dell'evoluzione di 
un dato iniziale permetterà di rappresentare la misura di volume su il come la misura di probabilità di 
un modello di spin unidimensionale con potenziali a decadimento esponenziale. A tal fine, introduciamo la 
nozione di probabilità condizionata. 

Definizione 4.1. (Probabilità condizionata) 

Sia m una misura di probabilità definita sui sottoinsiemi boreliani dello spazio delle sequenze con q 
simboli {1, q}^, e siano J — {ji, js} C Z, CTj = (cj^ ....(Tj^ ) e {1, q}"^ . 

Consideriamo le a-algebre B{J'^) generate dai cilindri C^, , con base J' C J*^ = Z \ J ; dato q_j possiamo 
definire su B{J'^) le misure 



(78) E -^m'{E) = m{EnCl^), E fn{E) = m{E) , 

per E e B{J'^). Dunque la misura m' è assolutamente continugp^ rispetto a ni, ed è possibile definire la 
derivata di Radon-Nikodym di ni' rispetto a ni 

^^È possibile infittire abbastanza la partizione in modo da escludere l'eventualità che sotto l'azione di una iterata S due 
punti giacenti sulla stessa varietà instabile e nello stesso 5-rettangolo si spostino in un S-rettangolo comune. 

■^^In particolare, dobbiamo infittire la partizione evitando che si formino dei rettangoli troppo "lunghi e stretti". 

■^'^Stiamo escludendo l'ipotesi che sotto l'azione di S un rettangolo si allunghi abbastanza da poter intersecare due volte lo 
stesso rettangolo. 

'^'^Due misure v, /i definite in Q sono assolutamente continue se e solo se i'(A) = -^=> /^{A) = 0; il teorema di Radon 
Nikodym afferma che per due misure assolutamente continue è possibile trovare una funzione misurabile / tale che 



(80) 



dm' 
dm 



(ct') = m{g_j\q^j,)] 



l'espressione è chiamata la probabilità condizionata di g_j rispetto a (ij. 
Osservazione 4.1. La probabilità condizionata m{a_j\a^ja) deve verificare l'ovvia condizione di compatibilità 



(81) 



^m{aj\q^j.) = 1; 



quindi, grazie alla, {81), per sapere chi è 'rn{q_j\q^ja) è sufficiente conoscere i rapporti 



(82) 



mia j\(T 



Come abbiamo già detto introducendo i sistemi di Anosov, ha senso considerare la mappa S come 
una mappa S : W"{x) — )• W^{Sx), oppure come S : W^{x) — >■ 1^^(5*2;); chiamiamo coefficiente locale 
di espansione e coefficiente locale di contrazione il modulo del determinante della matrice jacobiana di S 
rispettivamente lungo la varietà instabile e stabile. Nel seguito considereremo per semplicità solo il caso 
bidimensionale . 

Definizione 4.2. (Coefficiente di espansione locale e coefficiente di contrazione locale) Sia {il, S) un sistema 
di Anosov bidimensionale, e sia Q = {Qi, ...,Qq} un pavimento markoviano. Sia T la sua matrice di 
compatibilità, e X il codice della dinamica X : {1, q}^ — > fi. Definiamo il coefficiente di espansione locale 
A„(a;) e il coefficiente di contrazione locale Xs{x) come 



(83) 



A-'(x) 



dS' 



d^ 



-(x) 



Xs{x) 



dS 
di 



(x) 



dove con i pedici u, s indichiamo che la derivata di S è calcolata sulla varietà instabile o sulla varietà 
stabile, e ^ rappresenta la parametrizzazione della varietà; inoltre, definiamo il tasso di contrazione locale e 
il tasso di espansione locale come 



(84) Au{x) = logA„(a;), As{a) ^ -\og\s{x). 

Osservazione 4.2. (1) La definizione di sistema iperbolico liscio implica che Xu{x), \s{x), e dunque 



anche i logaritmi (84-), siano funzioni Haider continue di x. 



(2) Nel caso del gatto di Arnold i coefficienti di espansione e di contrazione coincidono con V autovalore 
instabile A+ e con V autovalore stabile A_. 

(3) Come vedremo tra poco, partendo da Au{X{a)), As{X{a)) potremo definire i potenziali a molti 
corpi del modello di Ising unidimensionale associato alla misura di volume. Inoltre, grazie alla 
Haider continuità di \u{x), Xs{x) e del codice X, A^ e As sono funzioni Haider continue nella 
spazio delle sequenze; questa fatta implicherà il decadimento esponenziale dei potenziali associati 
alla misura di volume. 

(4) Nel seguito, quando scriveremo Au{(j), As{g() ci riferiremo in realtà, con un leggera abuso di 
notazione, a yl„(X((T)), As{X{<j)). 



(79) 



u{A) = [ fd,,. 

JA 



Con l'osservazione [3] abbiamo anticipato il seguente risultato. 

Proposizione 4.1. (Sviluppo in potenziali di funzioni Hòlder continue) 
Le funzioni Au{<7), As{a) introdotte nella definizione 



(85) Aa^ C0St+^^2n+li<^-n,--,<^n), a^U,S, 

n=0 

dove {3*2n+i} decadono esponenzialmente, ovvero esistono C,k> tali che 

(86) |$^„+i(a_„,...,a„)|<Ce-^-". 

Osservazione 4.3. Le funzioni ^2n+i sono dette cilindriche, dal momento che dipendono da un numero 
finito di simboli (T_„, (T„. 

Dimostrazione. La dimostrazione di questa proposizione è molto semplice, e si basa su un argomento 
piuttosto ricorrente nella teoria dei sistemi di Anosov. Consideriamo le sequenze a; e {!,..., q}^, e sia z il 
tempo di mescolamento della matrice di compatibilità T, ovvero tale che T^^, > per ogni a, a' € {1, g}; 
quindi, per ogni cr, a' esiste almeno una stringa z9(cr, a') formata da z simboli compatibili tale che la successione 
a-d{a,a')a' è compatibile. 

Costruiamo le sequenze o;" come, indicando con o_^^ ^ la frazione della stringa a_ compresa tra gli indici 
a, h (inclusi), 

(87) O;" = CT[_oc^_„__j„i]?9(CT_„_^_i,Cr°,Jo:"_„„]l?(CT°,CT„+2 + i)CT[„+^+l_+^^ 

ponendo Aa{g^) — cast possiamo scrivere la funzione Aa{o) come 

oo oo 

(88) AJa) = cost + ^ (Aa(a"+i)-Aa(a")) =cost + ^$^„+i(a_„,...,a„). 

ri=0 n=0 

Infine, grazie alla Hòlder continuità di Aa{a) esistono due costanti C, fc > tali che 



4-2, possono essere scritte com 



(89) |<l>L+i| = \Aa{ct+') - A,(a")| < Ce-'=", 

e ciò conclude la dimostrazione. 



□ 



Abbiamo ora gli strumenti per capire cosa diventa la misura di volume di un sistema di Anosov se 
"codificata" in dinamica simbolica. 

Proposizione 4.2. (Codifica della misura di volume) 

Sia (fi, S*) un sistema di Anosov bidimensionale, e Q — {Qi, ■■■,Qq} un pavimento markoviano con 
matrice di compatibilità T e codice X. Sia /i" la misura di volume in ft associata alla metrica definita in Q, 
e TTiQ la sua misura "immagine" nello spazio delle sequenze {l,...,q}^ tramite X (ci. osservazione^ nella 



nizione 3.2). I rapporti di probabilità condizionate calcolate con la misura TOq verificano 



In generale ciò sarà vero per qualunque funzione Hòlder continua. 



(90) 



exp 



oo 

E 

k = l 



mo(crL„,...,g^|gj, \j\ > n) ^ sin y;(X(g')) 
mo(fT^„...,<|a„|j| >n)- sin^(X(a")) 

^„(rV) - ^„(rV") + ^.(r- V) " ^.(r" V") ) , 



delie </j(a;) è l'angolo tra la varietà stabile e la varietà instabile nel punto x, e a' , ct" sono sequenze in 
{1, ...,q}^ i cui valori con indice j coincidono per \ j\ > n. 



Osservazione 4.4. Grazie alla proposizione 4-1 possiamo riscrivere la (TPOl) in modo che assuma una forma 



reminiscente della teoria dei modelli di spin. Definiamo le funzioni ('&x)-f cz nel seguente modo: 



(91) '^xi^x) = a meno che X ^ {lj + l,...,l + 2n) 

(92) <ì>'k{ax) - nn+iilx) seX = {l,l + l,...,l + 2n), 
per a_x G {l,...,q\'^ e qualche l,n (£ Z; con questa definizione si ha che 



(93) Aa{a)= J2 '^xilx). « = ",5. 

X centrato Ì7i 

Quindi la può essere riscritta come, se a{X) = u per centro(X) > e a{X) = 5 per centro(X) < 0^ 



(94) 



^o{'y'{_„^„]\<Jj, \j\ > n) 



cxp 



E 

Xn[-n,n]5^0 



aiX) 
X 



i<lx)-'^>''^''\^x) 



B{a) 



dove il termine B è dovuto al rapporto dei seni e al fatto che nella (90) la somma su k parte da 1; quindi, 



poiché B{a) è Haider continuo può essere sviluppato a sua volta in potenziali a decadimento esponenziale 
centrati in 0. La {94) è a tutti gli effetti il rapporto delle probabilità di due configurazioni di spin a', a" 
di un modello di Ising unidimensionale con condizioni al bordo fissate al di fuori di un volume [—n,n\, e 
interagente attraverso potenziali a molti corpi diversi se il centro dell'intervallo reticolare X si trova a destra 
a sinistra dell'origine 0. 



Notare che nella (94), come nella (90), le somme all'esponente sono assolutamente convergenti perché 
le quantità sommate sono calcolate in configurazioni a, gf diverse solo in un intervallo finito [~n,n\, e le 
funzioni sono Haider continue. 

Dimostrazione. Come per l'esistenza delle partizioni di Markov, anche in questo caso la dimostrazione 
della proposizione diventa notevolmente piii semplice se visualizzata graficamente; per questo motivo, cer- 
cheremo sempre di illustrare i punti chiave con delle figure. Otterremo la (90) nel caso n = 0; il caso n > 



si studia con un argomento del tutto analogo a quello che stiamo per presentare. 

Consideriamo due sequenze qf,qf' S {!,..., g}f. con a'^ = a'^ = aj per |j| > 0, e poniamo X{a) 
X{qf) — y; assumiamo per semplicità che x,y ^ ^i=idQ', è chiaro che i due punti x,y appartengono alla 
stesse varietà stabili e instabili, dal momento che hanno storie uguali nel passato e nel futuro. Per il teorema 
di Doob, cf. [8 , si ha che, /i-quasi ovunque, 



(95) 



(r~i-N...-l 1...n'\ 



— lim 

WO {(t'Ì I CTj , I j I > 0) N^oo ^(j^ 



N...-1 I...N 
(T-N...{T-i{TL'ai...aN 



quindi, per dimostrare la (90) dobbiamo semplicemente calcolare il limite nella parte destra della (95). 



L'area dell'S-rettangolo ^" ^, g dunque la sua misura di Lebesgue mo, può essere approssimata 

con 



(96) 



A 



N...N 

-JV---0"'-, 



\S:\-\S^\sm^{x) + oiX 



2N\ 



dove (5^, (5g sono vettori tangenti alla varietà instabile e alla varietà stabile applicati in un vertice 
x' (non importa quale) dell'S'-rettangolo che congiungono x' ai punti x" , x'" (vedere figure |6j [?]), e ipjx) 
è il seno dell'angolo compreso tra S^, (5f; nella formula (96) abbiamo sostituito ip{x') con (p{x) perche^ 
|sin(/j(a;) -sin(^(x')| < CA^^, con c, C > 0. 




Figura 6. Rappresentazione schematica del cilindro C^j'^;;^o...ctjv ■ In questo disegno ab- 
biamo approssimato il cilindro con un rettangolo, e questo è possibile a meno di o(A^); 
vedere figura [7] per una rappresentazione piìi realistica della situazione. 



Lo stesso argomento può essere ripetuto per il cilindro C^^^"^^,, e quindi il limite nella parte destra 



della ( 95 ) coincide con 



(97) 



smip{y) N^oc \Sl \ • I ■ 



possiamo riscrivere i rapporti nella (97) come, poiché S o S — 1 



Ciò è dovuto al fatto che W'^{x), W^{x) dipendono in modo Hòlder continuo da x. 



Su 



i-N. ..N 



Figura 7. 



(99) 



m ^ Y\U^s{S'^x')\s-H^^ 



e poiché k è arbitrario poniamo k = ujN con uj G (0,1). Per — > oo i punti x', y' tendono a a;, 
y, i quali si trovano sulla stessa varietà stabile e instabile dal momento che x = X{q^), y = X(a") con 
g,- = cr" se |i| > 0. Quindi, grazie alla Hòlder continuità di A„(x), As(x) i rapporti di prodotti nelle (98 1, 



(|99|) tenderanno ad un limite finito; l'ultima cosa che rimane da discutere è il comportamento di Ig^^^yj . 



I S u 

|g-^iv^|[ per iV ^ co. I due vettori S'^'^ S^, S'^^ 61^ diventano "molto vicini" per N grandi, nel senso che la 

distanza tra x e y è 0(A"^) (inizialmente è 0(1)), e inoltre — 0(A^^~'^''^); quindi, poiché l'angolo 

ip{x) Hòlder continuo in x, 

Un argomento dello stesso tipo può essere ripetuto per 6^, S^, quindi il risultato finale è che 
e la (90) è dimostrata ponendo A„(ct) = log A„(X(ct)), As{(j) = —\ogXs{X{a)). 



Figura 8. Avvicinamento dei segmenti S"^^^^, S'"^5^. 

□ 

Quindi, come abbiamo già notato con l'osservazione [4^ possiamo immaginare la misura di volume codi- 
ficata in dinamica simbolica come la misura di probabilità di un modello di Ising unidimensionale interagente 
attraverso potenziali a molti corpi diversi a seconda che ci si trovi a destra o a sinistra dal centro del volume 
nel quale è confinato il modello. E abbastanza naturale chiedersi se nel limite in cui ci si sposta infinitamente 
lontano verso destra o verso sinistra, i. e. nel passato o nel futuro del sistema di Anosov, la misura di proba- 
bilità del corrispondente modello di spin sia determinata da potenziali d'interazione di un'unico tipo; nel caso 
ciò fosse vero, la misura non potrebbe che coincidere con la codifica in simboli della misura di probabilità 
invariante del sistema dinamico. 

Rispondere a questa domanda equivale a risolvere il problema dell'esistenza del limite termodinamico per 
un modello di Ising unidimensionale con potenziali a decadimento esponenziale; ma dalla teoria dei sistemi 
di spin sappiamo che questo limite esiste e che non dipende dalle condizioni al bordo, dunque nel nostro caso 
dai potenziali localizzati nel passato e nell'origine. Riprenderemo il problema in termini piìi precisi nella 
prossima sezione. 

5. Una misura invariante per i sistemi di Anosov 

5.1. Misure di Gibbs. Prima di tornare a discutere la questione dell'esistenza di una misura di proba- 
bilità invariante per un sistema di Anosov, diamo una definizione generale di misura di Gibbs; come vedremo, 
sia la misura invariante di un sistema di Anosov che la codifica in simboli della misura di volume apparter- 
ranno a questa classe di misure di probabilità. Prima di fare ciò, specifichiamo le condizioni "fisiche" che 
devono soddisfare i potenziali a molti corpi di una misura di Gibbs. 



Definizione 5.1. (Potenziali per la dinamica simbolica) 



Sia T una matrice di compatibilità mescolante per le sequenze in {0,...,n}^. Chiameremo B lo spazio 
delle sequenze $ — {"l'xlxcz parametrizzate dai sottoinsiemi finiti di Z, formato dalle funzioni 



(102) $x : {0,--,w}t 

che, ponendo ||$x|| = ^^^^^{o n}^ verificano le seguenti proprietà. 

(1) $ è invariante per traslazioni su Z, ovvero V cti, Up 

(103) $x(cri,...,crp) = ^rx{<yi,-,<yp) 

seX^ (Ci, Cp) e tX = (Ci + 1, Cp + !)■ 

(2) $ è "sommabile" : 

(104) ||$|| ^ 5] ^ < +00. 
Diremo che B è lo spazio dei potenziali su {0, ...,n}'^. La funzione 

(105) A,i^) = j: 

X30 ' ' 

sarà chiamata energia potenziale per sito oppure funzione energia associata a $. 

Osservazione 5.1. Possiamo prendere in prestito dalla meccanica statistica classica un semplice esempio 
di funzione energia, considerando ad esempio un sistema di spin interagenti su un reticolo bidimensionale 
con energia potenziale totale data da 

(106) E -.-.'J>(->^)-E E ^^^^^^4^; 

l'argomento della prima sommatoria è l'energia potenziale per sito, e axcry'^{x,y) = ^x,yW,'^') è un 
potenziale a due corpi. 

In questo contesto, possiamo dare una definizione generale di misura, o stato, di Gibbs. 

Definizione 5.2. (Misura di Gibbs) 

La misura m è una misura di Gibbs su {0, con potenziali $ = {$x}xez ^ Lì se è una misura di 

probabilità sui boreliani di {0, le cui probabilità condizionate verificano m-quasi ovunque la proprietà 



™(gA|gAc) 

^('ZaI'Zac) 

per og ni in tervallo connesso K <Z "L tale che |A| > a{T), se a{T) è il tempo di mescolamento di T (ci. 



definizione 3.1), e per — ((t^ct^c ) , ct" = (ct^ct^c) G {0, 



La proprietà (101), insieme alla condizione di compatibilità to((7^|(t^c) = 1, implica che 
(108) mi<j^\a^.) 



interpretando la (108) come zero se {0, ...,n}j 



Osservazione 5.2. (1) Il paragone del risultato (90) ottenuto per la misura di volume con la (101) ci 



permette di dire con cognizione di causa che la misura di volume è uno stato di Gihhs. 
(2) In altri approcci le misure di Gibbs vengono introdotte diversamente, ad esempio come soluzione 



di un problema variazionale; in quei casi la proprietà (101), che abbiamo assunto nel nostro caso 
essere la proprietà che definisce una misura di Gibbs, diventa un teorema. La (101) forma un 
insieme di equazioni per le probabilità condizionate che vengono chiamate equazioni di Dobrushin- 
Lanford-Ruelle, o più semplicemente equazioni DLR. 

5.2. La misura SRB. Concludiamo la nostra introduzione alla teoria dei sistemi di Anosov discutendo 
il problema dell'esistenza di misure di probabilità invarianti, anche dette statistiche del moto, per questa classe 
di sistemi dinamici; in generale detS' ^ 1, dunque la misura di volume non sarà invariante sotto l'azione 
dell'evoluzione temporale. Ciononostante, attribuiremo una particolare importanza alle statistiche del moto 
i cui dati iniziali siano scelti a caso con una misura di probabilità assolutamente continua rispetto al volume. 

Definizione 5.3. (Misura SRB) 

Sia (f2, S) un sistema di Anosov. Se no è una misura di probabilità assolutamente continua rispetto alla 
misura di volume su fi e per tutte le funzioni continue F definite in fi il limite lim^v-^oo ^J=o^ F{S^x) 
esiste iiQ-quasi ovunque e può essere scritto come 



N-l . 

(109) hm iV-i ^(^'^) = / l^idy)F{y), 

allora ii è chiamata la misura SRB di (fi, S). 

Osservazione 5.3. La proprietà di invarianza sotto S della misura SRB discende banalmente dalla sua 
definizione, infatti 

N-l 



I fi{dy)F{Sy) = lim N'^ V F{S^+'^:. 

In Af— i-oo ^ — ' 

N-l . 

(110) = hm iV-i J2 ^(^'^) = / l^idy)F{y), 

se F è limitata, e ciò è implicato dalla continuità di F e dalla compattezza di f2. 

Proposizione 5.1. (Esistenza della misura SRB) 

Sia (fi, S) un sistema di Anosov. Esistono due misure di probabilità fi^, tali che, per ogni F £ C^ifl), 
a> 0, 

N-l „ 

(111) hm nS^'^) = / l^±idy)Fiy) 

3=0 

quasi ovunque in x rispetto alla misura di volume /io su fi. 

Osservazione 5.4. (1) Notare che, a differenza del teorema di Birkhoff, cf. proposizione 2.2. la misura 

di probabilità con la quale vengono scelti i dati iniziali non è invariante. 
(2) In generale, la misura SRB nel futuro /i+ e la misura SRB nel passato /i_ sono diverse, ovvero per 
una generica funzione Haider continua F fi^{F) ^ fi^{F). 

Dimostrazione. Per una dimostrazione di questo teorema rimandiamo a [8]; la dimostrazione non è 
particolarmente difficile, ma si basa su alcuni risultati che, dato il carattere introduttivo di questi primi due 
capitoli, non abbiamo discusso. □ 



Dalla dimostrazione della proposizione |5.l| scopriamo che la misura SRB è uno stato di Gibbs determinato 
dai potenziali a molti corpi generati dallo sviluppo del tasso di espansione yl„(tT), per la statistica nel futuro, 

da Ag^a), per la statistica nel passato; ciò equivale a dire che, ad esempio nel futuro, le probabilità 
condizionate soddisfano le equazioni DLR 

(112) iiffey^-'^ClL -""''''^'')- 

dove a' = {(z'j^oi^c) , <z" = Ì^a^a^)- Questo fatto può essere capito facilmente in modo euristico, cf. 
appendice |Xj 

E piuttosto curioso notare come il problema apparentemente senza speranza della statistica dei moti di 
un sistema estremamente caotico, in un certo senso il paradigma di un sistema caotico, possa essere risolto 
attraverso lo studio di un modello di Ising unidimensionale, uno dei piii semplici modelli noti in meccanica 
statistica; in particolare, nelle appendici \C\ [P] dimostreremo alcune proprietà di analiticità della misura SRB 
attraverso una tecnica ben nota in meccanica statistica dei sistemi di spin. 

6. Un'applicazione: il gatto di Arnold perturbato 

Come abbiamo già visto nella sezione |3.1[ gatto di Arnold è un semplice esempio di sistema di Anosov; 
l'evoluzione di questo sistema avviene sul toro con la mappa 

(113) 5=(; i). 

La misura SRB di questo sistema di Anosov è banalmente la misura di volume, dal momento che 

1 detS"! = 1. Come vedremo, perturbando "poco" la mappa S il sistema perturbato sarà ancora un sistema 
di Anosov; in particolare, potremo calcolare esplicitamente i potenziali della nuova misura SRB. Tutta la 
nostra analisi sarà resa possibile dal seguente teorema, dovuto ad Anosov. 

Proposizione 6.1. (Teorema di stabilità strutturale) 

Siano (17, S) e {Q, S') due sistemi di Anosov. Se le mappe S , S' sono abbastanza vicine insieme alle loro 
derivate prime allora esiste un omeomorfismo H : Vi ^ Q tale che 



(114) SoH = HoS'. 

Osservazione 6.1. In altre parole, ciò vuol dire che è possibile mappare una dinamica nell'altra attraverso 
un "cambiamento di coordinate", cf. figura [9} generalmente non differenziabile. 

Tra poco costruiremo esplicitamente H nel caso del gatto di Arnold perturbato, e dimostreremo che H è 
analitica nella perturbazione e solo Hòlder continuo nel punto. Per una "idea" della dimostrazione generale 
di questo teorema, ovvero che non si basi su un caso specifico, rimandiamo a |10j . Nel nostro caso, seguiremo 
l'analisi svolta in [8], |12j . 

6.1. Costruzione della coniugazione H. Ci proponiamo di calcolare esplicitamente la coniugazione 



H partendo dalla proprietà (114), nel caso in cui l'evoluzione del sistema perturbato sia descritta da 



(115) = ^0^ + £/(^) 



dove 5*0 è data dalla (1131 e f{ip) è un polinomio trigonometrico, ovvero /(V') = SìyGZ2 \u\<n ^^^'^f ^'^ 

funzione analitica in ?/;• Per fare ciò, cercheremo una soluzione dell'equazione 



assumeremo che fj^ip) sia una 




Figura 9. Rappresentazione schematica dell'azione della coniugazione H ("diagramma 
commutativo"). Per semplicità, nella figura abbiamo disegnato una traiettoria continua. 



(116) 



HoS^S.oH 



che sia esprimibile come sviluppo in serie di potenze in e, ovvero H(ip) = tp + X]fc>i ^^hii'): e dimostre- 
remo che la serie è convergente per |e| < Eq > 0; quindi, in questo intervallo la coniugazione H sarà una 
funzione analitica di e. Inoltre, con la stessa tecnica otterremo anche l'Hòlder continuità in con esponente 
/? per e < £(/3) < £(0) — eo. Scriviamo la (115) come 



(117) 



ho So = Sq o h + ef o [Id + h); 



al primo ordine in e la (117) diventa 



(118) 



So o h'^^^ + f = h'-^K So. 



Chiamiamo v^, v_ i due autovettori normalizzati di So corrispondenti rispettivamente all'autovalore A+ 
e all'autovalore A_, e definiamo A = A^^ — A_; poiché 



(119) 



h{il:) = A^[iP)v^ + h_{il:)v_, 



le due componenti della (118) nelle direzioni i;_ sono 



(120) 
(121) 



Possiamo risolvere le equazioni ( 121 ) iterativamente, e il risultato è 



(122) hi'Hf) = -a ^ \-\^+'\'^USP,l) a = ±, 

dove Z4. = [0, 00) n Z e Z_ = (—00, 0) n Z. In generale, all'ordine fc-esimo si avrà che 



E 

ai, ...,as = ± 



ki + ... + ks ^ k- l,k,> l 



(123) 



dunque è chiaro che tutti gli ordini dello sviluppo di h possono essere costruiti in funzione di h^^^ ; inoltre. 



l'espressione ( 123 1 può essere rappresentata in termine di alberi, cf. figura 10 




Figura 10. Rappresentazione grafica della (123). 



La parte sinistra della figura 10 corrisponde a h^\ip), i rami della parte destra sono h)^^' , h)^^' ecc., 
mentre il vertice v equivale a 



(fel) I,(fe2) 




(124) ^aA-l^'+H" \^j[v,^ ■ d^j USl±)- 

le somme su ki, a, p, s sono sottintese. Un albero d con k nodi porta su un generico ramo / una coppia 
di indici ai, pi con e Z e a; e {— , +}, e sui nodi v due simboli a^, con Uy = ai^ e Pi, e Z^, tali che 



(125) p{v) =pi^_ = ^Pw, 

dove la somma è sui nodi che seguono v (l'ordinamento è specificato dalle frecce), e ly denota il ramo 
v'v uscente dal nodo f ; ad ogni albero assegniamo un valore 



(126) 



Val(^) = n n^-. f'^ASrf)^ 



dove da = Va' ^ l'insieme di nodi in e vi, sono gli s„ nodi che precedoncr^w. Se con 

Qk,a indichiamo l'insieme di tutti gli alberi con k nodi e con etichetta a associata alla radzcedell'albero, 
ovvero al ramo uscente dall'ultimo vertice, allora la quantità hai'ip) può essere rappresentata come 



(127) 



/i„(7A) = ^e^- J2 Val(^9), 

k=i tìG0k,c 



e l'ultima cosa che rimane da fare è stimare il raggio di convergenza della serie di potenze. Per fare ciò 



useremo la stima di Cauchy, grazie alla quale possiamo maggiorare le derivate nella ( 126 ) come, chiamando 



G il massimo di fa nella regione V = {ipi : llm-^il < rg, i = 1, 2} per un certo rg > e m-|_, m_ le molteplicità 
delle derivate 9+, 9_ 



(128) 



'77? ' 



< G- 



dal momento che mi + m2 = Sy; grazie alla (128) abbiamo chep^ 
(129) 



ivai(z?)i< n n ^i^"^""" 



v<£V{-a) 



vev(-d) 



dunque poiché il numero di possibili alberi con k nodi è < 2 



2fe 



(130) 



Val(?9) 



< 2^2 



fcr,2fc 



„fe-i 



1 - A 



Quindi, grazie alle (127), (130) il raggio di convergenza di è dato da 



33> 



Nel caso in cui v sia un nodo iniziale, i.e. in cui corrisponda ad una "foglia" dell'albero, nella 1 126 1 non c'è nessuna 

derivata 

■^'^Per ottenere l'ultima uguaglianza nella 1 129 1 abbiamo usato che J^uge Sk = k — 1. 



(131) 



G 1 



rol-AJ 

L'Hòlder continuità di H può essere dimostrata semplicemente notando cha^ per < /3 < 1, 



(132) 



infatti, la ( 132 1 implica che 



(133) 



\4' ~ V'' 

Quindi, per |e| < e(/3), con 



/3 



< 2^''2'' 



1 - Ai-zs 



2^^ 



(134) 



£(/?) 



1 



ro 1 - Al- 

la coniugazione H è Hòlder continua in i]} con esponente /? e analitica in e, dal momento che e(/3) < 
e(0) = £o- Notare che 77, in generale, non è lipschitziana, perché lim^_>.i e(/3) = 0. 

Infine, la mappa iJ è un omeomorfismo. Infatti, se ^"1 7^ "02 allora H^ipi) ^ H{^p2), perché ciò impli- 
cherebbe per ogni k S'^(iJ(V'i)) = S^{H{ìp2)) ^ H{Sq^i) = H{SQÌp2); ma ciò è assurdo, dal momento che 
Sq è mescolante e H è Hòlder continua. Inoltre per ogni ip esiste ip tale che H{'ij;) = ip; questo è facile da 
visualizzare sul toro T^, dove se ^p "gira" attorno al cerchio anche il punto H{'ijj) lo farà, e dunque per ogni 
fissato '0 esiste e una funzione continua k{ip) tali che = (^ + fc(<y9). Per d generico l'argomento è lo stesso, 
solo ripetuto d volte. 

6.2. Costruzione della misura SRB. Il prossimo passo consiste nel riuscire a costruire esplicitamente 
i potenziali di Gibbs della misura SRB del sistema perturbato; per fare ciò, abbiamo bisogno di calcolare 
esplicitamente il coeSiciente di espansione, se siamo interessati alla statistica nel futuro, oppure il coefficiente 
di contrazione, se vogliamo descrivere la statistica del sistema nel passato. Ingenuamente, potremmo pensare 
di parametrizzare la varietà stabile e la varietà instabile del sistema perturbato attraverso la coniugazione 
H, ovvero con delle equazioni tipo 



(135) (£(É) 

dal momento che, per a = ±, i — ^ - 
è destinato a rivelarsi infruttuoso, perché la parametrizzazione definita tramite H non può che essere in 
generale solo Hòlder continua in i, e dunque inutile ai fini della costruzione dei coefficienti di contrazione e 



tv^ è la varietà instabile o stabile di 5*0. Questo tentativo però 



di espansione, c/. definizione 4.2 

Chiamando w_^{1l>)^ W-l'ip) i vettori tangenti alla varietà instabile e alla varietà stabile del sistema 
perturbato nel punto ■(/;, questi vettori devono soddisfare l'equazione 



(136) 



■^^Per ottenere la ( 132 I abbiamo scritto la differenza nella parte sinistra come l'integrale di una derivata, abbiamo stimato 
la derivata con la formula di Cauchy e abbiamo usato il fatto che la massima distanza tra due punti sul toro è \/2tt. 



dove A± è in generale una funzione Hòlder continua di ■0, ed è ciò di cui avremmo bisogno per costruire 
i potenziali della misura SRB. 

Però, per delle ragioni che saranno chiare tra poco, saremo in grado di costruire esplicitamente queste 



quantità calcolate nel punto e non in ip; l'equazione ( 136 1 diventa, sostituendo ìp con H{ìp) e definendo 

= o H , X± — X± o H, 

(137) DS,{H{f))w^{Ìj) = X±{f)w^{Sof). 

Come nel caso della coniugazione H, possiamo provare a risolvere la ( |137[ ) perturbativamente, ovvero 
cercando una soluzione del tipo A±(V') = A± + SX±{ip), Wj.{il^) = Vj. + Swj.{Tp) con 

(138) SX±{f) = 7±(É) = E^±^W 

fe>i 

(139) Sw^if) = k^{t)v^ = J2k^^\t)v^, 



k>l 



(fe) , (k) 



O (e*')- Scrivendo la (137) nella base degli autovettori del gatto di Arnold Vj., ponendo 



DSei^j) ^So + eDfitP), la (137) implicache 



(140) 

dove a 



-X {edpf^i^^) + ek^{f')dM^") - k^iSot^Mt")) , 



±e /3 



a. Le equazioni ( 140 ) possono essere risolte ad ogni ordine in e in funzione di tutti 
gli ordini precedenti, e dunque alla fine del primo ordine, dato dsrj 



(142) 
(143) 



v(l) 



fc«(7/^) 



-(2n+l) 



(144) 



k^'\^) 



A;('"+'^a+/_(5o"("+'V)- 



Con una tecnica analoga a quella usata nel caso della coniugazione H si può dimostrare che esistono 
Eg, £'(/?) tali che gii sviluppi in e di 7±, k± convergono per |e| < Sg, e che le soluzioni 7±, k± sono Hòlder 
continue con esponente /3 in |e| < e'(/3) < e'(0) = eó, c/. [8], |12j . 

Grazie questi sviluppi perturbativi siamo ora in grado di dire esplicitamente chi sono i coefficienti di 
espansione e di contrazione e i vettori tangenti alle varietà stabili e instabili nel punto H{'ip). Se non avessimo 



composto sin dall'inizio con la coniugazione H avremmo trovato la (1401 con 5*0 sostituito da S^; dunque. 



per dimostrare l'analiticità in e avremmo dovuto prima verificare la differenziabilità in ijj, e viceversa. Ad 
ogni modo, il fatto di essere in grado di calcolare le quantità necessarie per la costruzione della misura SRB 



■^^Notare che se al posto del gatto di Arnold (113 i consideriamo la sua "radice quadrata' 



(141) 



1 1 
1 



i fattori di convergenza _\ saranno sostituiti da X,^^"'^^\ dal momento che in questo caso A_ = —A, 



solo nel punto H{'ìp) non rappresenta un problema, perché se Vq è una partizione di Markov per il sistema 
imperturbato (fi, So) possiamo introdurre una partizione di Markov Ve per (il, S^) come 



(145) Ve = H{Vo); 

dunque, può essere rappresentato con la stessa storia a di solo codificata con un nuovo codice 

Xe = H o X. Quindi, possiamo scrivere i valori medi SRB come delle somme sulle storie simboliche del 
sistema imperturbato codificate con il nuovo codice Xe, pesando ogni addendo con un fattore 

(146) , 

dove Au = Au o Xe è dato da 



(147) 



Au{ra) = log f [A+ + i+iXeia))] ^ 
\ I— - 



w4Xe{Ta))\ 
{Xe{a))\ 



Riassumiamo il contenuto di questa sezione nella seguente proposizione 



Proposizione 6.2. (Potenziali del gatto di Arnold perturbato) 

Data una partizione di Markov Vo = {Pi, Pn} in per So, sia a la rappresentazione simbolica 
rispetto a Ve = H{Vo) di un punto ip ~ Xe{a). Esiste s'{/3) tale che il coefficiente di espansione Xu{o^) di 
Se{g^) lungo la varietà instabile di ip è analitico in e nel disco \e\ < £'{P), e come funzione di a è Haider 
continua con esponente j3 e modulo C{(3); lo stesso vale per il tasso di espansione 



(148) A^ira) = log [A+ + j+{Xe{<7))] 



\w+iXe{Ta))\ 

k+(^e(a))| 



che genera (nel senso discusso nella proposizione 4.1) i potenziali (fl[-ri n]) dello stato di Gibbs 

associato alla misura SRB. 



Osservazione 6.2. (1) Dal momento che tutte le funzioni che appaiono nella (14^8) sono Haider con- 

tinue i potenziali '^x{q_x)' con X = [—n,n], decadono esponenzialmente, i.e. esistono C, c > tali 
che, se o![-n n] corrisponde ad una sequenza a "troncata" nell'intervallo [-~n,n\ e continuata con 
opportune condizioni al bordo, 

(149) |$h„,«](fl[-„,„])| = l^n(^Z[-(„+i).(„+i)]) - A.(ct[_„,„])| < Ce-^"; 

nella teoria dei .sistemi di spin unidimensionali, un potenziale di questo tipo viene chiamato un 
potenziale di Fisher. 

(2) Notare che, nel limite termodinamico, il rapporto \ ^!"|^ non modifica lo stato di Gibbs, perché 

\'^+ (^l)) I 

sia il numeratore che il denominatore sono funzioni Haider continue; dunque 

A/2-1 

(150) {W+iXeir'^'m - \w4Xe{T^a)\) ^ 0(1), 

fc=-A/2 



mentre X]fc=-A^/2 ^^(^'^— ) ~ ^i^)' P^fciò la (150) rappresenta un termine di bordo. Ciò vuol 
dire che lo stato di Gibbs del sistema è indipendente dalla particolare metrica scelta per misurare 



i moduli dei vettori, a patto che si tratti di scelte "ragionevoli" , ovvero che gli elementi di matrice 
del tensore metrico siano funzioni Haider continue. 



CAPITOLO 3 



Non equilibrio 

Il contatto tra la teoria matematica dei sistemi di Anosov che abbiamo appena discusso e la meccanica 
statistica dei sistemi fuori dall'equilibrio è contenuto nell'ipotesi caotica, che presenteremo nella sezione [2] 
prima di fare ciò però, discutiamo brevemente il significato dell'ipotesi ergodica e le sue conseguenze nella 
descrizione statistica dei sistemi all'equilibrio, 



1. L'ipotesi ergodica 

La meccanica statistica dell'equilibrio ha come scopo la deduzione delle proprietà macroscopiche di in- 
siemi di particelle che evolvono seguendo una dinamica hamiltoniana; ciò vuol dire che se H è la funzione 
hamiltoniana e p = (pi, ...,p3jv), q = {qi, ...,(73Ar) sono rispettivamente gli impulsi e le coordinate spaziali di 
N particelle, allora il sistema evolverà risolvendo le equazioni differenziali 

(151) Pi = ^~^^ 

dq, dpi 

Nella realtà fisica l'impulso e la posizione non possono essere misurati simultaneamente con precisione 
arbitraria; una misura "molto precisa" dell'impulso implicherà una misura "poco precisa" della posizione, 
nel senso che se indichiamo con Sp, 6q le incertezze su impulso e posizione allora la massima precisione 
raggiungibile sarà limitata dalla condizione 

(152) SpSq = h, 
dove h è un'opportuna costante. In generale avremo che 

(f53) SpSq > h, 



e, se h coincide con la costante di Planck, l'affermazione ( 153 ) è il contenuto del principio di indetermi- 
nazione di Heisenberg. 

Quindi, possiamo considerare lo spazio delle fasi come decomposto in tante piccole cellette numerate 
Ai, con i — l,...,Af, ognuna di volume h; ogni celletta contiene un continuo di possibili configurazioni 
microscopiche {p,q), ma tutte le configurazioni contenute all'interno di ogni si riferiscono allo stesso stato 
fisico, nel senso che le proprietà macroscopiche del sistema non cambiano in modo fisicamente apprezzabile 
se ci si sposta all'interno di una celletta. Nonostante il numero di possibili configurazioni microscopiche sia 
infinito, lo stato fisico della totalità delle particelle può variare solo in un numero finito di modi; ad ogni 
modo, per un sistema di 10^"^ particelle (e anche per sistemi molto piìi piccoli) questo numero è enormemente 
pili grande di qualunque scala astronomica. 

Supponiamo che l'hamiltoniana totale del sistema non dipenda esplicitamente dal tempo, ovvero §iH ~ 
0; ciò implica che l'energia totale sia uguale ad una costante, che chiameremo E. Dunque, possiamo pensare 



che la dinamica hamiltoniana ( 151 ) induca una mappa S sulla superficie dello spazio delle fasi ad energia 



costante la cui azione consista in una permutazione delle cellette A^; chiaramente, questa mappa rifletterà 



43 



tutte le proprietà della dinamica hamiltoniana, in particolare, per il teorema di Liouville, il volume delle 
cellette non varia sotto S. 

Nonostante questa schematizzazione, il problema della deduzione di proprietà macroscopiche di un siste- 
ma - e dunque della sua termodinamica - a partire dalla dinamica hamiltoniana rimane ancora inaccessibile; 
per riuscire nel nostro scopo abbiamo bisogno di sapere, o di assumere, qualcosa in più sulla natura della 
dinamica microscopica. L'idea di Boltzmann, tradotta nel nostro linguaggio, è la seguente. 

Ipotesi ergodica (lE): l'azione della trasformazione di evoluzione S, vista come permuta- 
zione di celle dello spazio delle fasi sulla superficie con energia costante, è una permuta- 
zione a ciclo unico; 



(154) SAk = Ak+i k^l,2,...,M 

se le celle sono opportunamente numerate (e Ajv'+i = Ai). 

L'ipotesi ergodica rende il più semplice possibile l'azione di S sulla superficie dello spazio delle fasi ad 
energia costantqj come vedremo in modo euristico, la conseguenza di questa ipotesi è che le medie temporali 
possono essere calcolate come integrali di volume sullo spazio delle fasi. 

Chiamiamo C l'insieme di cellette A che appartengono alla superficie ad energia costante, J\f il loro 
numero totale, e definiamo la misura di probabilità ^ come 



(155) 



m(A) 



1/Af se A e C 
altrimenti 



chiaramente fi è stazionaria, ovvero /x(5'A) = /x(A), dal momento che S permuta celle dello stesso ciclo. 
Dunque, assumendo (lE) per il teorema di Birkhoff (c/. proposizione 2.2) abbiamo che, scegliendo un "dato 
iniziale" A, 



1 _ 

(156) ^lim -^0(5^-A)=^0(Afe)M(Afc); 

inoltre, dal momento che h è piccolo possiamo approssimare la somma sulle celle con un integrale, 



1 .-^ L ,0(p.q)d^^pd^'^q 
(157) lim ^YO (S'^A) = ^-^^^ -, 

dove Z = d^^ pd^^ q è la normalizzazione. Quindi, l'ipotesi ergodica attribuisce una parti- 

colare importanza alle misure di probabilità assolutamente continue rispetto al volume, ovvero del tipo 



H (d'^'^ pd^^ q) = f{p,q)d^^pd'^'^q; l'integrale (1571, ad esempio, può essere riscritto come un integrale su 
tutto lo spazio delle fasi fi con una misura fj, ((F>rf^^g) = 6 {n = E) d^^pd^^q, se <5 è la delta di Dirac. 



In generale, per un sistema fisico l'integrale nella ( |157 ) è estremamente difficile da calcolare; ad ogni 



modo, ciò che è davvero fisicamente rilevante sono le relazioni tra valori medi di osservabili macroscopiche 



piuttosto che i valori medi in sé, e molto spesso queste relazioni possono essere dimostrate partendo dalla ( 157 ) 
senza dover necessariamente calcolare gli integrali. Quindi, accettando l'ipotesi ergodica la termodinamica 
di un sistema con un numero arbitrario di gradi di libertà diventa il riflesso della dinamica hamiltoniana, o 
meglio di poche sue proprietà. 



Chiaramente, questa ipotesi è manifestamente falsa nel caso in cui il sistema ammetta altre costanti del moto oltre 
all'energia totale; ciò implica che la dinamica, vista come permutazione di celle, non può essere un ciclo unico. In questi 
casi l'ipotesi va riformulata considerando l'azione di S come permutazione delle cellette sulla superficie dello spazio delle fasi 
determinata da H = E e da tutti gli integrali primi della dinamica. 



2. L'ipotesi caotica 



I sistemi fuori dall'equilibrio rappresentano una classe di sistemi di particelle per i quali l'ipotesi ergodica 
non fornisce le risposte giuste. Come abbiamo discusso nell'introduzione, questi sistemi sono sottoposti a 
delle forze non conservative, ovvero dissipano calore. In particolare, poiché la dinamica non è hamiltoniana 
la divergenza delle equazioni del moto sarà in generale diversa da zero; ciò rende il teorema di Liouville inap- 
plicabile, e con esso anche l'ipotesi ergodica. Per avere una comprensione microscopica della fenomenologia 
di un sistema fuori dall'equilibrio sembra essere indispensabile fare delle nuove assunzioni sulla dinamica 
delle particelle, necessariamente diverse dall'ipotesi ergodica, e verificare che queste ipotesi abbiano delle 
conseguenze osservabili. 

II fatto che il volume di una celletta A non sia conservato dall'evoluzione temporale implica che se il 
sistema è confinato in uno spazio delle fasi limitato la contrazione dello spazio delle fasi non potrà che essere 
in media minore o uguale di zero; nel caso in cui sia strettamente minore di zero il sistema, nel limite di 
tempi infiniti, evolverà in un insieme di misura di volume nulla. Dunque, ammesso che esista la misura di 
probabilità stazionaria che descrive lo stato del sistema non può essere assolutamente continua rispetto al 
volume. 

Nel capitolo [2] abbiamo visto che i sistemi di Anosov (c/. definizione 3.3) ammettono una misura di 
probabilità invariante - chiamata misura SRB - che ha dei potenziali diversi da quelli della misura di volume, 
dunque in generale è sm(?oZor^rispetto a quest'ultima; ci chiediamo se questi sistemi dinamici, oltre ad essere 
interessanti da un punto di vista matematico, abbiano qualcosa a che vedere con la Fisica. A prima vista 
non sembrerebbe: un sistema reale, ad esempio un gas composto da 10^^ particelle, è probabilmente caotico 
in qualche senso, ma è decisamente troppo restrittivo pensare che la dinamica sia strettamente assimilabile 
a quella di un sistema di Anosov. 

Anche l'ipotesi ergodica però, se presa alla lettera, è praticamente impossibile da verificare nelle appli- 
cazioni, e gli unici sistemi per i quali può essere dimostrata rigorosamente sono estremamente "semplici" , 
tanto da non avere nulla a che vedere con la realtà fisica studiata dagli sperimentatori. Ciononostante, le 
conseguenze dell'ipotesi ergodica, ovvero le relazioni tra valori medi di osservabili che formano la termo- 
dinamica, sono largamente verificate. Ciò suggerisce che non abbia senso, o perlomeno che sia superfiuo, 
verificare la stretta validità dell'ipotesi ergodica in ogni sistema fisico all'equilibrio che si vuole studiare con 
le tecniche della meccanica statistica; la cosa che ha senso fare, invece, è accettare che al fine di calcolare 
quantità macroscopiche il sistema si comporti come se fosse ergodico. 

Il nostro punto di vista sarà esattamente lo stesso; assumeremo che su scala macroscopica il fatto che 
il sistema non sia strettamente un sistema di Anosov non abbia nessuna importanza. Riassumiamo queste 
idee con la seguente ipotesi caotica, [4], [3], |13) . |14) . 

Ipotesi caotica (IC): Allo scopo di misurare proprietà macroscopiche, un sistema di parti- 
celle in uno stato stazionario si comporta come un sistema di Anosov reversibile. 

In (IC), per sistema reversibile intendiamo un sistema dinamico (fl,S) per il quale esiste una mappa 
(differenziabile) / - un' inversione temporale - tale che 



(158) IoS = S-'^oI con/2 = ±l; 

vedremo in seguito che la reversibilità giocherà un ruolo importante per la deduzione di conseguenze 
osservabili dell'ipotesi caotica. 

L'ipotesi caotica è un'assunzione che prescinde dal fatto che il sistema si trovi o meno all'equilibrio. Un 
sistema di Anosov conservativo, ovvero per il quale la contrazione dello spazio delle fasi è identicamente nulla. 



Due misure fi, u definite in fi sono singolari se esistono due insiemi disgiunti A, B £ fi con AU B = fi tali che fi è zero 
su tutti i sottoinsiemi misurabili di B mentre u è zero su tutti i sottoinsiemi misurabili di A. 



ammette come misura di probabilità invariante la misura di volume; dunque, l'ipotesi caotica all'equilibrio 
implica l'ipotesi ergodica. 

La validità di (IC) può essere decretata solo ed esclusivamente partendo da risultati sperimentali, e 
per avere qualcosa da verificare sperimentalmente dobbiamo dimostrare che (IC) abbia delle conseguenze 
osservabili. Il teorema di fluttuazione di Gallavotti - Cohen è un esempio (per il momento uno dei pochi) 
di conseguenza osservabile dell'ipotesi caotica; come abbiamo accennato nell'introduzione, questo teorema 
stabilisce una proprietà di simmetria sulle fluttuazioni del tasso di produzione di entropia. Ma, per il 
momento, non è chiaro cosa siano l'entropia e il suo tasso di produzione in un sistema fuori dall'equilibrio; 
discuteremo questo punto nella prossima sezione. 

3. Il problema della definizione di entropia 

In generale, un sistema fuori dall'equilibrio in uno stato stazionario dissipa calore ad un regime costante; 
dunque il sistema produce e cede entropia, nel senso che il valor medio del tasso di produzione di entropia 
è strettamente positivo e indipendente dal tempo. Quindi, l'entropia totale del sistema non può essere una 
quantità ben definita. 

Nelle prossime due sezioni discuteremo l'identificazione del tasso di produzione di entropia con il tasso di 
contrazione dello spazio delle fasi; ci baseremo prima su un particolare modello di sistema fuori dall'equilibrio 
molto importante nelle applicazioni, [2], |15j . |16| . e in seguito faremo vedere che, se la misura di volume 
ammette un limite debole per tempi infiniti, questa identificazione è del tutto generale, |17j . |13j . 

3.1. Un modello. Indichiamo le posizioni delle particelle all'interno di un contenitore Cq con x = 
{xi, ...,Xj^j), chiamiamo V{xi, Xpf) — V{x) l'energia potenziale d'interazione e F^{x;G), i — 1,...,N, le 
forze esterne non conservative, dove G — (Gi, Gs) sono i parametri che regolano le intensità delle F^; dun- 
que, F{x] 0) = 0. Come abbiamo accennato in precedenza, per evitare che l'energia aumenti indefinitamente, 
ovvero per garantire l'esistenza di uno stato stazionario, è necessario che il sistema sia posto a contatto con 
un termostato; quindi, se {— sono le forze dovute alla presenza del termostato, i — 1, ■■■,N, le equazioni 
del moto delle particelle all'interno di Cq possono essere scritte come, indicando con le masse delle N 
particelle. 



(159) m,x, = -d^y{x) + F^ix;G)-d,, i = l,...,N. 

Si può pensare ad un tipo di termostato come ad un insieme di serbatoi infiniti, ad esempio contenenti un 
gas, messi a contatto con il sistema S e asintoticamente in equilibrio termico; ogni serbatoio, che chiameremo 
Ca, dove a = l,2....,m, sarà caratterizzato da una temperatura Ta = cost. Questi termostati, però, sono 
ovviamente intrattabili da un punto di vista sperimantale, i.e. numerico. Dunque, possiamo immaginare 
un'altra classe di termostati pensando che il sistema Co sia a contatto con un numero m di serbatoi finiti 
Ca, e che le forze non conservative agiscano solo su Cq; sui sistemi {Sa}, contenuti nei {Ca}, faremo agire 
delle forze {'i?a} in modo da imporre il vincolo anolonomo 



(160) Ka^Y.1 

3 

Definiamo la temperatura del sistema Ea come 



^maXj — cost. 

3 



iNaKB 

dove fcs è la costante di Boltzmann e ÌS!a è il numero di particelle contenute in Ca \ accetando la definizione 



(1611, le forze {i^a} mantengono i sistemi a temperatura costante. Notare che, a rigore, ogni modello di 
termostato finito deve essere considerato "non fisico", perché in generale è impossibile separare il sistema 



dal resto dell'universo; ma d'altro canto non ci rimane scelta, se vogliamo simulare l'evoluzione di questo 

modello su un calcolatore. Ad ogni modo, è un'opinione diffusa che, nel limite termodinamico, il modo in cui 

venga estratto calore dal sistema sia inessenziale, ovvero che lo stato stazionario sia indipendente dal tipo di 

termostato utilizzato; il fatto che lo stato fisico di un sistema sia indipendente dal tipo di termostato scelto 

darebbe luogo ad una "equivalenza degli ensemble" per la meccanica statistica fuori dall'equilibrio, |18) . 

Se le particelle contenute in Co interagiscono mediante un potenziale ip{q. — q .), dove q , q . rappresentano 

—i —j —% —j 

le coordinate spaziali di due particelle, l'energia potenziale d'interazione del sistema Sq sarà data da 



(162) 



analogamente, indicando le coordinate delle particelle contenute in Cq con j ^ 1, Na, si ha che 



(163) 



K(x)=^^„(x^-x,") 



mentre l'energia potenziale dovuta all'interazione tra le particelle all'interno di e quelle contenute in 
Co potrà essere espressa come 



N 



(164) 



Waiq,x'')^Y.T.'^-k^-^'ì^' 
i—1 j<i 



dove Wa{q^ — Xj) rappresenta l'interazione tra due particelle contenute rispettivamente in Ca e in Co 

nelle posizioni x^, q .. Infine, ipotizziamo che i sistemi {Sa} non interagiscano direttamente tra di loro, ma 

—j 

solo attraverso Co- 



Grazie alle (162), (163), ( |164[ ), le equazioni del moto del sistema possono essere scritte come (assumiamo 
per semplicità che sia in E che nei {Ea} le masse siano tutte uguali): 



(165) 
(166) 



mq . 



= -dA ^0 (9) + E a;'^) + (9; G) 



^a^j = -d^^{Va{x'') + Waiq,xf))-l';, 



e, per realizzare la condizione ( 161 ), imponiamo che le forze {^j} verifichino il principio di minimo sforzo 
di Gauss, [13^ . Definiamo la forza e lo sforzo Za relativo alle accelerazioni {a^} come 



(167) 
(168) 



Bt = -~di^{Vaix) + Waiq,xf)) 



si può dimostrare che, cf. |13j . la forza da aggiungere alle che minimizza lo sforzo Za tra tutte le 

possibili accelerazioni {a°} compatibili con un vincolo 'ipa{x,x) = è data da 



(169) 



■dUt,x) = d±éaix,x) 



T,j % ■ d^,Ì'a{x,x) + :^B°j ■ d±^1pa{x,x) 



Quindi, se ipa — Ka la ( 169 1 diventa 



e La, Va rappresentano, rispettivamente, il lavoro esercitato dalle particelle di E su quelle appartenenti 
a Sa e una derivata totale, ovvero 



(171) 



(172) 



La = -^%iya(g,a;°)-i| 



N„ 



Va 



La quantità di calore prodotta per unità di tempo Q da un sistema allo stato stazionario è naturalmente 
definita come il lavoro che le forze del termostato esercitano sul sistema S per unità di tempo, ovvero 



(173) 

dal momento che 



(174) 

la quantità Q è data da 

(175) Q = Y.Q- 

a 

Introduciamo la contrazione dello spazio delle fasi dovuta al termostato a-esimo Zla(^, x) come 

(176) aa (i, x)=J2 % ■ ' 

j 

come vedremo, è ragionevole identificare questa quantità con il tasso di produzione di entropia nell'a- 



esimo termostato. La divergenza (1761 Ca = {5Na ~ l)a°- (nella (1761, z?^, è data dalla (170)) può essere 
scritta esplicitamente come 



(177) 



Va 



ksTa ksTa 



{iNg - 1) 



dove La corrisponde al calore ceduto per unità di tempo da S a Eq, e è una derivata temporale totale; 
assumendo che Va sia limitata, V non contribuisce in valor medicj^ 

Quindi, ponendo La = Qa possiamo definire una nuova contrazione dello spazio delle fasi come 



^Infatti, 



(178) 



lim — 



TJo 



dtV 



lim 

T->oo 



Vt-Vo 



se y è limitata. Ad ogni modo, è possibile trattare anche il caso in cui V non è limitata, cf. |19| . 



(179) 



dove stiamo trascurando di sommare una derivata totale e dei termini O (^^j- Chiaramente, l'espres- 
sione nella ( 179 1 è diversa dalla vera contrazione dello spazio delle fasi a; ma, per Na "grandi" , in valor 



medio la differenza è trascurabile. 



Per questa classe di termostati la ( 179 ) giustifica l'identificazione tra contrazione dello spazio delle fasi ed 



entropia prodotta per unità di tempo. Notare che non stiamo cercando in nessun modo di definire l'entropia 



per un sistema dinamico fuori dall'equilibrio: la (179) è uguale all'entropia ceduta per unità di tempo al 
termostato, che è all'equilibrio. 

Oltre a mettere in luce una relazione tra contrazione dello spazio delle fasi e tasso di produzione di entro- 



pia, i termostati gaussiani isoemetici rendono le equazioni del moto (1651, (1661 invarianti sotto V inversione 
temporale i — )■ —x, dal momento che = a^i^ è pari sotto questa trasformazione, e sia le forze che i 



potenziali d'interazione sono puramente posizionali; inoltre, per le definizioni (179), (171), la quantità (179) 



cambia segno sotto l'azione deW inversione temporale, naturalmente identificata come l'inversione dei segni 
delle velocità di tutte le particelle. Come vedremo in seguito, la reversibilità giocherà un ruolo importante 
nella deduzione di conseguenze osservabili dell'ipotesi caotica. 



3.2. Un'interpretazione generale. Se il sistema tende ad uno stato stazionario l'identificazione tra 
tasso di contrazione dello spazio delle fasi e tasso di produzione di entropia può essere dedotta in modo del 
tutto generale, i.e. senza senza riferirsi ad un modello particolare, |17j . Consideriamo le equazioni del moto 



(180) 



X = f{x), 



le quali inducono una mappa St sullo spazio delle fasi, e immaginiamo che al tempo t = lo stato del 
sistema sia descritto da una misura fio(dx) assolutamente continua rispetto al volume, ovvero della forma 
fj.o{dx) — po{x)dx; al tempo t la misura evolverà in Utidx) = pt{x)dx 



(181) 



Pt{x) 



dx 



Po{x), 



dove Jt{x) 



dx 



Gibbs al tempo t come 



è il determinante jacobiano in x della mappa x — S-tx. Definiamo l'entropia di 



(182) 



che diventa, per la ( 181 ), 



^(t) = - / Pt{x) log pt{x)dx, 



(183) 



Po{S- 


-tx) 


Po{S- 


-tx) 


Po{S- 


-tx) 



dS.t{x) 



dx 
dS^t{x) 



log po{S-tx) 



dS.t{x) 



dx 



dx 
dS.t{x) 
dx 



log [pois -tx))dx 
dS^tix) 



log 



dx 



dx; 



dx 



il primo termine nella parte destra della ( 183 ) non contribuisce a £ perché ponendo S-tX ~ y diventa 
uguale alla costante — J po{y) \ogpQ{y)dy, mentre il secondo può essere riscritto come, con lo stesso cambio 
di variabile, 



dS^tiSty) 



d{Sty) 



(184) _ypo(y)lo, 

dove l'uguaglianza nella ( 184[ ) discende dall'identità 

dS_t{Stx) dStx 



'^y = j poiy) log 



dSty 



dy 



(185) 



= 1. 



dStx dx 

Dalle equazioni del moto x — f{x) si ottiene che, se x = Sty e cr(x) — — div/(2;). 



(186) 



d 


dSty 




dSty 


dt 


dy 




dy 



<y{Sty), 



(187) 



e il tasso di creazione di entropia E diventa, grazie alla (186), 

è = - J po{y)cr{Sty)dy 

PQ{S-tz) (j{z)dz 
fit{dz)a{z) ^t^oo - / p{dz)a{z), 



ammettendo che fj,t{dx) converga debolmente a fi(dx). 

Quindi, se il sistema tende ad uno stato stazionario descritto da una misura di probabilità invariante 



p (che accettando l'ipotesi caotica coincide con la misura SRB nel futuro interpretiamo la (1871 come 
l'entropia creata dal sistema e ceduta al termostato per unità di tempo nello stato stazionario. Infine, c un 
risultato noto, cf. |20J, che per un sistema che verifica l'ipotesi caotica (a) > 0; in particolare, se la misura 
stazionaria non è assolutamente continua rispetto al volume si ha che (cr) > 0. 

4. Il teorema di fluttuazione: introduzione 

Accettando l'ipotesi caotica, il teorema di fluttuazione di Gallavotti - Cohen afferma che, chiamando 
a{x) il tasso di contrazione del volume dello spazio delle fasi e (cr)+ il suo valor medio con la misura SRB, 
introducendo la quantità £r come 



(188) Srix) = 

le fluttuazioni di Er verificano la proprietà 

1 



1 



(189) 



lim 



a{x(t))dt, 



■log ^+("-=^\ =1, 



se p+ (e-r = p) è la probabilità SRB dell'evento {x : £r{x) — p]. Nella prossima sezione vogliamo presen- 



tare una dimostrazione rigorosa del risultato (189) nel caso di un sistema di Anosov discreto {VL,S)] il caso 
nel continuo può essere ricondotto al discreto mediante la tecnica della sezione di Poincaré, cf. [6] . 



5. Dimostrazione del teorema di fluttuazione 



In questa sezione riportiamo una dimostrazione, tratta da 5j, della (189) nel caso di sistemi dinamici 
discreti {fl,S). Grazie alle tecniche esposte nel capitolo [5] vedremo che il problema sarà ridotto allo studio 
di un modello di Ising unidimensionale con potenziali a decadimento esponenziale. 

5.1. Notazioni. Prima di iniziare, ricordiamo qualche definizione. Se (fi, 5) è un sistema di Anosov, 
chiamiamo W^{x), W^{x) le porzioni connesse di varietà instabile e varietà stabile passanti per il punto x 



e contenute all'interno di una sfera di raggio 7 (c/. definizione 3.3), A(a::) il determinante jacobiano della 
mappa S : Q CI, e indichiamo con A„(x) la quantità 

(190) K{x)= H X{S^x); 

infine, con A„, Xs ci riferiremo ai determinanti jacobiani di S vista rispettivamente come mappa sulla 



varietà instabile e sulla varietà stabile e, analogamente alla (1901, 



(191) XuA^) = n >^u{S'x) 



3 = ~2 

¥-1 



(192) A,,„(a;) - [] Xs{S^x). 

Dunque, A„ può essere scritto come 



(193) A„(a;) = Xu,n{x)Xs,n{x)xn{x)., 

sin Q 2 X 
sina^S 2" a: 

5.2. La misura SRB approssimata. Introduciamo la misura iit.t come, per T > |, 



dove XnKx) = — — \ _„ / è il rapporto dei seni degli angoli tra W" e nei punti S^^ x e S^^x. 



(194) 



(195) 



Plt,t{F) = J fiT,T{dx)F{x) 

f(X(g^))A-UX(g^)) 



^<T^e{i,...,n}i,-^'^i 

dove (Jrp è una sequenza cr_T, ctt con ai e {1, n}, {1, n}p ' è l'insieme delle sequenze compa- 
tibilj^in [— r,T], X è il codice della dinamica simbolica e A^^t — Yl,j=-Tj2 = Y^'^j=~T/2 logA^oXo^^, 

se 61 è l'operatore di traslazione verso sinistra nello spazio delle sequenzerl le funzioni F{X{u_rp)), Xu{X{u_rp)) 
vanno pensate calcolate nella configurazione che si ottiene "prolungando" in modo compatibile le sequenze 



^c/. definizione 3.1 per evitare confusione con il "volume" T della misura /jy ^ indichiamo l'insieme delle sequenze di n 



simboli compatibi in un insieme A con {1, n}^ e non con {1, ra}y . 

^In questa sezione chiamiamo 9 e non r l'operatore di traslazione verso sinistra nello spazio delle sequenze per non creare 
confusione con il tempo t. 



finite all'infinito. In questo modo abbiamo raffinato la partizione di Markov di (£7,5), cf. 



definizione 3.2 ovvero abbiamo diviso Q in piccoli rettangoli identificati dalle sequenze a^-, e assegniamo 
ad F, Xu un valore per ogni rettangolo scegliendo un punto in ognuno di essi; ciò può essere fatto in modo 
completamente costruttivo, ovvero senza ricorrere all'assioma della scelta, dal momento che il numero M 
di rettangoli è finito, JV < n^^ . Inoltre, se T è molto grande questa scelta influenzerà "poco" il risulta- 
to del valor medio perché F, A„, X sono funzioni Hòlder continue. Sempre grazie alla Hòlder continuità, 
nel linguaggio della meccanica statistica dei sistemi di spin la scelta della continuazione delle sequenze Orp 
equivale a fissare le condizioni al bordo di un modello di Ising unidimensionale con potenziali a decadimento 
esponenziale; è noto che questi modelli non presentano transizioni di fase, cf. appendice D, ovvero non limite 
T — oo lo stato di Gibbs non è sensibile alla scelta delle condizioni al bordo. Quindi, possiamo definire il 
valor medio SRB nel futuro di qualunque osservabile Hòlder continua con il limite termodinamico 

(196) hm ^JiTAF) = lJi+{F); 

ci riferiremo alla misura /iT,r; per T, r fissati con T > ^, parlando di una misura SRB "approssimata" . 
5.3. Sistemi dissipativi reversibili. Assumiamo che il nostro sistema sia dissipativo, ovvero che 

(197) f ti+{dx) log X-\x) = {a)^>0, 
e che esistc0 un'isometria differenziabile / con P — ±1 tale che 

(198) /o 5 = 5-^0/; 



con la (198) stiamo assumendo che la dinamica sia reversibile. Considerando per semplicità il caso 
bidimensionale, chiamando Wu{y), Ws{y) due vettori normalizzati tangenti alla varietà instabile e alla varietà 
stabile nel punto y, possiamo definire A„ e Ag mediante le relazionQ (indicando con DS il differenziale di 5) 

(199) DS-\x)wuix) = \-\x)wu{S-'^x) 

(200) DSix)wsix) = X,{x)ws{Sx); 

poich^H/" (Ix) = IWix) e quindi 

''È importante verificare che sistemi di Anosov dissipativi reversibiU esistano; ciò può essere fatto ad esempio notando che 
se il sistema {fi, S) è dissipativo allora anche il sistema {fi' , S') con fi = fi X fi e S' (x,y) = {Sx, S~^y) lo sarà, e I o S = I o 
se I{x,y) = (y,x). 

cf. definizione 

*^Questa proprietà è una diretta conseguenza delle definizioni di W^{x), Wix), cf. definizione 3.3 infatti, se d è una 
distanza in fi, 

(201) = {y (^Q.\d[S-''x,S-''y) <y'd{x,y)} 

(202) W (x) = {?/ e n |d {ff'x, S"y) < X"d(x, y) } , 
e poiché / è un'isometria, i.e. non altera le distanze, ponedo y' = ly e assumendo che = 1, 

W"ilx) = {y en\d{S-"Ix,S^"y) < X"d{Ix,y)} 

= {ly' 6 n |d {lS"x, IS"y') < X"d{Ix, ly') } 

(203) = IW'ix). 



4.2 



(204) Wu{Ix)^DI{x)ws{x), 
grazie alla (198), 



(205) 



DS-\lx)wu{Ix) = DS-\lx)DI{x)ws{x) 
= D{S-^ o I){x)ws{x) 
= DI (Sx) DS{x)w,{x). 



Dunque, per la (205) la (199) calcolata in Ix diventa, moltiplicando la parte destra e la parte sinistra 
per [DI{Sx)]-\ 



(206) DS{x)ws{x) ^ \-\lx)ws{Sx), 

e questa equazione confrontata con la ([200) implica che 



(207) 



Nel caso d-dimensionale il ragionamento è analogo, con l'unica differenza che i vettori w„, Ws saranno 
sostituiti da due sottospazi con dimensioni o?„, ds, con du + dg = d. 

5.4. Dimostrazione del teorema. Introduciamo la media adimensionale su un tempo r del tasso di 
contrazione dello spazio delle fasi come 



t/2-1 



(208) 



£r{x) 



logA^^(a;); 



j = -T/2 



per quanto discusso nella sezione [sj la quantità definita nella (2081 può essere interpretata come la media 
adimensionale del tasso di produzione di entropia. Se /ìq è una misura di probabilità assolutamente continua 
rispetto al volume, per /iQ-quasi tutti i punti di D, si ha che 



(209) 



£t)+ = lim - V EriS^x) = / fi+{dx)er{y) = 1 



(210) 



lim er{x) = 



dal momento che nel limite r -> oo la quantità £r{x) non dipende piìi dal punto (c/. proposizione 5.1) 
e, per le pM, (pOTl), 



(211) 



lim £t{x) 



lim er{Ix)] 



nella (211) l'argomento del limite nella parte destra differisce da quello nella pa rte s in istra per una 
quantità che può essere maggiorata con se i?2 è una costante positiva (c/. formule (193), (208)). Siamo 
pronti per dimostrare il seguente risultato. 



Teorema (di fluttuazione) . Esiste p* > 1 tale che la misura SRB jj,^ verifica 
(212) p - ò < lim , , log -. — f < p + ò 

^i+{srix) e-[p-S,p + S]) -"^ 

Vp, \p\ < p* . 

Dimostrazione. Se X è il codice della dinamica simbolica (c/. capitolo[2j proposizione 3.1 ), la funzione 
er{x) può essere convertita in una funzione sulle configurazioni di spirj^ 

t/2-1 



(213) Sria)^- J2 ^(^'^)' 

fe=-r/2 

dove 9 è l'operatore di traslazione verso sinistra nello spazio delle sequenze compatibili e L{(j) ~ 
j^p^ log {X (a)) . Poiché L è Hòlder continua può essere espansa in potenziali a decadimento esponenziale, 

cf. proposizione |4. !( 



(214) Li<j) = Y,h-n.n][^[-n,n]); 



n>0 



dunque, chiamando X = (a;,a; + 1, ...,x + 2 m) u n generico intervallo reticolare connesso composto da 



2m + 1 simboli e X il suo centro, riscriviamo la (216) come 



(215) ^rilì^l _E Ixi^x)- 

Jf:Xe[-f,i-l] 

Le seguenti proposizioni sono equivalenti, cf. [6], |21j . 
Proposizione 5.1. Esiste p* > 1 tale che il limite 

(216) - log n+{er (a) e [p- S,p + S]) ^r^ao max -({s) 

T se[p-s,p+s] 

esiste e la funzione C(s) è analiti ca e strettamente convessa per p £ {—p*,p*); inoltre, la differenza tra la 
parte destra e la parte sinistra della (216) tende a zero maggiorata da ^ , per un'opportuna costante D > 0. 



Proposizione 5.2. La funzione A(/3) definita come 

(217) A(/3) = lim -log(e'^"<"> + ("--i)\ 

T— >oo T \ / + 



ovvero 



(218) A(/3)=max[/3(a) + (s-l)-C(s)], 

S 

è analitica in /3 ^ (—00,00) e asintoticamente lineare^ i.e. esiste > 1 tale che 

(219) lim r'A(/3) = (a>+(p*-l). 

p— ^00 



Con un abuso di notazione indicheremo con lo stesso simbolo le funzioni Sr, Et o X. 



Osservazione 5.1. (1) valore p* è il supremo dei valori possibili di EtÌ^x) nel limite t ^ oo, ovvero 

P* = sup^gf^linisup^e^(x), [5]. 



(2) Quindi, per \p\ < p* il limite nella (212) esiste; al di fuori di questo intervallo nulla può essere 
detto. 

Rimandiamo all'appendice [C] per una dimostrazione dell'analiticità di A(/3). Introduciamo la funzione 
èr come 



(220) 



X) 1 



jfc[-5,5-i] 



dove con Ix {(Lx) indichiamo ciò che rimane sottraendo ai potenziali Ix {{Lx) parte dovuta al rapporto 
dei seni, cf. formula (193|. La (2201 "approssimerà" Sr nel senso che esiste C > tale che 



(221) 



■èri < 



C 



notare che ora e o I = —e. Quindi, poiché 



(222) 

avremo 



(223) 



un discorso analogo può essere fatto per la misura approssimata ^J.T,T■ Per dimostrare il teorema sarà 
sufficiente provare che, se Ip s = [p — 3,p + S], 



(224) 



1 



log 



fJ-+{£TÌ^) ^ Ip^STvir)) / <p+S + ri'{T) 



{a)^ T /i+(er(o:) e I-p^5±n{T)) l >P-S-t]'{t) 



con ri'{T), rj{T) — ^r-s-oo 0. Per fare ciò useremo il fatto che le probabilità calcolate con la misura SRB 
possono essere stimate con la misura /iT,r, se T > ^; i n seguito considereremo il caso T — ^. 

Per quanto detto in precedenza, cf. formula ( 195 ) e discussione successiva, la probabilità SRB di una 



configurazione di spin = (t_z...ctz è data da, se $x sono i potenziali che si ottengono sviluppando 
nel senso della proposizione |4.1[ 



(225) 



^+(cr_r...crr) = lim 



dove l'asterisco sulla sommatoria indica che stiamo fissando la configurazione cr_x...(Tr, e sommando 

O 9 9' 



sulle altre compatibili. La (225) può essere riscritta come 



(226) 



^l+{a^r...ar) = 



con, se a è il tempo di mescolamento della matrice di compatibilità (cf. definizione 3.1 1, 



J2 



.X) 



(227) 



5Z *x(fTx) 

Js:^[-5,f],xn[-f-a,5+a]=0 

X^[-5,5],Xnhf-a,5+a]#0 



poiché ^„ è hòlder continua esistono b,b' > tali che (c/. proposizione 4.1) |$x {qLx)\ — j 
dunque 



(228) 



x^[-5,5],xn[-5-a,5+a]#0 



E 



< 2asup V |$x (o:x)l = log^- 



Grazie alla (228) possiamo maggiorare la (226) come, ricordando che ogni "spin" può assumere n valori, 



(229) 



(226) < 



L 2 ' 2 ' ^ 



dove con il doppio asterisco indichiamo che stiamo sommando sulle sequenze g^^_^ _t_^-^^j^t_^^ 

indipendenti da a^; quindi, i potenziali relativi ad intervalli reticolari che non intersecano [— | — a, | + a] 
fattorizzano tra numeratore e denominatore. Per la minorazione il discorso è analogo a quello che abbiamo 
appena fatto, quindi si ottiene che 



< 



1 *x(^!:x) 



(230) 



^+(cr_x...crr) 



"^xc[-f ,f 1 *x(iix) - 



> 



xc[- 



" 2 ' 2 ' 
I-X (H-X ) 



-xc[- 



-] *x(^!:x) 



Dalle (223), (230), dal momento che ir dipende da un numero finito di simboli, è chiaro che possiamo 
trovare una costante Bi tale che, per ogni T > ^, 



(231) < BiHT^r (èr e IpJ+aj < BiflT,r (sr G /p,5+2£ 

dove Et è ciò che si ottiene sottraendo a e,- la parte relativa al rapporto dei seni; in questo modo, 
£t = r(<T)+ ^u)r^J,T — ° I- Una Stima analoga può essere fatta per la minorazione, dunque 



(232) 



l M+fe, G/„5^2c) f < ]^ + 1 log t^T.Ae.^I,,s) 
T (cr) , ^ r ^ I > '°g^i J ^ Ino- MT.^(g^gJp.^) 



e otterremo la (212) mostrando che 



(233) 



1 fJ'T,TÌ£r e Ip,6) ( <P + S 



Il rapporto di probabilità nella ( 233 1 può essere scritto esplicitamente come 



(234) 



. ^uMXia^)) E,.,, ^ Xu!r{IX{a^)) A.,.(X(c7^)) ^ 



(235) 



Ee.gj,., A„^U^(gT)) f <max^^A-i(X(a^))A-i(^(flT)) 



E: 



\sAX{<Lt)) l >min^,A-UX(aj,))A-i(^(^ZT)) 



poiché il massimo e il minimo nella (235) sono calcolati rispettando la condizione 

1 



^^logA-UX(a^))A-^(X(a^)) e 7^,5 , 



la ( 233 ) è immediata; ciò conclude la dimostrazione del teorema. 



□ 



6. La reversibilità condizionata 

Consideriamo un'osservabile F con una parità ben definita sotto l'azione dell'inversione temporale /, 
ovvero 



(236) 



F{Ix) = rìFF{x) 



con rjp = ±1, e introduciamo una funzione continua ^{t). Vogliamo dimostrare che, nel contesto 
dell'ipotesi caotica e considerando sistemi dissipativi, 



(237) 



1 eC/^,e. e/p., 

— -nr — 7^ i 



^r— >-oo -L; 



dove U-y, fjpU^ sono dei "tubi" di larghezza 7 attorno a (p{k), rip(p(—k). In quel che segue nulla cambierà 
se al posto di F consideriamo un insieme di osservabili {Fi, Fn} con parità definita, ognuna con una 
traiettoria ipj] per comodità di notazione ne considereremo una sola. 

La dimostrazione della (237) è praticamente identica alla dimostrazione del teorema di fluttuazione 
riportata nella sezione [5] l'unica differenza è che la ( 234 ) è sostituita da 



(238) 



^e^eip,s-F(S''X(a^))eu^ A„ :^(X(o;y)) 

^e^eI-p,s,F{S''X{a^)eT,FU-, \i.\iX{aj,)) 
_ J2e^€lp,s,F{S''X(!J^))€U^ A„^^ (X(a-r)) 
^s-r<£lp,6,F{Si-X(g:^))eU., ^n)riIX{aj,)) 
_ J2e^eIp,s,F{S''Xia^y)£U^ KAX{aj,)) 
^e^eIp,s,F{S''Xia^y)<£U^ As,t(-^(CTt)) ' 

e da qui in poi tutto procede come già visto in precedenza. 



Possiamo riscrivere il teorema di fluttuazione (237) nella forma approssimata 



(239) 



ovvero 



fi+({FiS''x)}l^_^ eu-„ereip. 



H+{er G I-p,s) 



(240) 



^l+(^{F(s''x)}l_^^eU-, 









la formula (2401 ha un interessante significato fisico. 

Pensiamo di osservare la quantità £r su un intervallo temporale T ^ r ^ 1, in modo da averne ^ 
diverse misure; questi valori saranno molto spesso vicini ad 1 (dal momento che 1 è il valor medio di Et), e i 
valori lontani da 1 occorreranno in modo essenzialmente casuale. La frequenza con la quale si manifestano 
i valori —p è legata alla frequenza di p dal teorema di fluttuazione (212); ad esempio, la frequenza con la 
quale il valore p = — 1 si manifesterà rispetto alla frequenza del valore p = 1 è proporzionale a e~^^'^^+. 



Ciò 

valori p = 1 avremo un valore p — —1^ 
e in corrispondenza di questo dato, per il teorema (237), la probabilità che un'osservabile fisica segua il 



vuol dire che nella nostra raccolta di dati sperimentali per ~ e 



cammino invertito nel tempo sarà asintoticamente uguale alla probabilità che segua il cammino "corretto" 
con p — 1] per questa ragione ha senso dire che in un sistema che verifica l'ipotesi caotica la freccia temporale 
è intermittente, [22) . 

7. Altri risultati noti fuori dall'equilibrio 

Per concludere il capitolo, in questa sezione riassumiamo brevemente qualche altro risultato noto fuori 
dall'equilibrio. In realtà, per "non equilibrio" in questi casi si intenderà una cosa diversa da quella che 
abbiamo discusso finora; il primo risultato che presenteremo, noto come identità di Evans - Searles o teo- 
rema di fluttuazione transiente, descrive il comportamento di sistemi fuori dall'equilibrio in uno stato non 
stazionario, mentre gli altri due, la formula di Jarzynski e i teoremi di fluttuazione a molti indici di Bochov e 
Kuzovlev, trattano sistemi che vengono portati fuori dall'equilibrio attraverso trasformazioni hamiltoniane, 
i.e. senza dissipazione. 

Nonostante la loro generalità dunque, questi risultati (a parte il primo), non avranno la pretesa di 
descrivere il non equilibrio come lo abbiamo inteso finora, ovvero come stato stazionario di un sistema 
sottoposto a forze non conservative; anche se a prima vista potrebbero esserci delle similitudini con il teorema 
di fluttuazione di Gallavotti - Cohen, e bene tenere in mente che la situazione fisica che si sta prendendo in 
considerazione è completamente diversa. 

7.1. L'identità di Evans e Searles. In un articolo successivo ([23]) al lavoro "pioneristico" [2], 
Evans e Searles dimostrarono che, assumendo semplicemente che il sistema sia reversibile^ chiamando E±p 
l'insieme delle condizioni iniziali delle traiettorie lungo le quali il volume dello spazio delle fasi si contrae di 
gTPf^+T un tempo r, con ct_|_ pari al valore medio di a, allora 



(241) 



IJ'LiEp) 



IJ-lÌE-p) 

dove è la misura di Liouville nello spazio delle fasi. 



La dimostrazione della (241) può essere riassunta nel modo seguente: 



^ '^Affermando che un sistema è reversibile Evans e Searles intendono dire che se le equazioni del moto ammettono come 
soluzione la traiettoria 7r(to)— »r(ti) allora ammetteranno anche la soluzione invertita nel tempo 7r{ti)^r(to) > dove T{t) è lo 
stato del sistema al tempo t, specificato da tutti i valori di posizione e impulso delle particelle. 



__ p-dEp) 



(242) 



UlÌStEp) 

f^LjEp) 

fiL{Ep)e-^P''+ 



dove è stato usato il fatto che iil{IA) = I-Ll{A) (| det-D/| = 1). 

Le differenze principali con il teorema di fluttuazione di Gallavotti-Cohen sono le seguenti: 

• la misura è assolutamente continua rispetto al volume, al contrario della misura SRB la quale, nel 
contesto dell'ipotesi caotica, per un sistema dissipativo è concentrata in un insieme di misura di 
volume nulla; 



la ( 241 ) non ha termini di errore, ovvero è un'identità esatta a tutti i tempi r. 



Al contrario di quanto sostennero Evans e Searles in seguito, non è ovvio che la relazione (241 1 implichi 



il teorema di fluttuazione di Gallavotti-Cohen; secondo Evans e Searles, l'argomento per cui ciò dovrebbe 



essere possibile è il seguente. Scriviamo la (241) applicando l'operatore St alla la misura di Liouville /i^^, 

-Ti 



(243) 



ST{^iL){Ep) ^ ^ìl{S ^ Ep) ^ r(pa++0{^)) 

SrifiLÌiE^p) ^ìl{S-^E^p) ^ ' 



e effettuiamo il limite r — >■ cxd prima del limite T — >■ oo, in modo che il termine di errore all'esponente 
vada a zero; infine, nel limite T — >■ oo la misura S't(/^l) dovrebbe tendere alla misura stazionaria del sistema 
fuori dall'equilibrio, ovvero alla misura SRB accettando l'ipotesi caotica. 

La critica piti spontanea a questo argomento è che non è detto che sia possibile invertire i limiti in 
questo modo (nel teorema di Gallavotti-Cohen si effettua prima il limite che definisce la SRB e poi quello 
su r), e sicuramente non lo è senza fare prima delle assunzioni particolari. Da un punto di vista fisico si 
potrebbe obiettare che dovrebbe essere sufficiente prendere T "molto grande", in modo tale che S't(mì) sia 
"vicina" alla misura SRB; ma questa affermazione è decisamente poco ovvia, dal momento che la misura 
SRB è singolare rispetto alla misura di volume. 

Il fatto che il teorema di fiuttuazione e l'identità di Evans e Searles siano due cose ben distinte è 



ulteriormente evidenziato dalla presenza di esempi di semplici sistemi per i quali la ( 241 ) è valida, ma non 
il teorema di fluttuazione, |24 . 



7.2. La formula di Jarzynski. La formula di Jarzynsky, [25], permette di calcolare la differenza 
di energia libera AF tra due stati all'equilibrio collegati da una generica trasformazione hamiltoniana (un 
protocollo), tramite successive misure del lavoro W effettuato dal sistema sull'ambiente esterno al termine 
della trasformazione; il risultato di Jarzynski è il seguente: 



(244) 



-I3AF ^ (^g-PW^ 



con f3 = jr^, dove T è la temperatura dei termostati (che supponiamo costante) e il valore di aspettazione 
(...) è calcolato con la misura e~^^^(ipd(7, se coni/i, H2 indichiamo l'hamiltoniana del sistema rispettivamente 
allo stato iniziale e allo stato finale. 

Per una dimostrazione della (244| rimandiamo alla letteratura originale, |25j . Le ipotesi essenziali 
affinché la (244) sia vera sono due. 



• La trasformazione che collega i due stati deve essere hamiltoniana, ovvero se immaginiamo che il 
processo sia governato da un parametro A(t) ad ogni tempo l'hamiltoniana del sistema sarà data 

da Hx^t), con ì/a(o)=Hi e H^^tj^^^i^) = H2. 

• L'accoppiamento tra il sistema e il termostato deve essere debole: ciò vuol dire che se con G\(t) 
indichiamo l'hamiltoniana totale di sistema e termostato, i.e. 



(245) Gx(t){q.^) = Hx(t){q) + H{x) + h{q,x), 

dove H{x) rappresenta l'hamiltoniana del termostato e h{q^ x) l'interazione tra termostato e 
sistema, si deve avere che ^ 1-^(9)1 ^Q^x': in questo modo, l'interazione h{q,x) viene 

trascuratcEE • 

La quantità W che compare nella (244| merita una discussione; in realtà questa quantità è identificata 
con la variazione di energia totale di sistema e termostato, ovvero 



(246) M^ = GA(t,„)-GA(o), 

dal momento che sistema e termostato sono considerati isolati dal mondo esterno. Quindi in un espe- 



rimento, per avere le "risposte giuste" dalla (244), bisogna essere in grado di controllare completamente gli 



scambi di energia che avvengono con l'ambiente. 

7.3. Due precursori. Prima di passare a studiare le conseguenze del teorema di fluttuazione vogliamo 
concludere questa breve rassegna di risultati in meccanica statistica fuori dall'equilibrio discutendo un lavoro 
di G. N. Bochov e Y. E. Kuzovlev del 1981, cf. [26], |27) , che, ben 12 anni prima delle simulazioni numeriche 
di Evans, Cohen e Morriss, conteneva già alcuni degli elementi delle teorie sviluppate successivamente. 

Consideriamo un sistema classico che evolve con hamiltoniana Hq — H(){q,p), dove q — (gi, gs^r), 
p ~ (pi, ...,p3Ar) sono le coordinate e gli impulsi di N particelle; indichiamo un insieme di m osservabili 
macroscopiche con Q{q,p) — {Qq('?,p)}, dove a = 1, ...m. Supponiamo che il sistema sia all'equilibrio, e che 
dunque possa essere descritto con la misura canonica 

g-/3ffo(9,p) 



J^dpdqe-^"«^^-'P^' 



(247) po{q,p) 

JQ' 

dove con fi indichiamo lo spazio delle fasi. 

Al tempo t = accendiamo una perturbazione, che si manifesta con un cambiamento di hamiltoniana 



(248) Hiq,p,t)=Ho{q,p)~hiq,p,t), 



(249) hiq,p,t) ^ X^{t)Q^{q,p) = Xit)Q, 

dove gli indici ripetuti sono sommati; notare che, in questo modo, l'evoluzione del sistema rimane 
hamiltoniana, e quindi l'elemento di volume dello spazio delle fasi può contrarsi, ma solo mantendendo il suo 
volume costante. 

Ipotizziamo che il nostro insieme di osservabili Q = {Qa} abbia parità definita sotto inversione tempo- 
ralcQ ovvero 

^^Non è importante che \h\ <S dal momento che alla fine AF non dipenderà dal termostato. 
^"^L'inversione temporale è naturalmente identificata con il cambiamento di segno delle velocità delle particelle. 



(250) 

e inoltre che 



Qa{q,-p) ^ £aQ{q,p), = ±i 



(251) 



Ho{q, -p) = Ho{q,p); 



nel seguito, con (...) j^-^g-j indichiamo i valori medi calcolati con la misura canonica relativa all'hamiltoniana 



(2481, dove h{p,q,t) dato dalla (249), e la perturbazione X{0) è stata accesa al tempo 9. Definiamo /„, E 



come 



(252) 



Uq{t),pit)) 
E 



di 



Q{q{t),p{t)) 

t 

X{T)I{T)dT 



X{T)Q{T)dT 



(253) 



dT 



H^{q{T),p{T))dT- 



introducendo una funzione di prova u{t) tale che lim(_>.-|-oo u{t) = 0, e considerandcj^ delle perturbazioni 
X(t) caratterizzate da limf_j.±oo X{t) — 0, nel loro lavoro |26j gli autori dimostrano che 



(254) 



/ x{e) \ 



:X{-e) 



Notare che X(t)I{t) è l'incremento di energia interna per unità di tempo dovuto alle forze esterne 
(conservative), dunque la quantità E che appare nella ([254 1 rappresenta l'incremento totale di energia interna 
nel sistema. Introduciamo il funzionale Dg[u{T); X{t)] come 



(255) 



Dg[u{T);X{T)] = In (^e^-'^ uir)I(r)dr\^ 



x{e) 



è possibile dimostrare che, partendo dalla (254), 



(256) 



De[u{T) - PX{t)-X{t)]= D^e[-eu{~T)-eX[-T)]. 



Grazie alla (254), o equivalentemente alla (256), Bochov e Kuzovlev riescono ad ottenere delle relazioni 



simili alla generalizzazione del teorema di fluttuazione che abbiamo discusso nella sezione |6j i. e 

(257) P9[{Q„(T)};X(0)]e-^^ = Pe[KQ„(-T)} ; eX(-0)] 

(258) Pe[{/„(r)};X(0)]e-'^^ = Pe[K/„(-r)} ; eX(-0)], 

dove P0[A{t)] X[t)\ è la probabilità di A al tempo t quando il sistema è sottoposto ad una perturbazione 
X{t) accesa al tempo 9. Una conseguenza diretta delle (257), (258) è che 



(259) 



\^ /x{e) ' 



^■^Stranamente, Bochov e Kuzovlev nel loro articolo 1261 non discutono affatto il problema della convergenza degli integrali 
per i — >■ oo; effettuano i limiti senza fare assunzioni speciali su X(t) e su u{t). 



quindi, poiché per la convessità dell'esponenziale (e ^) > 
la condizione 



, dal momento che /3 > la (259) implica 



(260) 



(E) 



x{e) 



> 0. 



Infine, è possibile dimostrare che la (256) implica delle relazioni tra funzioni di correlazione a tutti gli 



ordini nella perturbazione. Indichiamo con {■■■) x(e) ^'^ aspettazioni troncate rispetto alla misura di probabilità 
canonica con hamiltioniana ( 248 ); se j^j^ è la derivata funzionale rispetto a u{t) allora V aspettazione troncata 
di ordine n è definita come 



(261) 



(/„,(ti).../„„(i„)) 



T _ (5" 

^(») " 6Ua,,{ti)...SUo.^{tr,) 



De[u{T);X{T)], 



i.e. Dg è il funzionale generatore delle aspettazioni troncate di I(t). Grazie alle (2561, (261 1 Bochov e 
Kuzovlev dimostrano che 



T 



51 1 (^"i(^i)----^",.(^«)^/9i('ri)--//3,(Tfe))L0. ]^X^^.(Tj)dri...drfe; 



(262) 



considerando il caso n = 1, ti = t, la, causalità implica 



(263) 



{n-t))ex(-e) - {H-t))o = {Qi-t))o = per i > 0, 



e sempre per la causalità nel caso n = l,ti = t gli integrali nella parte destra della (262) sono 



identicamente nulli se r,; > t. Con queste considerazioni la (262) si riduce a 

(-l)'=-i/3'= 



(264) 



(Ut)) 



E 



X(9) ^ ^, 
fc=l 



^a(i)^ft (t-i)...//3, (rfc))^(g^ J]^ (rfe)dTi...dTfe 

i=i 



ovvero è possibile determinare completamente la risposta non lineare del sistema. Al primo ordine in X 
abbiamo che 



(265) 



(Ut)) 



X(0) 



{Ut)I,(T))^X^{T)d7 



ovvero l'usuale teorema di fluttuazione-dissipazione. 

Il resto del lavoro è dedicato alla discussione di applicazioni delle relazioni (262) nel caso di modelli 
particolari, e al tentativo di ricavare equazioni simili per sistemi inizialmente in uno stato stazionario fuori 
dall'equilibrio; purtroppo però, la discussione diventa molto oscura. 

Ad ogni modo, i risultati da loro ottenuti ricordano molto nella forma quelli ottenuti a partire dalla metà 
degli anni '90; in particolare, vedremo che relazioni simili alle (262 ) (chiamate da Bochov e Kuzovlev teoremi 
di fluttuazione a molti indici) potranno essere dedotte dal teorema di fluttuazione di Gallavotti-Cohen, cf. 
capitolo |4j sezioni [T] [2j 



CAPITOLO 4 



Applicazioni del teorema di fluttuazione 



In questo capitolo mostreremo che, accettando l'ipotesi caotica, il teorema di fluttuazione di Gallavotti- 
Cohen implica delle relazioni tra funzioni di correlazione che includono la formula di Green-Kubo e le relazioni 
di reciprocità di Onsager, e in generale estendono la prima al di fuori della teoria della risposta lineare. 



1. Risposta non lineare 

Consideriamo un sistema di Anoso\|^ {Q, S), e definiamo Sr come la media adimensionale su un tempo 
T del tasso di contrazione del volume dello spazio delle fasi cr(x), ovvero 



(266) e.(x) = -i— J2 ^('^'^); 

il teorema di fluttuazione afferma cheS 



(267) — ^ log ^+/^" ^\ ^r^oo 1 se\p\<p\ 

dove /i+ è la misura SRB nel futuro, e /i+(er — p), (o')^ sono rispettivamente la probabilità dell'evento 
{x : £{x) = p} e il valor medio di (j{x) calcolati con la misura SRB. Nel seguito studieremo esclusivamente 
sistemi dissipativi, ovvero per i quali (cr)_|_ > 0; diremo che il sistema è dissipativo se G ^ 0, dove G sono 
i parametri che regolano le intensità delle forze non conservative. In generale, assumeremo che iy(x) sia 
analitica in x e in G, e che dunque abbia la forma, sottintendendo una somma sugli indici ripetuti, 



s 

(268) a{x) = J2g.J^°\x) + ^ G., ...G,,a,,...,, (x). 

i=l k>2 

1.1. Conseguenze del teorema di fluttuazione. Introduciamo la funzione X{f3) come 



(269) A(/3) = lim -log/e^"("-(")-i)<">+\ ; 

r— yoo T \ / + 

si può dimostrare che X{f3) è analitica in /3 e che è possibile portare le derivate in /3 dentro al valor medio 
SRB, cf. appendici \C\ [P] Dunque, definendo la quantità Ck come 



In generale, possiamo pensare che {fi, S) sia il sistema dinamico discreto che si ottiene osservando ad intervalli di tempo 
opportunamente prefissati un sistema dinamico che evolve in tempo continuo, cf. [6j, i.e. attraverso la tecnica della sezione di 
Poincaré. 

■^Per comodità di notazione, nel seguito con l^+{sr = itp) indicheremo sempre fi+{er £ Ip,s)i se Ip s c un intervallo 
[p — S,p + 5] con (5 > arbitrariamente piccolo. 



63 



(270) 



Qk 



= lim (r^-^aX [e,(a;) - l]'^ 



/3=0 



se (7+ = + 5 6 dove con (...)^ indichiamo l'aspettazione troncata, i.e. 

(271) (r''-i4 [£,(x) - 1]'^)^ = 9^A,(/3), A,(/3) = ^ (e^^(^-(^)-i)<^>+ , 

lo sviluppo in serie di Taylor della ([269| sarà dato da 



(272) 



/c>2 



fc! 



Il seguente lemma sarà ciò di cui avremo bisogno per studiare la risposta non lineare. 



Lemma 1.1. Il teorema di fluttuazione (261) implica che 



(273) X{P)^\{-l-p)-2{a)^p-{a)^. 

Dimostrazione. Se non effettuiamo il limite r ^ oo il teorema di fluttuazione diventa, ponendo 

TTrij)) = fl+{er ^ p) , 



(274) 



quindi, grazie alla (274) 



A(/3) 



lim -log / e^<'"> + (P~i)'^7r^(j3)dp 

lim - log / e^<''> + (f-i)^+^f<'"> + 7r^(-p)dp 

r^oo r J 

lim -log [ e^<''h'-P~'^^^~^''^Kr{p)dp-2(3{a)^^ (a)^ 



(275) 

e il lemma è dimostrato. 



r— >-cx3 T 

A(-l-/3)-2/3 {a),-{c7)., 



□ 



Ponendo /3 = la ( 273 ) diventa 
(276) 



M+ = A(-i) = E 



fc>2 



(-l)'Cfc 
fc! ' 



come vedremo, a partire dalla (2761 al secondo ordine, ovvero (cr)^|^'' = si può dedurre la formula 



di Green-Kubo, risultato ben noto in teoria della risposta lineare. Concludiamo la sezione dimostrando una 



generalizzazione della (276) 



Lemma 1.2. Il teorema di fluttuazione per il funzionale generatore delle aspettazioni troncate A(/3) \273 ) 
implica che 



(277) 



k>0 



(-l)''+"Cfc+„ 
fc! 



n>2. 



Dimostrazione. Scriviamo la relazione ( |273[ ) usando lo sviluppo di A(^) ( |272[ ), 



(278) 

poiché 



fc>2 



k>2 



k\ 



(279) 



la (278) diventa 



(280) 

e dunque, per n>2, 



k\ 

k>2 



E 



fc, E • )/3'^-^"-2(a)+/3-(a)+, 

k>2 ' j"=0 ^ 



(281) 
(282) 
(283) 



Qn 



fc>r 



fc- 



E 
E 

fe>0 



(-l)"Cfc 
(fc-n)! 

(-l)^+"Cfc+„ 
fc! 



□ 



1.2. Conseguenze del teorema di fluttuazione generalizzato. Adesso, consideriamo una generica 
osservabile O antisimmetrica sotto l'azione dell'inversione temporale /, ovvero tale che O o I — ~0, e 
chiamiamo rjr la sua media temporale adimensionale su un tempo r, ovvero 



(284) 



r(0)_ 



J2 OoS^. 



Se -K-T-ijì^q) = fJ'+i^r ~ PtVt — q) è la probabilità congiunta delle variabili p,q, con una dimostrazione 
praticamente identica a quella del teorema di fluttuazione ( 267 ) ( cf. capitolo [sj sezione |6]) si ottiene che 



(285) 



hm — log — ^ = 1. 



Analogamente al caso precedente, ponendo = (cr)_|_, 0+ = + i introduciamo il funzionale genera- 
tore delle aspettazioni troncate delle variabili rcr+ [e^ — 1] , tO+ [rjr — 1] come 



(286) 



A(/3i,/?2)= lim llog/e-(^i(-^(-)-i)<->++'3^(''^(-)-i)«^)+)\ 

r— >-oo T \ / - 



Dal momento che A(/3i,/32) è una funzione analiticcj^ per Pi, G (—00,00) può essere scritta come, 
sommando gli indici ripetuti. 



(287) 



00 „ 

A(/?i,/32) = ^^^/3ai.../3a,, «. = 1,2; 



k=2 



se ni, n2, ni + 712 = A;, sono rispettivamente il numero di ai = 1, aj = 2 in {ai, i coefficienti 

Cai...ak sono dati da 

(288) = dp^...dp\{Pi,P2) = lim (r'=-i<i [e, - 1]"^ O!^^ [r/, - 1]"^)";, 

dove 



(-289) = (^^r{0i(er(x)~l){a)^+P2(rìr(x)-l){O)^)\j 



In modo del tutto simile al caso precedente, dimostreremo che il teorema di fluttuazione (2851 implica 
una relazione per il funzionale generatore delle aspettazioni troncate A(/?i,/32)- 



Lemma 1.3. // teorema di fluttuazione [285) implica che 



(290) A(/3i,/32) = A(-l - /3i, -^2) ~ 2/3i - - 2/^2 . 

Dimostrazione. Come per la ( 273[ ), questa relazione si dimostra partendo da 



(291) 



7r,(p,g)=7r,(-p,-g)e^<->+feO(i); 



infatti, grazie alla (291 ) si ha che 



(292) 



A(/3i,/32) = lim -log /e"('^i(^'-^)<"^ + +^^(«-i)<'^>+)7r,(p,g)dpd(7 

r^oo T J 

^ lim - log [ e^i^^(P-^^<''h+P--(.i-^)('^) + )+-P(^) + 



XTTr{—p, —q)dpdq 



^00 T 

X7Tr{p^ q)dpdq 

lim ilog /e^((-'3i-i)(p-i)(->+-fc(9-i)(o>+) 



y,e<~^P^(-h-^M<^) + -i-) + )Trr{p,q)dpdq 
A(-l - (3i,-p2) - 2 /3i - 2 ;32 - (a)_ 



□ 



^Comc por il caso precedente, questa affermazione può essere dimostrata grazie alla tecnica illustrata nell'appendice [c| 



Ponendo /3i = — 5 , /32 = — 5 nella ( 290 ) si ottiene che 



(293) 



ovvero, sostituendo nella (293) A(/3i,^2) con il suo sviluppo (2871, 



(294) 



(295) 



E E 

k>2 a\...ak 

-E E 

k>2 Ql...Qfc 



A:! 



k\ 



2^ 



k>2 ^ ' Qi...afc 



A; Il 2 



_ ^ifc-(2i + l)2(2i + l) 

2'-' (2/ + l)!(fc-(2; + l))!^ 



se con ||a|| indichiamo la parie intera di a, e 



(296) 



lfe = 1,1,1, 



volte 



Dunque, mentre all'equilibrio il valor medio di un'osservabile antisimmetrica sotto / è identicamente 
zercj^ fuori dall'equilibrio diventa una somma di quantità O (G*^) per fc > 1; in particolare, (O)^ può 
essere espresso come somma di aspettazioni troncate con un numero dispari dj^a^ — 2. Questo fatto non è 
sorprendente, perché la proprietà O o I = —O implica che una qualunque rappresentazione di (C)^ debba 
necessariamente cambiare segno sostituendo O con O o J; ciò non potrebbe accadere se nello sviluppo ci 



fossero dei termini con un numero pari di = 2. Nella prossima sezione vedremo che la (295) all'ordine 



pili basso diventerà la formula di Green-Kubo; dunque, accettando l'ipotesi caotica la (295) può essere vista 
come una generalizzazione di quest'ultimo risultato. 



(2941, 



Notare che ponendo O = ct, la (295) si riduce alla (2731, poiché in questo caso Cai...ak = Ck e, nella 



(298) 



{}- (-1)^^=1 ^""•) - 2'= - (-1)^= (-1)"^.. Y ^ 



ai...Qfc Ql = l Q!fe = l 

Infine, per concludere la sezione dimostriamo il seguente risultato. 



Lemma 1.4. L'identità [290) implica che 



infatti, poiché |detD/| = 1 (detDZ è il determinante della matrice jacobiana di /), indicando con (0)0 la media di O 
calcolata con la misura di Liouville, 



(297) {O>o = - (O o /)o = - (O>o ^ = 0. 

^Ricordiamo che nella notazione Cc<i...Qj Qì = 1 indica la presenza di (t^_ \sr ~ 1], mentre Oi = 2 corrisponde a Oj^ \r\r — 1]. 



(299) 



Ci 



ia,....,2,2.2,, 



l volte n-l volte 



E 

k>0 



(-1) 



n-\-k 



kì 



-Ci. 



1,1 2,2.2., 



n>2. 



fc + i ì'olte Ti-Z volte 



Dimostrazione. La (290l diventa, sviluppando X(j3i,f32) aspettazioni troncate, 



nLi /^Qj p 



fc>2 



/c! 



(-i)^nti/?. 

fc! 



-Cr, 



2/3i (^t) 



(300) 

dove /3i = (1 + /3i) e 132 = P2; poiché 



k>2 

-2/32 (0>, - (a), , 



(301) 



E ^=E( ■ )(i+/3iF^2"'ci,i 



ai...Qfc i=l 



2,2... 

j volte k—j volte 



la ( 300 ) può essere scritta come 



- E^È( ■ )(i+AFr^ci,i..^, 



fe>2 



(302) 



k>2 i = l 

k>2 j = l i=l ^ 

-2/3i (a), -2/32 (0)+-('T>. 



j volte 



k—j volte 



dovcl^ = (l,l,l....),2,._^--(2,2,2... 
Quindi, per n > 2, 



/ volte n—l volte 



= Ci 



k>2 



E 



(-1)' ' f k-n + l 
fc! 



k — n + l 



l\{n-iy.C 



1,1..., 2,2,2... 

Z iioZte Ti-i volte 



1,1.1... . 2.2,2. 



fc-Ti + Z volte n-l volte 



(303) 



E 

k>0 



(-1) 



fc! 



-C i^i^i... ^ 2,2,2... • 



fc + i volte n-l volte 



Ponendo O = a, la (3031 si riduce alla (277). 



2. Risposta lineare 



AfEnché la teoria sia consistente è necessario verificare che i risultati appena ottenuti implichino all'ordine 
pili basso l'insieme di relazioni tra funzioni di correlazione che formano la teoria della risposta lineare, [1]; 
come vedremo tra poco, dal teorema di fluttuazione sarà possibile ricavare due risultati notevoli come la 
formula di Green - Kubo e le relazioni di reciprocità di Onsager. Il collegamento tra la teoria della risposta 
lineare e il teorema di fluttuazione è stato notato per la prima volta in |15j . ed approfondito nei lavori 
successivi [ 28], [29] , 

Per la (276), il teorema di fluttuazione all'ordine diventa 



(304) 



v(2) 



a 



(2) 



2 ' 



vogliamo mostrare che la relazione (304) contiene la formula di Green-Kubo, risultato meglio noto come 
teorema di fluttuazione-dissipazione. 

Poiché a{x) = J2i=i GìJI^\x) + 0{G'^) la (304| può essere riscritta come 



(305) 

dovcEl 



+ Ig=o 



G=0 



(306) 



a{x)fj,^{dx) 
dG,dGj(j{x)y.+ {dx) 



G=o 



G=o 



+ / dG,(7(x)dGM+{dx) + [i ^ j) 



+ ì '^{x)dG,dG,fJ-+{dx) 



G=0 



G=0 



Il primo addendo della (306) è zero perché a = —a o /, |dct9/| = 1 e la misura ^+ per G = diventa la 
misura di volume; poiché cr{x)\Q^Q = anche il terzo termine è nullo, dunque 



(307) dG^dG, (a) = i^dG, (jP)^ + Og^ (jf)) J 

Definiamo la corrente Ji{x) e il coejficiente di trasporto Lij come 



G=o 



(308) J^{x) = BgMx) 

(309) = dG,{-H^))+\G=r 

è chiaro che jf\x) = Mx)\^=^, e poiché Mx) = Y.k>Q con jf^ = O (G'=), jI''\Ix) = -jf\x), 



Questa espressione è chiaramente formale, dal momento che stiamo scambiando l'integrale con la derivata, e non è ben 
chiaro cosa sia la derivata della misura. Rimandiamo all'appendice |D] per una discussione rigorosa di questa questione. Notare 
che se il sistema imperturbato è un sistema di Anosov per G abbastanza piccoli il teorema di stabilità strutturale, cf. 1101 . 
afferma che anche il sistema perturbato lo sarà; in particolare, con una tecnica analoga a quella esposta nella sezione |6] del 
capitolo 2 si può dimostrare che i potenziali della misura SRB sono funzioni analitiche in G. 



da, U 



G=0 



(310) 



k>l 



G=0 



G=0 



Grazie alle (307), (310) la (3051 diventa (moltiplicando per ^ entrambi i membri) 



(311) 



Lij + Lji 



1 °° / 

=1 E (U"\s^-)jr\- 



k— — oo 



- .L 



7(0) 



G=0 



che, almeno per i — j, è la formula di Grecn-Kubo. La (295) al secondo ordine diventa 



(312) 



2 ' 



che può essere riscritta esplicitamente come 



G=0 



(313) 



fc — — OO 



r(0) 



G=0 



1 E ((^r(5'vr(-))„-(4°').(^. 



7(0) \ / 7(0) 



indicando con {■■■)(, la misura invariante all'equilibrio, i.e. la misura di Liouville. La (313) è il teorema 



di fluttuazione - dissipazione, e poiché la parte destra è simmetrica sotto lo scambio i o j con la (313) 
abbiamo ottenuto anche le relazioni di reciprocità di Onsager 



(314) 



Quindi, se il sistema ammette un'inversione temporale / il teorema di fluttuazione all'ordine più basso, 
il secondo, equivale alla formula di Green - Kubo. 

3. Una verifica del teorema di fiuttuazione 

Per concludere, applicheremo le relazioni ottenute nelle sezione [T] di questo capitolo per verificare il 
teorema di fluttuazione in un caso semplice. Consideriamo un sistema dinamico definito sul toro T^, ottenuto 
identificando i lati del quadrato [0, 27r] x [0, 27r], che evolve ad intervalli di tempo discreti con la mappa S 
data da 



(315) 



S 



1 1 
1 2 



questo semplice sistema dinamico è chiamato gatto di Arnold, ed è un tipico esempio di sistema di 
Anosov, cf. capitolo [2] 

Perturbiamo S in modo da avere che l'evoluzione sul toro sia descritta da 



(316) 



S,^p ^Si^ + elif) mod 2tt; 



il teorema di Anosov afferma che il sistema perturbato, per e sufficientemente piccoli, è ancora un sistema 
di Anosov dal momento che per e < Eq esiste un omeomorfismo H (un "cambiamento di coordinate") tale 
che 



(317) SeoH = HoS, 

cf. capitolo [2j In generale, questo omeomorfismo è analitico in e in un disco di raggio Eq (dove può 
essere espresso come H = l + J2k>i con /i''''' — 0{e^)), e solo Hòlder continuo in (p. 

Come abbiamo visto nel capitolo |3j affinché il teorema di fluttuazione sia vero è necessario che il sistema 
(ri, S) sia reversibile, ovvero che esista / tale che 

(318) /o5 = S'"^o/, 
e la mappa I deve essere differenziabile; indicando con 



(319) 



1 



r/2-1 



la media su un tempo t del tasso di contrazione del volume dello spazio delle fasi a, l'esistenza di 

se Au, As sono rispettivamente 



un'inversione temporale differenziabile implica che, cf. capitolo [3j sezio ne [5 



il logaritmo del coefffciente di espansione e di contrazione, cf. formula (2071 e definizione 4.2 



(320) 



A„ o / 



-A.. 



ovvero 



(321) 



Er O I 



o 



Nel caso del gatto di Arnold imperturbato l'inversione temporale esiste ed è data da 



(322) 



Io 



-1 
-1 1 



ed è (banalmente) differenziabile; ma in questo caso il teorema di fluttuazione diventa vuoto, dal momento 
che la contrazione dello spazio delle fasi a = — log |det£)S'| è identicamente uguale a zero. Nel caso della 
mappa perturbata S^, invece, questa quantità è una funzione non banale di e, x; inoltre, per e abbastanza 



piccoli la perturbazione non rompe la simmetria per inversione temporale (3181, dal momento che possiamo 
definire una mappa con la proprietà o — o ponendo 



(323) I, = HoIaoH-\ 

Ma come per la coniugazione H, Ie{'>P) sarà in generale analitica in e per |e| < e^ e solo Hòlder continua 
in V', dunque l'esistenza di Ig non basta per soddisfare le ipotesi del teorema di fiuttuazione. Inoltre, con un 
calcolo esplicito al primo ordine in e abbiamo che, cf. capitolo [2] sezione [6j 

(324) Ai^Hf) = X-+'d+f+{±) 

(325) 4i)(^) = -\z}d.f.{^). 



quindjj^M o Iq ~As . Perciò, ci aspettiamo che nel caso del gatto di Arnold perturbato il teorema 
di fluttuazione sia falso. 

3.1. Risultati numerici. In realtà, dai risultati di alcune simulazioni numeriche, figure [l] [2] non 
è affatto chiaro che in questo caso il teorema di fluttuazione sia violato, dal momento che l'andamento 
teorico (retta verde nelle figure [T] [5]) è consistente con la dispersione sperimentale; riportiamo i dettagli delle 
simulazioni. 

• Le dispersioni di punti delle figure [T] [2] sono state ottenute effettuando per N = 100 volte un numero 
T = 10^ di iterazioni della mappa perturbata rispettivamente con 



(326) l^if) ^ ( ) l^it) = ( ^"^('Z'^) 

scegliendo ogni volta delle condizioni iniziali casuali con un generatore di numeri random. 

• In ogni punto l'errore coincide con la deviazione standard della media di 100 misure. 

• Il tempo r in £^ è pari a 100. 

• La media {cr)rp equivale a ^ Y^J=o ° ^ = 10^- 




Figura 1. Andamento ^ log Num. occorrenze di{e.=p} gr T = 10», r = 100, e = 

T{cr}j, Mum. occorrenze di{ex=-p} 
0.05, perturbando con una singola armonica; ogni punto è il risultato della media di = 100 
simulazioni. 



'"All'argomento di A^^^ compare la mappa Iq e non perché nel sistema perturbato i coefficienti di espansione e di 
contrazione sono definiti nel punto H{'tp), cf. capitoloj^j sezione|6j in realtà, la funzione Auiip) corrisponde a j4u(J?(</))), se Au 
è il logaritmo del coefficiente di contrazione definito nella l|l36|. 




e e = 0.05, perturbando con due armoniche; ogni punto è il risultato della media di = 100 
simulazioni. 



Ispirati da questi risultati numerici, poiché i valori medi di funzioni analitiche in ìp e e sono analitici in 
e, cf. appendice |Dj proviamo a verificare ordine per ordine il teorema di fluttuazione mediante le relazioni 
tra funzioni di correlazione ricavate nella sezione [l] di questo capitolo. 

3.2. Primo caso: perturbazione composta da una singola armonica. Perturbiamo il gatto di 
Arnold con 



(327) 




con questa perturbazione, il tasso di contrazione dello spazio delle fasi nel punto (p — H(xIj) gT'^ è dato 
da, se con |I?S'e| indichiamo il modulo del determinante della matrice jacobiana di S^, 

(328) aif) = - log \DSe\ = - log(l + 2e cos(V'i)). 

Sia (T che la coniugazione iJ e il tasso di espansione sono funzioni di s analitiche in |e| < Eq, dunque 
esprimibili come 



(329) a{f) = 

k>l 

(330) Hif) = V. + ^/i('=)(V^) 

fe>i 

(331) A^il) - log(A+)+^AW(V^), 



fc>i 



con a ^^\h^^\j^u'^ = 0(e'=). In quel che segue dovremo calcolare esplicitamente le derivate di valori 
medi relativi alla misura SRB; la derivata coinvolgerà l'osservabile contenuto nel valor medio e la misura 
SRB, intesa come stato di Gibbs di un modello di Ising unidimensionale. Quindi diremo che, c/. appendice 

oo 

(332) - J2 , 

k— — oo 

e la sommatoria converge perché la misura SRB è mescolante, cf. [8] . 

3.2.1. Secondo ordine. Il teorema di fluttuazione all'ordine equivale alla relazioni 



.(2) 

(333) 



v(2) _ C'2 



V (2) (2) 

le quantità , C2 che vi compaiono possono essere scritte esplicitamente come 

(334) (a)f = (a«C/(i))^ + (a(2)+9„a(i)/»«)^ 

k— — oo 

00 



(335) = J2 , 

A;— — 00 

dove con (...)g indichiamo che i valori di aspettazione sono calcolati con la misura di Lebesgue fj,o, 
[/(•) — — J2T=~oo ° ^ indici ripetuti a = 1, 2 sono sommati. Le quantità cr^^), cr'-^-', AÌ^-*, ha'' sono 
date da (c/. capitolo [2] sezione |6]) 

(336) Ai'Ht) = A;ia+/+(^) = ^^^cos(^i) 

OO 

h^^\t) = -^a;(^+^)/+(5^V:) 

(337) = -E m^ ^, sin(5V-ei) 

®Como abbiamo visto nella sezione |2] la relazione (333l contiene la formula di Green-Kubo. 



(338) 

(339) 
(340) 



p=i 

= y = rSÌn(5-PV-ei) 

è(A- + l)^ - 

Cr(l)(V^) = -2cOs(V'l) 

a'^^^ìp) = 2cos^(V'i). 



È facile vedere che i valori di aspettazione nelle ( 334 1 , ( 335 ) si riducono essenzialmente a degli integrali 



del prodotto di due funzioni trigonometriche, calcolate in (pi e in {S'^f) ■ ci = S^^ipi + 812^2 (ej = (1 0)) 
con k arbitrario; l'unico contributo diverso da zero è quello dovuto a fc = 0, dal momento che 5*^'^^° = 1 e 
Sii è una funzione strettamente crescente di |fcQ Un semplice calcolo mostra che 



(341) 
(342) 



c: 



(2) 



2£^ 



dunque la relazione (<t 



(2) _ 



è verificata (e con essa il teorema di fluttuazione al secondo ordine in 



3.2.2. Terzo ordine. Al terzo ordine, le relazioni da verificare sono 



(343) 
(344) 



A3) 



C. 



(3) 



c. 



(3) 



- 0, 



dove 



(345) 



\(3) _ 



2 \ /o 



^Per A; > ciò è ovvio, dal momento che 5 ha elementi di matrice strettamente positivi. Per A; < l'afFermazione è vera 
perché = 5^ , e per k negativi è una funzione strettamente crescente di 



(346) 



(347) 



C. 



(3) 



— — OD 

oo 

CO 
OO 



k— — oo 
oo 



A;^ — oo 



(3) 



^ ((a o S'^ - (a) + )(a o - (a) + )(a - 

k.j — — oo 



(3) 
+ 



(1) 



Oltre alle quantità già discusse al secondo ordine, nelle (345l-(347) appaiono anche (c/. capitolo [2j 
sezione |6]) 



(348) 



(349) 



(350) 



oo \-(p+l) 

^„(A+ + l)è(A„ + l)è ' ^ ) ^) 

OO 

^Ar'a„/_(5-f^)/iW(5-''V^) 

oo \P^^ 

E , ,.r., , cos(5-pv: • is-p±) 



^(A_ + l)l(A„ + l)l 

^ A;ia„a+/+(v.)M.') + A;ia_/+(v)fci'^(V') 



A 



(A+ + 1)(A„ + 1) 



-E 



A 



-2(p+l) 



p=0 



(A+ + 1)(A_ + 1) 



cos 



(V'i)cos(5-(f+iV-ei) 



A^ 



2(A+ + 1) 



cos(7/;i 



(351) 



^(3)(^)^_^COs3(^l)- 



Analogamente al caso precedente, i valori di aspettazione presenti nelle (3451, (3461, (3471 consistono 
essenzialmente in integrali come 



(352) 



2ir 



difi dip2 



cos{S^ip ■ ei) cos{S^ip ■ ei) cos(9?i) 



dipi dip2 ^^^(^gk^^^ _,_ s'^2f2) cos{SliLpi + S{2f2) cos(<^i), 



27r 27r 

che danno un contributo diverso da zero solo per coppie di indici k, j tali che 



aiS\i + a2S{i + «3 = 
«15^2 + a2S{2 = 0, 



(353) 
(354) 

con tti = ±1; affinché l'integrale (352) dia un contributo diverso da zero almeno una coppia di queste 
equazioni deve essere soddisfatta. Come dimostreremo facilmente tra poco, ciò sarà impossibile. 

L'equazione (354) con ai = —«2 ha una soluzione ovvia per k = j, mentre se ai — «2 la soluzione sarà 
k — ~j, dal momento che 



(355) 



nk 



-s. 



-fc 



per 



le equazioni del tipo (354) non possono essere soddisfatte da nessun'altra coppia di indici, perché \S^p\ 
con a ^ j3 è una funzione strettamente crescente di Per j = k la (353) diventa ±1 = 0, mentre per 

j ~ —k l'equazione da studiare è 

(356) S-fi + S'fi'^ + «3 = 0, 

la quale non ha soluzione perché S^i > 1 e S^i = S22 > 1- 

Quindi C3 = 0, e la relazione (cr)]^ = — |— diventa = 0; anche al terzo ordine il teorema di fluttuazione 
è (banalmente) verificato. 

3.2.3. Quarto ordine. Proseguiamo nella nostra ricerca di una violazione del teorema di fluttuazione 
provando a verificare le relazioni 



(357) 
(358) 



C. 



(4) _ 



c: 



(4) 



\(4) 



a 



(4) 



(4) 



c 



(4) 



6 



24 



(4) / C, 



per dimostrare che a questo ordine il teorema di fluttuazione è violato è sufficiente trovare che C3 7^ 
Le quantità che appaiono nella (357) sono date da 

00 

^(4) _ 



(4) 



^ {{a o - {a)^)ia o - {a)^){a (a) + )) 

A;,j — — 00 



(4) 



A:,j — — 00 



(359) 



f2 (aWo5Vi)o5VWc/«)^ 



k.j — — oo 



^Pcr A; > l'afFermazionc è ovvia, vedere nota ijo]); per < ci si riconduce al caso A: > usando la proprietà i 355 \. 



(4) 



oo 

k,j.l— — oo 

3 ((a o 5^ - {a)^)ia o 5^' - {a)^))^ {{a o - - (a)^))^) 



(4) 



k,jJ — -~-oo 



oo 



k,j,l——oo 



(360) 



discutiamo le (|359|), ( |360[ ) separatamente. Nella ( |360[ ) compaiono sia degli integrali tipo 
(361) 



L 27r 2^ 

che dei prodotti di integrali come 



cos(S' 'ip ■ ei) cos{S^ip ■ ei) cos(S' ip ■ ei) cos(i^ • ei), 



(362) 



271 



cos(5 «^•ei)cos(S^V-ei) / cos(5 • ei) cos(v3 • ei), 



e dobbiamo capire per quali terne (fc, j, l) le ( 361[), ( 362 ) danno un contributo diverso da zero. Ripetendo 
l'analisi fatta per il secondo ordine, è chiaro che la (362 ) è diversa da zero solo per k = j, l = 0, mentre per 



la (361 ) il discorso è piià complesso; le equazioni da studiare sono 



(363) 
(364) 



ai^ii +02 5^1 +a3S'Ìi +«4 = 



dove Gì = ±1. Le equazioni (363 1, (|364|), scegliendo opportunamente gli a^, sono sicuramente vere per gli 



insiemi di indici {k,j,lGZ:k= j, l = 0}, {fc,j, ^ e Z : fc = l, j = 0}, {k,j,l G Z : j = l, k = 0}; vogliamo 
dimostrare che queste scelte sono le uniche possibili. Per fare ciò basta notare che, se S" è data dalla (315), 



(365) 



S- 



k+i 



= 2S'l2 + > 25i*^2 sek>0 



(366) = -Sr2' <-2S;l'+'^ ^231+' 

quindi, per due generici indici m, n G Z tali che \n\ > \m\ 



se k < 0, 



(367) 



\S['2\ > 215^21. 



Grazie alla (3671 è facile rendersi conto che l'unico modo per soddisfar e la (|364[ ) è con una combinazione 
tipo k = ±j, / — 0; tra queste, il caso k = — j , Z = è da escludere perché la ( 363 1 corrispondente diventerebbe 



(368) 



Sii + Sii + 03 + 04 0, 



che, poiché S^i,S^^ > 1 per ogni fc e Z, può essere risolta solo imponendo fc = 0, e questa scelta è 
già compresa nell'insieme {k,j,l E 1, : k = j,l ^ 0}. Dunque, l'insieme delle terne (fc, j, l) per le quali 



l'integrale (361) è diverso da zero è determinato dalle condizioni k = j, Z = 0, e da quelle analoghe ottenute 
permutando k, j, l. 

Passiamo adesso alla ( 359 ) ; dal momento che 



(369) 



(2) 



indicando a o con cr^ e sottintendendo le somme su fc, la (359) può essere riscritta come la somma 
delle seguenti tre quantità: 



(370) 
(371) 
(372) 



3(aW.«.(^) 



- 3 



- 3 



Per quanto appena discusso, la parte destra e la parte sinistra della (371 1 sono diverse da zero solo per 
k = j, mentre la (372), poiché U^^^ = —J^'iZ-oo^'i^^ o S'' = ~ (a++i) J2'il-oo coso5'', dà un contributo non 
nullo considerando gli indici {k,j, l E Z : k — j, l — 0}, {k,j, l E Z : k = l, j = 0}, {k,j, l E Z : j = l, k = 0} 
per la parte sinistra e {k,j,l E Z : k — j, l — 0} per la parte destra. 

Per quanto riguarda la (370), invece, scrivendo esplicitamente le quantità contenute nei valori di aspet- 
tazione si ottiene 



E E-TT^ ((cos(5'=.)cos(5^-)sin(-)sin(5M)o 

k,j=-oo i=0 ^ + ' 

-(cos(5'=.)cos(5^.))„(sin(.)sin(5M)o) 

oo oo 24A'~"'^ 

E (attt) ^i'^(^"'-))o 



E 



(373) 



(cos(5^-.)cos(^^.))o(sin(-)sin(5-'-))o); 



i prodotti di valori di aspettazione presenti nella ( 373 ) sono diversi da zero solo per combinazioni di 



indici che verificano le condizioni k = j, l — 0, ed è possibile verificare che queste condizioni sono le uniche 



che rendono non nulli gli altri termini presenti nella (373). Infatti, per quanto discusso in precedenza 



riguardo alla (360), le scelte di indici che rendono il valore di aspettazione del prodotto di quattro funzioni 
trigonometrichd^^l (calcolate rispettivamente in S^^p-ei, S^ip-ei, S^ip-ei, Lpi) diverso da zero devono avere due 
indici uguali e Fmdice rimanente uguale a zero; vogliamo dimostrare che nel caso dei valori di aspettazione 



del prodotto di due seni con due coseni presenti nella (373) le combinazioni utili sono quelle che rispettano 
k = jJ = G. 

Supponiamo che non sia così e consideriamo ad esempio una terna (fc, j, l) con k = 0, j = l, per la quale 
l'integrale da calcolare è 



(374) 



(cos(-)cos(5'-)sin(-)sin(5'-))o; 



la (374) può essere riscritta come 



Noi l'abbiamo dimostrato solo per quattro coseni, ad ogni modo il discorso può essere ripetuto, ottenendo gli stessi 
risultati, considerando in generale combinazioni di seni e coseni. 



(375) 



(sin(25'-)sin(2.))„. 



e la (3751 è diversa da zero solo per l — 0. Dunque, l'insieme di indici caratterizzato dalle condizioni 



k = 0, j = l contiene una sola combinazione utile, ovvero k = j = l — 0, la, quale però è già compresa tra 
quelle che verificano k ^ j, l = 0. 

In conclusione, grazie a tutte queste considerazioni è facile verificare che 



(376) 
(377) 



a 



(4) 



c, 



(4) 



6A,£^ 3 

(aT+T) + 2' 

3e\ 



quindi nel caso di una perturbazione composta da una singola armonica il teorema di fluttuazione è 
violato all'ordine e^. 



3.3. Secondo caso: perturbazione composta da due armoniche. Perturbiamo adesso il sistema 



con 



(378) 



sin(?/'i) + sin(2-ì/;i) 




e con questa perturbazione il tasso di contrazione dello spazio delle fasi è dato da 



(379) 



ct(V') = -log(l + 2£(cos(V'i) +2cos(2V'2))); 



è facile vedere che anche in questo caso il teorema di fluttuazione e vero all'ordine e^. All'ordine e'^, 
invece, le equazioni da studiare sono 



(380) 
(381) 



015^2+62^^2 = 0, 



dove bi può essere uguale a ±1, ±2. Per quanto detto in precedenza, le equazioni del tipo (381 1 possono 



avere soluzione solo per k = ±j; tra queste scelte, la (380) è soddisfatta solo con k = j = 0, nella forma 



(382) 

oppure 



S 11 ~t~ 1 — ^ 



(383) - S'fi - +2 = 0. 

Quindi, al contrario del caso caratterizzato da una perturbazione composta da una singola armonica. 



esiste una combinazione di indici che risolve le equazioni (380), (381) e dunque 



(384) C^/^= ('T^^^o^Vi^o^VW) =-12£V0; 

A;,j — — 00 

questo risultato implica che il teorema di fluttuazione sia violato già a partire dal terzo ordine in e. 



3.4. Conferme numeriche a posteriori. Possiamo provare a verificare numericamente i risultati ot- 
tenuti calcolando esplicitamente il funzionale di grandi deviazioni ({p) nei casi relativi alle due perturbazioni, 
e confrontando la differenza —Cip) + Ci~P) con le dispersioni sperimentali ottenute per vari valori di e. Il 
risultato è che la funzione C{p) è data da, al quarto ordine in e. 



(385) 



(p-i) 



C2 
4 



ip 



48 ^ 384 



ir 



Ci + O (e^) 



rimandiamo all'appendice|B]per i dettagli di questo calcolo. Notare che la presenza di potenze di {p— 1)'^ 
con k > 2 nella ( 385 1 implica che le fluttuazioni del tasso di produzione di entropia non sono governate da 
una legge gaussi ana. 

Grazie alla (3851, poiché {(t)_^ — O (e^) (c/. notai?! appendice [b| , 



(386) 
(387) 

(388) 



C(p) + C(-P) = 
A = 

B = 



p{a) 
1 



1 

24 



^A]+Bp^ + (e^) 

C2 C3 Ci 

'+ ^ T T ^ 48 

C4 

2 



+ O (e^) 



per verificare la violazione del teorema di fluttuazione studieremo l'andamento di A in funzione di e. In 
un andamento tipo quelli riportati nelle figure [T|2l la quantità A può essere misurata a partire dalla pendenza 
intorno ap — della retta che interpola la dispersione di punti sperimentali, e calcolando la deviazione da 1. 

In generale, gli andamenti sperimentali saranno affetti da un ulteriore errore dovuto al fatto che il tempo 
T di £7. è finito (nel nostro caso pari a 100); la correzione può essere stimata, |30| , e indicando con Cr(p) il 
funzionale di grandi deviazioni della probabilità Prob{er — p) si ottiene che 



(389) 



Cr{p) = ap) + o[^) 



Nel caso di una perturbazione composta da una si ngola a rmon ica ci aspettiamo che, a meno di errori 
dovuti alla finitezza di r, la quantità A definita nelle (386), (387) sia O (e^), mentre se la perturbazione 
contiene due armoniche A = O (e). Ciò è confermato 

dispersioni di punti sperimentali sono interpolate rispettivamente dalle curve 



agli andainenti riportati nelle figure [3j [4| in cui le 



(390) 
(391) 



f2{e) 



aie 



bi 



re 



a2£ + h2£ 



C2 
TE 



dove i parametri sono determinati minimizzando la somma dei quadrati dei residui. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

E 

Figura 3. Andamento di A in funzione di e nel caso di una perturbazione composta da 
una singola armonica; la dispersione di punti è interpolata dalla curva /i(e). 
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Figura 4. Andamento di A in funzione di e nel caso di una perturbazione composta da 
due armoniche; la dispersione di punti è interpolata dalla curva /2(£)- 



Conclusioni 



In questa tesi abbiamo discusso il problema della descrizione statistica dei sistemi fuori dall'equilibrio; 
in particolare, abbiamo studiato gli aspetti fisici e matematici di una possibile teoria microscopica degli stati 
stazionari di non equilibrio proposta da G. Gallavotti e da E. G. D. Cohen, [3], [4], |13j . |14| . Questa teoria 
si basa su rni^ assunzione, motivata da risultati sperimentali, che consiste nella seguente affermazione: 

Ipotesi caotica (IC): Allo scopo di misurare proprietà macroscopiche, un sistema di parti- 
celle in uno stato stazionario si comporta come un sistema di Anosov reversibile. 

La conseguenza principale di (IC) è l'esistenza di una misura di probabilità invariante, la misura SRB, 
che determina la statistica di dati iniziali scelti a caso con la misura di volume, c/. capitolo [2] sezione |5.2| 
in particolare, se il sistema dinamico è sottoposto a forze non conservative la misura invariante è singolare 
rispetto alla misura di Lebesgue, ovvero è concentrata in un insieme di misura di volume nulla. 

Una caratteristica molto importante della misura SRB e che può essere rappresentata come lo stato di 
Gibbs di un modello di Fisher, ovvero di un modello di spin unidimensionale con potenziali a molti corpi 
invarianti per traslazione che decadono esponenzialmente; questi sistemi possono essere studiati in dettaglio, 
e in particolare è possibile dimostrare che non presentano transizioni di fase, nel senso che se i potenziali sono 
analitici in un parametro (5 allora anche la pressione lo sarà, indipendentemente dalla taglia dei potenziali, 
c/. appendici [C] [Dj 

Il contatto tra la teoria matematica dei sistemi di Anosov e la Fisica è contenuto nel teorema di flut- 
tuazione di Gallavotti - Cohen; questo risultato è stato ispirato da una simulazione numerica in 2j, ed in 
seguito è stato verificato in molti altri esperimenti, cf. ad es. |15| . |30| . |31j . Consideriamo un sistema di 
Anosov (ri, iS) che evolve ad intervalli di tempo discretj^ e indichiamo con {cr)+ il valor medio SRB nel 
futuro del tasso di produzione di entropicT^ a e con Sr la sua media adimensionale su un tempo r, ovvero 



(392) ^r^^Y- 

J~ 2 

il teorema afferma che, ponendo Ip g = [p — 5,p + 8\ con 5 > Q arbitrariamente piccolo, e chiamando 
7r(eT- e Ip,s) la probabilità dell'evento {x : e-rix) £ Ip.s}, 

(393) ""/^^ ^ = e^P<^>++o(i) \p\ <p*>l. 

TT{er G I-p,S) 



In generale, possiamo pensare che (f2, S) sia il sistema dinamico che si ottiene osservando ad intervalli di tempo "oppor- 
tuni" un sistema dinamico che evolve in tempo continuo (un flusso, cf. definizione |1.2[ l; questa idea è alla base della tecnica della 
sezione di Poincaré. Il sistema dinamico risultante dalle "osservazioni" ha tutte le caratteristiche dei sistemi di Anosov, [6], e 
il teorema di fluttuazione può essere dimostrato in un modo completamente analogo a quello che abbiamo visto nel capitolo [s] 
sezione [s] 

^■^In realtà, il teorema è una relazione sulla fluttuazioni del tasso di contrazione dello spazio delle fasi; ad ogni modo, come 
abbiamo discusso nella sezione 3 del capitolo 3 questa quantità può essere interpretata come tasso di produzione di entropia. 
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Inoltre, se {Fi,...,Fn} è un insieme di osservabili con parità definita sotto l'azione dell'inversione 
temporale /, il risultato (393) può essere generalizzato nel seguente modo, cf. sezione [6] del capitolo |3] 



(394) 



■7r{Fj{S''x) ~ fj{k),er € Ip,s,j = l,...,n,ke [-f , f - 1]) 
TT{Fj{S''x) ~ Iipj{k),er € I-p,s,j = 1, ...,n,t e 



T T 

' 2 ' 2 



\P\ < P*, 



dove If{k) = rij(p{—k) se rjj è la parità di Fj sotto /, e con la notazione Fj{S''x) ^ fj{k) indichiamo che 
Fj{S'^x) è "vicina" a (p{k), nel senso che l'evoluto di Fj in un intervallo temporale [— |, ^ — 1] c contenuto 
in un "tubo" di larghezza 7 attorno a ipj . 

Nella parte originale della tesi abbiamo dimostrato che questi teoremi equivalgono ad un insieme di 
infinite relazioni tra funzioni di correlazione, che includono come casi particolari "vicino all'equilibrio" la 
formula di Green - Kubo e le relazioni di reciprocità di Onsager, cf. capitolo |4] In particolare, il valor medio 
di una qualunque osservabile O dispari sotto inversione temporale può essere scritto come 



(395) 



k>2 



fc-i I 
k II 2 



Ci 9 

j^fc-(2i + l):^(21 + l) 

(2; + l)!(fc- (2Z + 1))!' 



dove If. — 1, 1, 1, ||a|| è la parte intera di a i 



k volte 



(396) 
(397) 



(398) 



Ci 2 = \ìm d'^^d"p\\r{Pl,P2) 

—n—m T^OO ^ ^ 

A,(/3i,/?2) = l^e^(A(--i)(->++fe('*^-i)(o> + )^ 



r(0)_ 



in modo analogo è possibile dimostrare anche che 



(399) 



Ci. 



1,1,...., 2,2,2,., 



i volte Ti~i volte 



E 

fe>0 



k\ 



C 111.... 



2,2,2.... , 



n>2. 



k+l volte n-l volte 



Inoltre, cf. appendice [B] ammettendo che le forze non conservative siano caratterizzate da "piccole" 
intensità {G^}, abbiamo trovato un'espressione esplicita in serie di potenze di {Gì} del funzionale di grandi 
deviazioni C(p), definito come 



(400) 



C(p) = - lim T V(£^ g Ip^s); 



\k 



il nostro risultato è che Q{p) — X]a;>2 o^kip ~ 1) , dove i coefficienti Uk consistono in opportune combina- 
zioni di funzioni di correlazione a — a. Ad esempio, fino al quarto ordine in {Gì} si ottiene che, se Ck = Ci .0 



e Ci^^o è definita nella (396), 
(401) Cip) 



ip-l)' 



C2 
4 



(P - 1)- 
48 



C3 



384 



C4 + O (G^ 



'Come dimostriamo nell'appendice D nella 



i396i 



è possibile portare le derivate in /3i, /32 dentro al valor medio SRB. 



indicando con la notazione 0{G^) un termine che all'ordine dominante è proporzionale al prodotto di k 
variabili appartenenti a {Gi\. Dunque, poiché in generale 7^ se G 7^ 0, le fluttuazioni asintotiche di 
non seguono una distribuzione gaussiana 



Il teorema di fluttuazione ( 393 1 , nella forma 



(402) ~ap) + a~p)=p{<y)+ \p\<p\ 

implica che i coefficienti siano legati da delle precise relazioni; com'era prevedibile, abbiamo verificato 
fino al quarto ordine in {Gi\ che le relazioni che si ottengono sono esattamente le (395), ( |399[ ), sostituendo 
O con a. In particolare, queste relazioni non modificano il carattere non gaussiano di C,{p)i che al quarto 
ordine diventa 



(403) c(ri = ÒL^ft-<ii^ft(l+(^ + <i4Ì>!)+0(G^ 

Infine, abbiamo sfruttato i risultati ottenuti per verificare il teorema di fluttuazione ordine per ordine 
in un caso molto semplice, c/. sezione [3] del capitolo D in cui da alcune simulazioni numeriche il teorema 
sembrava essere verificato sebbene non ne fossero rispettate le ipotesi; la conclusione è che la violazione c'è, 
ma era troppo piccola per essere apprezzata nelle simulazioni. A posteriori abbiamo verificato numerica- 
mente questa violazione con delle nuove simulazioni, finalizzate allo studio della pendenza vicino a p = di 
p~^(o-)^^ [—C{p) + C(~P)]; le deviazioni osservate di questa quantità da 1, il valore atteso nel caso in cui il 
teorema di fluttuazione sia vero, sono decisamente consistenti con le nostre previsioni teoriche. 



'Ciò non è sorprendente, dal momento che la media Et è normalizzata con r e non con ^/t. 



APPENDICE A 



Derivazione euristica della misura SRB 



In questa appendice vogliamo presentare un semplice argomento euristico per dimostrare che la misura 
SRB può essere rappresentata come uno stato di Gibbs le cui probabilità condizionate verificano le equazioni 



(404) cxp f £ A„(rV") - A.(rV')ì , 

se ct' = {^jy^jyc), ct" = (cta^ac); seguiremo l'analisi svolta in [13j . 

Consideriamo un sistema di Anosov (il, 5), ed un suo punto fiss^O; un punto fisso O e 17 è caratterizzato 
dalla proprietà SO = O. Al tempo t = consideriamo una sfera B centrata in O contenente una massa 
unitaria distribuita con densità p; dopo un tempo T "molto grande" , l'insieme B sarà concentrato lungo 
la varietà instabile Wq passante per O. Un segmento 5 (consideriamo il caso bidimensionale per semplicità) 
giacente lungo la varietà instabile Wq n B centrato nel punto x corrisponderà ad un segmento S~'^5 G B 
centrato in S~^x, piìi "corto" di 6 di un fattore nr=o^ ^u^i^^'^^)^ dunque la massa che riveste S sarà data 

dea 



T-l 

(405) ^i{S)^phY[X-\S-''x)\S\, 

k=0 

dove h è i\ diametro della sfera S-y nella direzione stabile. Per grandi T l'evoluto S^B di B coprirà una 
porzione sempre più lunga di varietà instabile, quindi possiamo dire che la misura SRB, che è il limite per 
tempi infiniti della misura di volume, è concentrata sulla varietà instabile. Se consideriamo due segmenti S, 
5' appartenenti alla stessa porzione di varietà instabile Wy(a;) e centrati in due punti x, x' giacenti sulla 
stessa porzione di varietà stabile W^(x)^ in modo che i segmenti formino le basi di un "rettangolo" come 
illustrato in figura [l] il rapporto delle masse che li ricopre è dato da 



(406) 
(407) 



m(<5') 



Vur=lK\s-^x)v^l^^K\s^x)\s^ò\ 

v£r=l>^n\s-^x')v^^^^\-^\s^x)\sTb'\' 



la (407) diventa, nel limite T — ì' oo, 



Ad esempio, per il gatto di Arnold il punto i/j = (0,0) è banalmente un punto fisso. L'analisi euristica che stiamo per 
presentare rimane immutata se al posto di un punto fisso O consideriamo un •punto ■periodico O', ovvero un punto che è fisso 
per una certa iterata si può dimostrare, c/. [5], che per un sistema di Anosov (fi, S) i punti fissi sono densi in il. 

■^Se n è la dimensione della varietà stabile, dobbiamo moltiplicare per h" in modo da avere un risultato 

adimensionale. 



86 



6 



Figura 1. I punti a; e a;' si trovano sulla stessa varietà stabile, e i segmenti 5, S' sulla stessa 
varietà instabile. 



(408) 



lim 



dal momento che S'^S e S^S' tendono a coincidere. 



Il prodotto (408) converge perché i punti x e x' 



appartengono alla stessa varietà stabile e alla stessa varietà instabile, e le funzioni logA„(a;) sono Hòlder 
continue in x. Quindi, codificando i punti x e x' con le sequenze a' = (ct^^ac), a" = {g^'J^g^j^c) , con stessa 
continuazione in A'^ dal momento che i punti x, x' si trovano sulla stessa varietà stabile e instabile, la (408) 
diventa la (|404|). 



APPENDICE B 



Calcolo perturbai ivo di ({p) 



In questa appendice vogliamo discutere un calcolo perturbativo che ci permetterà di calcolare esplici- 
tamente, ad ogni ordine nei parametri perturbatori {Gì}, il funzionale di grandi deviazioni (^{p) definito 
come 

(409) Cip) = - lim - logM+(er = p)] 

r— foo T 



nella (409), fJ.+ {Sr — p) è la probabilità SRB dell'evento {x : Sr{x) = p} e, indicando con a il tasso di 
contrazione dello spazio delle fasj^a = — log \DS\, la quantità Sr è data da 



r/2-l 

(410) ^r = ^^ J2 



'+ ]=-t/2 

Alternativamente, C(p) può essere introdotta come la trasformata di Legendre di 

A(/3)== lim -log(e^''<'"> + (^--iA = lim -log / e^^^^^ + ^P-^K^iPÌàp 

r— >oo r \ / + r— >oo T J 

se TTrip) = M+(er = p), 



ovvero 



(411) Cip) = max [13 (a)^ ip - l) - A(/3)] . 

Dal momento che A(/3) è strettamente convesscj^e asintoticamente lineare, ovvero esist^p* > 1 tale che 

(412) lim r'A(/3) = (a)+(±p*-l), 
allora, per \p\ < p* , 

(413) /3(a)+(p-l)-A(/3) ^^^^o /? (p - p*) = -oo 

(414) /3(a)+(p-l)-A(/3) ^^^-oo /? (p + p*) = -c^; 

quindi se |p| < p* il massimo {3 — (3^ dell'argomento della trasformata di Legendre sarà raggiunto in un 
punto a derivata nulla, ovvero risolverà l'equazione 



^Comc al solito, ci limitiamo a considerare il caso di un sistema dinamico discreto (fi, 5). 

■^A{/3) può essere vista come la pressione di un modello di spin unidimensionale; è noto che, in assenza di transizioni di 
fase, l'energia libera è una funzione strettamente convessa nei suoi argomenti, |21| . |32| . 
^p* è definito come p* = sup^, limsup^ £t-(x). 
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(415) 



0, 



(416) 



che implica, sviluppando A(/3) in aspettazioni troncate, 



Co 



Ck 



k>3 



(fc-l)!C2 



Poiché in un sistema di Anosov "poco" perturbato i valori medi SRB di funzioni analitiche in {Gì} sono 
a loro volta analiticj^in {Gì}, e dal momento ch(|^la derivata in f3 di X{/3) è analitica in {Gì}, cf. appendici 
[cj[D]e [121, sotto la condizione|^ C2 ^ la quantità /3* è esprimibile come serie di potenze, e l'ordine n-esimo 
dello sviluppo può essere calcolato per iterazione sapendo ch^ 



(419) 



(420) 



pi'' = 



G. 



(2) 




Ck 



k>3 



[k - l)\G2 



in) 



n+2 , n—k+2 y „ s (n— m) 



fc=3 ^ ' m=" 

(a)+(p-l)n(") 



Co 



C2 

k=3 ^ '' m=0 {na:Vn,=mi=l 



ri+2 , n— fc+2 fc— 1 / \ (n— m) 



C2 



Questo fatto può essere dimostrato con una generalizzazione della tecnica che abbiamo usato per calcolare coniugazione 
e potenziali della misura SRB nel caso del gatto di Arnold, c/. |12 |. e grazie a quanto illustrato nelle appendici [d| [c| 
^CJ. nota|4] 

= implica, per la formula di Green - Kubo, che perturbando il sistema non si genera una risposta lineare; evitiamo 
di considerare questo tipo di "patologie" . A d og ni modo, questa condizione non è restrittiva perché se C2 = allora al 
denominatore avremmo C3, e la somma nella (416 \ partirebbe da 4 ecc. 

^Per dimostrare che {(7}_|_ = O (G?^) basta notare che se ^Q(dx) è la misura di volume normalizzata, introducendo il cambio 
di variabile x = Sy e ricordando che il tasso di produzione di entropia è identificato con il tasso di contrazione dello spazio delle 
fasi. 



(417) 



/ tJ.o{dx)= f Mo(dy) ^^^^ = / tio{d.y)e~''^y\ poiché 5-if1 = n; 



dunque, dal momento che per ipotesi (t(x)|q_q = 0, 



(418) 



= J2G,dG, (e- 



G^dG, (-'T(-)>o Ig=o = E G«9g. (-<t(-)>+ Ig=o. 



Grazie alle (419), (420) il punto /3* è calcolabile esplicitamente a tutti gli ordini e con esso anche C{p), 
essendo definita come 



(421) 



Cb)=/?* (a), (p-l)-A(/3,); 



in particolare, all'ordine n la (421) diventa 



71 — 2 ^1 

^^M )(p-i)_^i.(/3^^c.)'"^ 

m=0 fe=2 
ri-2 



^/3r' (c7)r"^(p-i) 



n -, n— /c 



(422) 



fc— 2 ^n— m ,n2 .. . ,n;^ ni — m « — 1 



Con la (422) possiamo, ad esempio, calcolare il funzionale di grandi deviazioni C(p) fino al quarto ordine 
in {Gì}, ottenendo che 



(423) 



Cip) 



C-2 



O" ^ 



Quindi, in generale la differenza — CCp) + Ci^p) sarà una serie di potenze tipo 



(424) 



k>0 



dove i coefficienti {afe} possono essere espressi in termini di Ci, (cr)^; il teorema di fluttuazione 

(425) -ap)+a-p)^p{^)+ 

implica che = per fc 7^ 1 e ai = (cr)_|_. Ciò sarà reso possibile dalle relazioni tra funzioni di 
correlazione che abbiamo ricavato con il teorema di fluttuazione per il funzionale generatore A(/3), c/. capitolo 
|4j sezione [ij Infatti, l'identità (4251 al quarto ordine equivale alle relazioni 



(426) 
(427) 



+ 



2 
2 



C3 C4 
6 24 



O 



grazie alle quali la ( 423 ) diventa: 



(428) 



ÌP-1) , ip-iy 



o(g5) 



Dunque, anche imponendo le (426), (427) la distribuzione di probabilità asintotica delle fluttuazioni di 
Er preserva il suo carattere non gaussiano. 



APPENDICE C 



Analiticità di A(/3) 

In questa appendice vogliamo discutere le proprietà di analiticità della funzione A(/3) definita come 

(429) A(/3) = lim T-Mog/e^'^<'"> + (^--iA = lini r^^ log /e^'^<'">+'^- \ -/3(a)+, 

T— >-00 \ / + T— i-OO \ / + 

dove Sr è data da 

5-1 

(430) er{X{a)) = ^ a {S^X{a)) ; 

-'~ 2 

per dimostrare che A(/3) è analitica in /3 G (— cx),oo) seguiremo l'analisi svolta in [33], [34], [8], |12| . 
Come abbiamo visto ad esempio nella sezione |5.2| del capitolo [Sj il valor medio SRB di un'osservabile O è 
definito dal limite termodinamico 

(43f ) = lim {O)^ , 

A— >oo 

se (...)a è la misura SRB "approssimata", i.e. uno stato di Gibbs contenuto in un volume finito A; 
dunque, 

(432) A(/3) + /3((t) , = lim lim r'Mog (e^''<''>+^-(^A . 

^ T— >oo A— >oo \ /a 



Il valor medio argomento del logaritmo nella (432| può essere scritto esplicitamente come, se {'l'jf} sono 
i potenziali della misura SRB, 

gr/3(a) + e,(X(^T^))-i:xcA*^(5;x) 

(433) /^r/3(o-) + e.M\ ^ 



y e-^A-cA*x(^;,) 

e poiché £r è sviluppabile in potenziali a decadimento esponenziale, c/. capitolo [3| sezione [5. 4[ possiamo 
generalizzare il problema studiando il rapportcj^ 



s:,cf-5,5ì['^*-(2-)+*-(^-)] .i*-(-x) 



e 



(434) 



^In realtà, nella (434i la somma sui potenziali '^x dovrebbe coinvolgere tutti gli insiemi X che intersecano [— J]; ma 



nel log B definito in i 436 i 



la differenza tra ques ta so mma e quella che appare nella (|434| è, se /3 è finito, stimata da una costante che può essere inglobata 



91 



dove con l'asterisco ricordiamo che le sequenze devono essere compatibili con (7\_t ti; inoltre, se a è il 



tempo di mescolamento della matrice di compatibilità (c/. definizione 3.1 ) 



(435) 



J2 ^xi^x) = '^xi'Lx) 

X!Z:[-i,i],Xn[-a-f,f+a]#0 



e il secondo addendo nella parte destra della ( 435 ) può essere stimato come, grazie al decadimento 
esponenziale dei potenziali. 



(436) 



< 2asup Y \^x [g_x)\ < logS > 0. 



X3^ 



Quindi, dal momento che le configurazioni di simboli con X [—a — |, ^ + a] =0 sono indipendenti 
da quelle relative agli intervalli contenuti in f , f ] , la (434| è maggiorata da, se {1, n} è lo spazio degli 
stati di un singolo spin. 



(437) 



(4341 < 



■^xrf-S sl['3*x(^x)+*x(^^ 



La minorazione può essere fatta in modo completamente analogo, e il risultato è che 



(438) 



(4341 { 



2 ' 2 J 



> 



che implica, poiché -B < oo, ponendo [/* = ExcA^^(— -f)' 



(439) 



lim lim T-Mog/e-''^*\ = lim A^Mog/e-*^^* 



in conclusione, il problema dell'analiticità di A(/3) può essere generalizzato studiando (il segno di /? è 
irrilevante) 



(440) 



AP(/3) = ^lim log {^e^*^* ^ ^ 



La ( 440 1 può essere riscritta come 
(441) 



AP(/3) = hm log - A"^ log 



dove 



(442) Z^^^e-^-cA[^*A-(^,)+*x]^ 

se sono i potenziali della misura SRB. Poiché sia '^x che $x decadono esponenzialmente, è la 
funzione di partizione di un modello di Fisher con potenziali a molti corpi 

(443) axi^x) ^ l3'fx{^x) + '^xi<Lxy^ 

per modello di Fisher intendiamo un sistema unidimensionale di spin interagenti attraverso potenziali che 
decadono esponenzialmente, e che sono diversi da zero solo per intervalli X connessj^ Dunque, il problema 
iniziale è ricondotto allo studio delle proprietà di analiticità della pressione di questo modello. 

In generale, questo tipo di problema può essere affrontato grazie alla tecnica della cluster expansion, 
[33j . [§], [12j . Questa tecnica però richiede che l'interazione iniziale sia debole; nel nostro caso invece /? è 
arbitrario, e inoltre non tutte le sequenze a sono ammissibili: ciò equivale a considerare dei potenziali con 
un cuore duro. Risolveremo il problema grazie ad una tecnica di gruppo di rinormalizzazione, [34j . 



1. Decimazione: prima parte 

Consideriamo un volume A C Z suddivisibile in intervalli {Bi}'^0' {-^i}'=o seguente modo: 

(444) A = So, Ho, Bi,Hi, Bi^i, Hi^i, Bi, 

dove 6 (Bi) T e N e i5 (Hj) = h G N. Quindi, stiamo assumendo che A contenga l{h + t + 2) + t + 1 
elementi; con la suddivisione (444), una generica sequenza di simboli S {1, ...,n}^ può essere rappresentata 



come 



(445) = {l3Q,rìo,Pi,r]i, ry,_.i, /?,) , 

dove ogni Pi, r]i contiene rispettivamente r + 1, h + 1 simboli. Poiché lo stato di Gibbs è unico, |21j . [8], 
siamo liberi di rappresentarlo chiamando spin /3j, rji. 

Come vedremo, integrando le variabili r/^ otterremo un sistema di spin /3j che vive su un reticolo rinor- 
malizzato, dove diremo che le distanze tra Bi e è pari ad uno, interagente attraverso un potenziale 
rinormalizzato a'"'^"; il nostro scopo è verificare che questo nuovo potenziale abbia tutte le caratteristiche 
necessarie per applicare la cluster expansion. Per risolvere il problema della compatibilità tra gli spin /3i 
sottintenderemo sempre che, se oq è il tempo di mescolamento della matrice di compatibilitèj^ T, h = nao 
con N 9 n > f . Una volta fissate le configurazioni Pi, le configurazioni rji dovranno rispettare dei vincoli di 
compatibilità con le prime; ad ogni modo, come sarà chiaro nel seguito, ciò sarà dal tutto irrilevante ai fini 
delle nostre stime. Possiamo scrivere l'energia totale del sistema come 



^Nel nostro cas o gli intervalli saranno formati da un numero dispari di siti, cf. proposizione 



'^cf. definizione 3.1 capitolo|2] 



4.1 



u- 



A _ 



(446) 

dove 
(447) 
(448) 



ax ÌQìx ) 

XCA 
l l-l 

Y.u'' (A) + E[^(/5-'?-^^+i)] 

i=0 i=0 



W(A,77.,ft+l) 



XcBi 

Y (^xÌqlx); 

X(ZBi\jHi\jBi+i:XnHi^(ts 



dunque, l'ultimo termine nella (446) coinvolge potenziali che intersecano più di un Hi, i.e. sono i 



potenziali a molti corpi nell'iZ-reticolo. Quindi, possiamo riscrivere la funzione di partizione come 



(449) = E E n«"^^''"' 

/3o,...,/3i i)o,--.,r)i-i \i=0 / 

dove con l'asterisco indichiamo la condizione X BiU HiU iJ^+i Vi. Definendo 



(450) 



,W(ft,,K,ft+l 



moltiplicando e dividendo per ni=o ^h. {Pi, A+i) la (449) diventa 
Za = 



i \ /i-i 



Po,. ..,13, \i=0 



Vi=0 



(451) 



nto^^(/3«.A+i) 



per ricondurci ad un modello di spin j3i interagenti attraverso un potenziale effettivo dobbiamo essere 



in grado di stimare il rapporto nella (451). I blocchi di spin Hi interagiscono attraverso i potenziali a molti 
corpi 

(452) aH,,...,H,+,{rì^,...,rì^+q) = ^ axi^x)'-. 

XnHk^9,k=i,...,i+q 

in particolare, ci interessa stimare la quantità 



(453) sup J2 sup|aff,,...,H, (77j/,,...,?7Hj|e5('= 

^ ■HBH:\H\>1 ^« 

dove H — {Hi, Hk) è un "intervallo" neh' iJ-reticolo e ij^ — {rjHi-, ■■■^'HHk)- Possiamo maggiorare la 



(4531 nel modo seguente: 



(453) 



E 



sup 

" H3H:\H\>1 ^ 



sup 



E 



< sup 

" ■H3H:\-H\>1 XnHi^(l>,i=l 



E 

xnffi5^0,i=i,...,|«| 



"X {q_x) 
sup |ax (o:_y)| e 



e 2 



f|fc-i| 



K<5(X)g-f5(Jf) 



(454) 



< 



fc>2 



1 



dove 



maxo- \ax (^x)!- La stima (454) è ciò di cui 



avremo bisogno per poter usare le tecniche della cluster expansion, che introdurremo nella prossima sezione. 

2. Cluster expansion 

La cluster expansion è una tecnica che permette di studiare somme tipo 



(455) 



dove {<Ì>y} sono potenziali che in generale possono non essere invarianti per traslazioni, e A è un sottoin- 
sieme finito l/] come abbiamo visto piìi volte, le funzioni di partizione dei modelli di spin rientrano nella 



classe di somme (455). La nostra tecnica si basa su un risultato fondamentale dovuto a Ruelle, |21) . che 
enuncieremo e dimosteremo nella prossima sezione. 

2.1. Polimerizzazione. Chiamiamo polimero qualunque insieme connesso finito ^ 7 C A, e indichia- 
mo con F una configurazione di polimeri F = {71, ...,7„}; non imporremo nessuna restrizione sui polimeri 
che formano F, i.e. permetteremo che i polimeri 7^ si sovrappongano e anche che coincidano. 

Definiamo attività di un polimero una qualsiasi funzione ^(7) definita sull'insieme dei polimeri 7 C A, e 
V attività di una configurazione di polimeri C, (F) come 

(456) C(r)-nc(^)^^"^ 

dove F(7) è la molteplicità del polimero 7 G F; quindi C, (F) è una funzione moltiplicativa, 

(457) C(ri + F2) = c(ri)C(r2), 

e imporremo che C((3) = 1- Diremo che una configurazione F è compatibile se non ci sono sovrapposizioni 
tra polimeri, i.e. se 7^ 07^ =0 per ogni 7^, 7j e F. La funzione di partizione delle configurazioni compatibili 
è data da 



(458) z{K)^ ^(r); 

rcA,7in7j=0 

definendo <i5(F), xv(r) tramite le proprietà 

(459) (p(F) 

(460) XV {V) 



j 1 se F è compatibile 

ì altrimenti 

r 1 se ogni 7 g F è contenuto in V 

ì altrimenti 



la (458) diventa 



(461) z(A) = ^<^(r)c(r)xA(r). 

r 

Introduciamo l(r) come 

e sia 5* una funzione sullo spazio dei polimeri del tipo 

(463) ^'(r) = l + *o(r) con*o(0)=O; 
possiamo definire il logaritmo e V esponenziale di ^(F) come 

(464) (Expvi/)(r) ^ i(r) + E^ E *(ri)...*(r„) 

n>i ■ ri+...+r„=r 

(465) (Logvi/)(r) = J2— — E *o(ri)...*o(r„). 

n>i ri+...+r„=r 

Quindi, chiamando V'^(r) il logaritmo di </3(r) si ha che, poiché XA(r), C(r) sono due funzioni moltipli- 
cativi 

(466) Za = E (E^P^"') (r)C(r)xA(r) = E [i(r)C(r)xA(r) 

r r 

+ Eà E '/''^(ri)C(ri)xA(ri) + ... + ^^(ri)C(ri)xA(ri) 

n>i ■ ri+...+r„=r 

(467) ^ e^rV^(r)C(r)xA(r)^ 
Siamo pronti per enunciare e discutere il seguente importante risultato. 

Proposizione 2.1. (Polimeri e cluster expansion) 

Sia A C Z'' un insieme finito e sia Cil) l'attività dei polimeri che verifica la condizione 

(468) |C(7)l<&^^o'^'e^'''"'^^'^^(7)', 
per opportune costanti òo,t'o,Ko > 0. Introducendo 

(469) fc(d,Ko) = 2e E e-""^^^, 
supponiamo che si abbia 

(470) [{l + Bo)e''"''° +k{d,Ko)] <1 



con Bq — 2bo. Allora la somma Z\ nella {458) può essere scritta come 



Notare che nella I467i il termine l(r)<^(r)xA(r) è diverso da zero solo per F = 0, e per questo valore è uguale a 1. 



(471) ZA = exp( ^^^(r)c(r)) , 

VrcA / 



con 



(472) sup V(^)\\C{l)\<B^voe 



^^'^r9C,(5(r)>r 



Dimostrazione. Chiamando F l'insieme dei polimeri distinti appartenenti a F = {71, ...,7„}, ponendo 
F! = n^ef -'^(t)' ^ chiamando G il grafo formato dai nodi {1, n} e dagli archi (i, j), j, j € {1, n}, posti 
in corrispondenza di una sovrapposizione tra due polimeri 7^,7^ G F, la funzione tp^ (T) che compare nella 



( 471 1 può essere scritta come 



1 

(473) ip^ ir) — — ^_2^(iiumero di archi in C) 

^' CCG 



La (473) prende il nome di formula dei coefficienti di Mayer, e l'asterisco indica che la somma coinvolge 
tutti i sottografi connessj^ di G che visitano tutti i nodi; se G non è connesso la somma non ha nessun 
addendo, dunque in questo caso 'p'^(r) = 0. 

Per convincerci della (|473[) dimostriamo che vale l'implicazione (473) 'p'^(r) = (Logiy9)(F). Per 



fare ciò, introduciamo una generica configurazione di polimeri F = {71, ...,7„}, dove ogni polimero ha una 
molteplicità F(7i), i = 1, ...,n; consideriamo distinte le F(7i) copie del polimero i-esimo aggiungendo ad ogni 
copia un nuovo indice 1, 2..., F(7i). Introducendo la funzione gij come 



(474) 5., 



-1 se 7i n 7j 7^ 
altrimenti 



poiché Lp{T) è definita dalla proprietà (460), possiamo porre 



(475) ^(F)^n(l+ff^^) 



i<3 



è chiaro che </?(F) = ^^yt^ '^^^ momento che se F! > 1 allora ip{T) = 0. 
Quindi, 

i<j k Ci...Ck non ordinati. Ci + ...+Cfc— F j 

= Eè E n^n^.V(Q) 

k ' Ci,... ,Ck ordinati, Ci + ...+Ck=r ™' j 

(476) = F!(Exp^^)(F)^ (Log^)(F) = V'^(F); 

nella prima somma l'ordine non dà luogo a nuovi addendi, i.e. le permutazioni di {Ci,...,Cfc} non 
vengono sommate, mentre nella seconda l'ordine conta (fattore ^) e con il doppio asterisco indichiamo che 
smettiamo di considerare distinti polimeri uguali (fattore y\^c ! ) • Definiamo la derivata di una fimzione 
definita sui polimeri come 



^Un grafo connesso è un grafo nel quale è possibile arrivare in qualunque nodo partendo da qualsiasi altro spostandosi 
sugli gli archi. 



(477) 



(r + g)! 

HI 



se r! è la molteplicità della configurazione T, i.e. il prodotto delle singole molteplicità che formano la 
configm'azione di polimeri; indicando con * l'operazione di convoluzione tra funzioni di polimeri, i.e. 



(478) 



(*i**2)(r)= J2 *i(ri)*2(r2), 



ri+r2=r 



è possibile dimostrare che valgono le proprietà 



(479) 
(480) 



D^(*i * *2) 



Ponendo 



(481) 



^(r) ^ C(r)^(r) 



che implicg[^Log(/3(r) = ^"^(T) = <!^{r)ip'^ (T), grazie alla definizione (477) introduciamo la funzione Ar 
come, indicando con vj/^i Vinversa di 'I', i.e. la funzione tale che '5^^ * "^i^) — l(r), 



(482) 



Hi+H-2=H 



è importante notare che Z?rv(-f^2) = '^{J^ + H2), dal momento che Lp{T + H2) è identicamente uguale 
a zero se {T + H2)\ > 1, quindi quando Tp è diversa da zero abbiamo che — 1. Quindi, poiché Tp 

si annulla se F! non è compatibile, i.e. se due polimeri che appartengono a F si sovrappongono, la (482) 
implica che 



(483) 



Ar(i?) = se F non è compatibile. 



Dalla (475) sappiamo che V9(7i, ...,7„) = n/<, (1 + 9{lt^li))i dunque 



^(7 + F + i/2) - ai)W + H2) n (l+ff(7,7')) 



(484) 



C(7)'?(F + i/2) (-1) 

se Ha 



N{S) 



se con l'asterisco indichiamo che stiamo sommando su tutti gli insiemi di polimeri S incompatibili con 
7, e dove la somma di polimeri 7 + H2 + F deve essere compatibile (altrimenti ^ = 0). Se 7 è compatibile 
con F dalla definizione (482) otteniamo che, ponendo H2 = S + H3 e grazie alla proprietà (484), 



Questo fatto è una conseguenza della definizione di logaritmo, cf. formula (465l, e della proprietà delle funzioni 
moltiplicative |457l. 



(485) 



SCH Hi + Hs^H-S 

SCH 



È importante notare che grazie all'espressione (4851 possiamo determinare una funzione di N{H) + 
N{T) + 7 = 771 + 1 polimeri attraverso quantità che dipendono da N(H) + A^(r) — m polimeri; grazie a 
questa proprietà saremo in grado di fare la seguente stima. Definiamo come 



(486) 

quindi, 



sup 

7i,...,7,,m>(}>l 



E |A{ 



H:N{H)=m-q 



(487) 

dove 



/„,+! = sup E |A^+r(H)|z(7 + r)-i 

r,l<|r|<m H:N(T)+N(H)=m 

< sup E E l^r+s(ff-^)lk(7 + r)rMC(7)l 
ra<|r|<m /i-.7v(r)+A'(_H")=m sqh 

< sup z(7) E MS) < z{j)i,n E E ^(^') 

r.i<|r|<m s:jv(s)<m-w(r) n>o y7'n7#0 / 

< z{j) exp E E - ^('^)-^™ ItI '^^P e ^("^'^ 



(488) 



^1 



(7) = feQj.lTle-«o5(7)elTl™P£«<i 2:^3; 2(7'), 



La quantità all'argomento dell'esponenziale nella (4881 può essere scritta esplicitamente come, c/. for- 



mula ( |468[ ) 
(489) 



sup E ^(7') - sup E 6o^^i^'e-''«^(^); 

7'9C 



è facile dimostrare che 



(490) 



sup E &o^i'''e-'=«'(^' < Bovo- 



Infatti, consideriamo in generale 



(491) 



sup E bov'^o^e-^"'^^'^- 



dal momento che il numero di alberi con n + l nodi di cui uno è una radice, i. e. è fisso, è maggiorato 
dsj^e^n!, chiamando ^ i, ^p) la lunghezza dell'albero 1} che ha per nodi i punti ^,^1, ■■■,£,p e ^ come 
radice, possiamo stimare la (491) come, indicando con l'asterisco la condizioncj^(5(^, ^1, ^g) > r: 



sup ^ bo^o 



l7lg-Ko5(7) 



™pE^ e - E E' 
™pE^E e - e ' 



-K0<5fl(?,?l:---,?<,) 



^-KO<5tì(?,?l,---,59) 



(493) 



q>o \ o^ceZ'' 



9>0 



dove nella ultima disugua glian za abbiamo usato che k(d, Ko)'^a < 1, cf. formula (470), e abbiamo posto 
So = 2Ò0. Inoltre, grazie alla d470| abbiamo che i^oe^""" < 1, 2boVoe^"''° < 1 dunque 



(494) 



^1(7) < 



1 



inoltre, sempre grazie alla (470), nel caso m = 1 la (486) diventa 



(495) 



h = sup |A^(0)| z(7)-i - sup ^ < boi^o < l, 
7 7 ^(7) ^ 



quindi, in conclusione la (487) implica che 



(496) 



per ogni m > 1. 



Grazie alla proprietà ( 480 ) abbiamo che 



(497) A,(r) = (^-1 * D,^) (r) = D.y^^iT) = ^^{j + r)^l±i^, 

quindi, 

''ciò può essere dimostrato nel modo seguente. Consideriamo n + 1 segmenti orientati indicizzati 1, 2..., n + 1, chiamiamo 
origine v di un segmento il suo punto iniziale e punti finali gli altri estremi. Tra questi, scegliamo un segmento c chiamiamo il 
suo punto finale radice; questo segmento è fisso, i.e. non parteciperà alla costruzione del grafo. In seguito, "attacchiamo" uno 
dopo l'altro i segmenti rimanenti ai punti finali dei segmenti già disposti. Chiamando "nodi" i punti iniziali di ogni segmento 
tranne quello associato alla radice, ogni nodo v avrà d„ segmenti entranti; quindi, il numero Af{n) di possibili grafi distinti può 
essere maggiorato con 



(492) 



n ri,* -. 

Min) = n\Yl J2j^<e'^nl 

«=idi=l 



dove con l'asterisco imponiamo la condizione J^ILi '^i — "1 '^^^ il fattore ni teniamo conto della permutazione dei vertici 
e i rapporti ^ servono a non contare le permutazioni equivalenti di rami che insistono sullo stesso nodo. 



'Ricordiamo che 5(7) = min^ 5^ (7)- 



(498) El^''(^ + r)| < E E \^,ir)\<zi^)J2i,-n = zh) 

r m>l r:A'(r)=m-l m>l 



e dunque, per le (493), (498), se Q e Z è un insieme che dista r dal punto p: 



(499) Yl \^^i^)\< E El^'^(^ + r)|< E ^(7) < So^oe-"»'-. 

r3p:rnQ#0 79p:7nQ#0 r 79p:7nQ5<i0 

□ 

Come vedremo nel prossimo corollario, le proprietà che abbiamo appena dimostrato con la proposizione 
|2.1| implicheranno l'esistenza del limite termodinamico per un gas di polimeri con attività invarianti per 
traslazioni. 

Corollario 2.1. (Analiticità e espansione in polimeri) 



Sotto le ipotesi della proposizione 2.L se ^(7) è invariante per traslazioni e analitica in un parametro /3 
allora il limite 

(500) P= lim Pa= lini A" ^ log Za 

A— fc50 A— foo 

esiste ed è analitico in j3. 

Dimostrazione. La pressione P\ dei polimeri con attività 7 — >■ ^(7) può essere scritta come 

(501) Pa = log Za = A^^ ^ ¥'^(r)C(r), 

rcA 



e sotto le ipotesi della proposizione 2.1 sappiamo che 

(502) sup |¥'^(r)||C(r)| <Boi^oe-''°'-. 

^^'^r95,(5(r)>r 

Se e è invariante per traslazioni, 

(503) hm A- Y ^"(r)C(r) = E "^^^W^^ 

rcA rgo ^ ' 

dove niV) è il numero di punti contenuti in F; inoltre, chiamando A la costante indipendente da (3 
necessaria per stimare i termini al bordcjjè facile vedere che 



^La costante A può essere stimata come 

(504) A<sup5:^!:(EM£^<Bo.o; 

nelle applicazioni alla meccanica statistica, in generale la costante vq sarà proporzionale a ell*lli, e # dipende da /3. Ma 
questo non rappresenta un problema, perché nei casi in cui sapremo ricondurci alla proposizione [2?T] la quantità ||3>||^ sarà 
limitata, in particolare piccola. 



E 



^^(r)C(r) 



< 



< 



E 

r9C:|r|>A 



n(r) 

^^(r)C(r) 



n(r) 



E 

r90,|r|<A 
A 
A 



'^^(r)C(r) 



n(r) 



A 

A 



E E 

p=or9e,(5(r)>A+p 



A 



(505) 



p=0 



A 



uniformemente in /3, c/. nota [9] Infine, se ^(7) è invariante per traslazioni 



|PaI < Ek^(r)|IC(r)| + ^<^ ^ k^(r)| |C(r)| + ^ 



(506) 



< BoUq ^ ( 

r>0 



A 

r>0 r90,<5(r)>r 

A = A/ < 00; 



allora, poiché P\ è analitica in j3 per il teorema della convergenza di Vitali, cf. 
P = limA^oo ^'a lo sarà. 



anche il limite 



□ 



2.2. Cluster expansion. In questa sezione discutiamo come "polimerizzare" lo studio delle proprietà 
di analiticità di un sistema di spin in una dimensione, in modo da rendere applicabili i risultati della sezione 
precedente; l'analisi che stiamo per fare può essere resa molto piià generale, [5], ma al prezzo di imporre 
condizioni piìi restrittive sulla forma dei potenziali d'interazione. 

Proposizione 2.2. (Cluster expansion per un modello di Fisher) 

Consideriamo un sistema di ^ spin su Z senza interazioni di tipo "cuore duro", ovvero senza vincoli 
di compatibilità. Sia $ un potenziale non necessariamente invariante per traslazioni; ponendo H'ì'xll = 
maxo-^ l'I'x (f^A:)! definiamo, per k > 0, 



(507) 



|cì>||^^sup \\^x\\e-'^''\ 

^ X3,\X\>1 



dove 6{X) è la lunghezza dell'intervallo X, i.e. 6(X) — \X\ — 1. 

(i) Indipendentemente dal fatto che $ sia invariante per traslazioni, la funzione di partizione = 
e~ ^^cA *x(o:^) py^fj essere scritta formalmente come 



(508) 



zt 



CeA 



exp[^^^(X)C(X)], 



dove ziX) = X^CT ^ *«('^)^ X = {71,.. .,771} è una configurazione di polimeri, Lp^{X) è la funzione 
combinatoria definita nella (-173), e C(^) = n7GX C(7) pe'^ certo ^(7) che verifica 



(509) 



IC(7)I < (■ 



e2(|l*IL+i)g-i 



l7l 



-1-5(7) 



La l[50^ diventa un'identità se la somma all'esponente è assolutamente convergente, 
(ii) Se i potenziali $ verificano la disuguaglianza 



(510) 



e^ll^lLe-J'^K-^ < e per < e < 1, 



allora la [509) è sostituita da 



(511) 



IC(7)I < e-^°l'^le-5'^*('^), 



con e " = ec per un'opportuna costante c; inoltre, la (510) implica che 



(512) 



sup 



E 



|^^(X)||C(X)||<i?ee"t' 



X3Ì,5[X)>r 



(Hi) Infine, se la |5jO| ) è verificata e se le funzioni (^(7) sono analitiche in $, nel caso invariante per 
traslazioni la misura di Gibbs /i$ è analitica in 



Osservazione 2.1. Come vedremo, grazie alla tecnica di decimazione che abbiamo introdotto nella sezione 
[7] di questa appendice, in un modello di Fisher sarà sempre possibile soddisfare la la ( 5l0j) . 



Dimostrazione. Innanzitutto, dobbiamo precisare cosa intendiamo per polimero in questo caso. Chia- 
miamo polimero qualunque insieme 7 C A C Z'' tale che per ogni coppia di punti ^, 77 € 7 esiste una 
"catena" Yi, ...,1^^ di sottoinsiemi di 7 con ^ G Yi, 77 € Ym, Yi n Y^+i 7^ 0, e per ogni Yi esiste almeno una 
configurazione CTy. con '^Yì {Q-y^ 7^ 



a 



Detto ciò, la funzione di partizione può essere scritta come 



<Ta6{1,...,«}A 



(513) 



n [( 



-*y(<Il,) _ 



1+1 



ycA,|y|>i 



e possiamo sviluppare il fattore binomiale HrcA l^^i "^^Ì-y) — 1^ + l 
relativi a polimeri compatibili, ovvero tali che 7^ H 7^ = 0; quindi la (513) diventa 



raggruppando i termini 



Zi = ^(A) E n 

7i,...,7p 1=1 



1 



(514) Ill^ 

dove 2(A) = n^eA Scr l'asterisco e il doppio asterisco indicano rispettivamente che i polimeri 

non si intersecano e che gli insiemi |Yj^*''| devono contenere una catena che connette due qualsiasi punti di 
7i. Grazie alla (514) è naturale introdurre le attività dei polimeri 7 e delle configurazioni di polimeri X come 



10 



Nel senso che la pressione limyv^txj ^ ^ '°S-^A ^ analitica in <& 
^Questa definizione in una dimensione è banale; lo è di meno per d > 1. 



C(7) 



(515) 



(516) 



C(^) ^ 



E E n 

U,>',=7 

nc(7)'''^\ c(0)^i, 

75^ 



e l'ultima cosa che ci rimane da fare è verificare che l'attività (^(7) così introdotta verifichi la stima (468) 
necessaria per applicare i risultati della proposizione |2.1| e del corollario |2.1[ Ciò può essere fatto nel modo 
seguente: 



IC(7)I < 



< 



|7| 9 . , . 

E E nU"-^^^-^^ -i 

l7l 9 / / 

^ ^ n (e-^'^^^^-^^-i 



hi q I 



4 g=ii=iY;35 

l7l 

9=1 

< (e2||*IL„+iy^lg-co5(7) 



(da |e^ - 1| < |a;|el^l) 



(c/. (507)) 



(e2||*||,„+iy'g-co5W (|x|<el-l), 
e, poiché per d = \ 5(7) = I7I — 1, quindi per polimeri con 5(7) > 1 (5(7) > \\^\^ 



(517) 



IC(7)I < (■ 



e2(ll*IL+i)e-fV''e-t^(^). 



Ponendo 6q = 1, ì^q = e^"*"+^ 2 ^ il risultato finale segue per la proposizione 2.1 e per il corollario 
sotto la condizione 



2.1 



(518) 

con Bq 



3e^°'"'+fc(l,|)]t.o<l 
= 2 e (c/. proposizione 2.1 formula (469l) 



(519) 



fc(l,Ko) = 2e J2 e^^ol^l < 46"'^° (1 + Kg 1) < 



?6Z,|C|>1 



e 2 



si può dimostrare che la condizione (5181 c implicata dalla (510). Infatti, definiamo e come 



(520) 



e grazie alla (510) sappiamo che esiste una costante c > tale che 



(521) 



e-^o = ce; 



inoltre, è facile vedere che i^q < e ", dunque la (510) implica che 



(522) 



i/q < ce. 



Infine, poiché per le (520), (521) 



(523) 



6 2 e 4 
< 



è possibile trovare un e > tale che 



(524) 



3e 



2uo 



16- 



.e 2 



ce e 



1 1 

i/Q < C2£2 



-2||<I>||,-1 



< ce, 



3e 



1 1 

2c2 e 2 



16ce 



< 1. 



□ 



3. Decimazione: seconda parte 

In questa sezione vogliamo mostrare che il nostro modello decimato, cf. sezione [Tj soddisfa le ipotesi 
della proposizione |2.2[ Avevamo riscritto la funzione di partizione come 



Za 



\i=0 



(525) 



/3o,...,A \i=0 



dove l'asterisco ricorda la condizione X (/i BilJ Hi\J S^+iVi, i.e. coinvolge intervalli che intersecano 



almeno due elementi dell'iJ-reticolo, e abbiamo dimostrato che (c/. equazione (455)) 



(526) 



E 



sup 

^ H3H:\H\>1 ^ 



sup 



< m IlalL. e 2' 



Adesso, se il potenziale iniziale aveva un raggio d'interazione k~^, il potenziale rinormalizzato tra i 
blocchi Hi che si ottiene fissando le variabili f3i avrà un raggio d'interazione ordine (tk) ^; dunque, possiamo 
riformulare la condizione (5101 come, sostituendo ||a|| con m ||q;||^ e~ 2^^^ 



(527) 



(tu) < e. 



Senza perdita di generalità, supponiamo che t — cosi ■ h: ciò implica che per h opportunamente grandi 



la (527) sarà sicuramente verificata. Il rapporto nella (525) ha la familiare forma di una media sugli spin ì]ì, 
-w{fìi,rìi,fii+i) . q^jQfjj^ possiamo ripetere l'argomento esposto nella dimostrazione della proposizione 



con pesi e 



2.2 e dunque, se F è una configurazione di polimeri appartenenti all'iJ-reticolo, il rapporto nella (525) può 
essere riscritto come 



(528) 



exp^^^(r)c(r). 



dove la somma all'esponente è assolutamente convergente per la proprietà (527), cf. proposizione 2.2 
grazie alla quale le attività C(7) verificano la stimEpn 



(529) 



IC(7)I < e"5'«'^lTle~J''^''('') < e~5''^l''l"5«^'5(T) 



se le distanze S{^) sono misurate suU' iJ-reticolo (quindi essenzialmente in unità t + h + 2). Definiamo 
i nuovi potenziali efficaci 



(530) 



*B„...,B,+,(/3„...,ft+p)= '/''^(r)C(r), 

rnBi#0,i=o,...,p 



che, lontano dai bordi del _B-reticolo, riflettono l'invarianza per traslazioni di ax- Grazie alla (527), 
dalla proposizione |2.2| sappiamo che 



(531) 



sup k^(r)|IC(r)| <i?o^^oe~''°^ 



r35,<5(r)>r 



con kq = cf. stima (529); dunque 



mfiyi\^B,,...,B,+^if3j,:;l3j+p)\ < ^(g + l)sup ^ 



k^(7)|IC(r)| 



q=0 



r9?,<5(r)>/i(p-i)+9 



< £(g + l)i?oi'oe-^("(P-^' 



+9) 



9=0 



(532) 



9=0 



dal momento che |C(7)| < ^o'''''^ «o^It) con — e "e , = cf. proposizione 2.1 Quindi, il 
potenziale efficace ^' può essere reso arbitrariamente piccolo. Chiamando X un intervallo nel B-reticolo, 



grazie alla ( 532 ) abbiamo che 



(533) 



1*11^ = sup J2 Il*x||e3i^(^' 

^ X3B.\X\>1 
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Anche qui, consideriamo H abbastanza grande da poter tenere conto degli altri fattori all'esponente semplicemente 
scegliendo opportunamente i coefficienti, in questo caso pari a j , ^ . 



ovvero anche il raggio d'interazione può essere reso arbitrariamente piccolo. Dunque, la funzione di 
partizione (525) può essere riscritta come: 



i-i 



(534) 



-[0.....Ì1 



dove è arbitrariamente piccolo se \X\ > 1 e U^{l3i), ^^.{/3i) rappresentano dei potenziali a singolo 
corpo; tutto è invariante per traslazioni lontano dai bordi del i?-reticolo. Infine, se Wj^, sono due generiche 



sequenze compatibili possiamo tener conto del prodotto di funzioni di partizione nella ( 534 1 sostituendo 
U^{/3i), ) con i nuovi potenziali 



(535) 

(536) 
(537) 



- ttB 



t/«(A)-l0^ 



XS, (R \ ]nr Zh{f3^,^^+l)Zh{0Jl,Ul2) 



x> 



se \X\ > 2; 



con questa definizione, la ( 534 ) diventa 



(538) Zjy^ e-2:x,|x,>i*x(^,)g-ELo^"(ft)e-i:U*«.(ft). 

^[0....,, 

Abbiamo fatto questa scelta perché, come vedremo nella prossima sezione. 



(539) 



hm sup log — — ^ - 0; 



dunque, scegliendo h abbastanza grande i potenziali ^'x verificheranno le stesse stime che abbiamo 
ricavato per eventualmente modificando in modo opportuno le costanti che vi appaiono. Con ciò, il 
problema è definitivamente ricondotto allo studio di un modello unidimensionale caratterizzato da potenziali 
a più di un corpo con taglia e raggio d'interazione piccoli a piacere, quindi siamo nelle condizioni della 
proposizione |2.2[ in particolare, grazie al corollario |2.1| la pressione del nostro mo dello di Fisher iniziale 



P — lim;_>.oo A(Z) logZw.s è analitica in /3. Per concludere, dimostriamo la proprietà (539). 



4. "Fattorizzazione" di una funzione di partizione in una dimensione 



Per dimostrare la proprietà ( 5391) ci serviremo di un risultato di Ruelle ben noto in teoria dei sistemi 
di spin unidimensionali, |36j , |37 . |8| . Se (ct) è l'energia potenziale per sito di una misura di Gibbs con 
potenziali {<ì>x}, cf. definizione 5.1 introduciamo la sua restrizione A,^{a) su Z"*" {i.e. sugli interi > 0) come 



(540) 



X30.XCZ,+ 



dove $x S B sono i potenziali della mis ura d i Gibbs su Z; diremo che m è una misura di Gibbs su Z"*" 
se le probabilità condizionate {cf. definizione 4.1) verificancj^ ponendo ct' — (o:ó^^(a,oo])^ 2^" = {!z!,i,<L{a,oo])j 



^■^Poiché A è definita in {1, n}y , la somma argomento dell'esponenziale nella 1 541 1 coinvolge al massimo a + 1 termini. 



(541) 



77l(a(,',...,<Va+l-) V 



Proposizione 4.1. Introduciamo /'operatore di trasferimento come 



(542) 



(£*/) (aoai...) = ^ e-^(--«--)/(aaoai...); 



(T = 



se f è una funzione continua nello nello spazio delle sequenze in allora, indicand O COTI 

raggiunto di C^, esiste > tale che le equazioni 



(543) 



ammettono un'unica soluzione h^, m^, con /i$ continua in {!,..., e definita in 
Inoltre, indicando con to+(/) valore di aspettazione di f con questa misura. 



(544) 



lim sup \X^''Clf - m<[.(/)/ì| = 0. 



/c->-oo 

Dimostrazione. Rimandiamo a [36], |8]. 
Consideriamo la funzione /^^ definita nel seguente modo 



□ 



(545) fl^^ia) 



-^^^".^€10 l^-ll.xc{z~u[o,...,,.-l]} *x(ft<T)x ^ _ ^ ^ ^^^^ Compatibili 

altrimenti 



allora, possiamo scrivere la funzione di partizione Z^iPi, fii+i), cf. formula (4501, come, se a;(/3i+i) è 
una sequenza semi-infinita compatibile con /3i+i, 

(546) ^/.(A,A+i) = (>C^/5a) (A+i^(ft+i)). 

Quindi, se cji, lj2 sono due generiche sequenze compatibili. 



(547) 



Zh{(3^^3^+l)Zh{ull,u2) _ 'C„/„,^.(A+i,w(A+i))/:^/^_^^(cj2) 



Z^{f3,,ui,)Z^{ui„/3,+^) /:^/^_^^(a;2)/:S/5,^^(/3.+i,^(ft+i))' 
e grazie alla proprietà ( |544 ) 



(548) 



>Ca/a,ft(/3.+l,^(ft+l))>C^/^.^^(u^2) ^ ^ 



La proprietà (544| è l'analogo del teorema di Perron - Frobenius per l'operatore di trasferimento £; se 



è l'autofunzione corrispondente all'autovalore massimo (l'unico positivo secondo il teorema (4.1 1), possiamo 



decomporre una qualsiasi funzione continua g definita in {1, n}rp nella base delle autofunzioni di come 



(549) g = ciha+g, 

dove Ci è la proiezione di g suha, ovvero per la ( 543 ) Ci — fn^ [g) ; applicando X^'^C'^ alla ( 549 ) otteniamo 

che 



(550) X-''£''^g = m+{g)h^\^l 

se A è il secondo autovalore in modulo piìi grande appartenente allo spettro di C. Quindi, poiché il 




coefficiente m'^{g) non dipende dal punto in cui era inizialmente calcolata g (al contrario di ha), la (548) è 
immediata. 



APPENDICE D 



Analiticità di un valor medio SRB 



Considerando una generica osservabile O analitica in un parametro e e in x, ci chiediamo quale sia il 
comportamento del valor medio SRB nel caso in cui anche i potenziali della misura dipendano analiticamente 
da e. In particolare, se (C')a ^ 1^ derivata di un valor medio calcolato con la misura SRB approssimatsj^ 
ovvero 

A/2-1 

(551) - ° 5')a - (0)a (a.A. ° ^')a) , 

fe=-A/2 

vorremmo dimostrare che per A — )■ oo il limite esiste ed è uguale a 



(552) 



5e(0) , = lim d,{0) 



A— i-c 



La strategia che seguiremo sarà quella di sfruttare le tecniche illustrate nell'appendice [C| Grazie all'in- 
varianza per traslazioni temporali della misura SRB possiamo riscrivere il valor medio di O come 



(553) 



quindi, per la (553 ), indicando con {'l'x} i potenziali a decadimento esponenziale invarianti per traslazioni 
della misura SRB, con un argomento del tutto analogo a quello esposto all'inizio dell'appendice [C] o nella 
dimostrazione del teorema di fluttuazione si ottiene che 



= lim hm ;^a^log?^ 



E, 



ExcA*x(<Ix) 



Ox{<Lx) 



13=0 

T Ti Cx(£x) 



^lim -a.log- 



e i 2 • 2 \ 



/3=0 



(554) 

Nella (554), {Ox{<Lx)} sono i potenziali che si ottengono sviluppando la funzione O, cf. proposizione 
|4.H anch'essi invarianti per traslazioni ed esponenzialmente decrescenti. A questo punto usariamo le tecniche 
esposte nell'appendice [C] (decimazione e cluster expansion) per scrivere la ( 554 1 come 



^Cf. capitolo [s] sezione 



5.2 
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(555) 



1 —>-oo 1 ' ^ 



/3=0 



ed è facile dimostrare che la nuova attività di una configurazione di polimeri F 9^C/3(r)|^^o P^^ essere 
stimata in modo analogo a C(7)- Infatti, dopo la procedura di decimazione avremo che 



(556) 



|CM7)l<e-^«l^le-''''(^\ 



per Lq,k > opportuni, cf. proposizione |2.2[ quindi, poiché (c/. formula (515)) 

(557) cm7)=(e e n(^"'^^"^^^-i))' 

dove $^ è il potenziale analitico in che si ottiene dopo aver decimato il modello, usando una stima di 
Cauchjj^per maggiorare le derivate in f3 di dei potenziali con i loro moduli possiamo trovare]^ C > tale che 



(559) 



dunque il fattore I7I non pone alcun problema, dal momento che può essere controllato con e 2 



Grazie all'invarianza per traslazioni dei potenziali e alla loro analiticità in e, dal corollario 2.1 sappiamo che 
il limite T — >■ 00, ovvero il valor medio SRB di O, è analitico in e. 

Infine, la derivata del valore di aspettazione con la misura SRB approssimata {0)j, può essere scritta 
come 



(560) 



i E v'^(r)9eKMr)l^=o' 



re 



X zi 

2 ' 2 J 



e per le (557), (559), in un modo completamente analogo a quello usato per ottenere la stima (5591, 



abbiamo che, se C è un'opportuna costante positiva. 



(561) 



■^Come abbiamo visto nella dimostrazione del teorema di Anosov per il gatto di Arnold, cf. sezione [g] del capitolo [2] se 
T> = {fi : Im/3 < /3o} con /3o > 0, dal teorema di Cauchy sappiamo che possiamo stimare la derivata in /3 di una funzione 
analitica f{P) come 

(558) dfifiP) < 



/3o 



Il fatt ore C I7I tiene conto anche delle derivate sulla "misura" con la quale è calcolato il "valor medio" nella I557i, cf. 
proposizione 2.2 



quindi, grazie al risultato (561 ) e all'invarianza per traslazioni dei potenziali la (560) converge uniforme- 
mente in e per T — >■ oo, e il limite è analitico. Questo risultato ci permette di scambiare la derivata rispetto 
a e con il limite T ^ oo che definisce la SRB, ovvero 



(562) 

e dunque 



de lim {0)j, ^ lim de (0)-^ , 



(563) 



OO 
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